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GEORGIY  IVANOVICH  PETROV 
Cfor  his  60th  birthday) 


On  May  31 ,  1972,  the  eminent  Soviet  Mechanical  Scientist 
Georgiy  Ivanovich  Petrov  celebrated  his  60th  birthday  and  37 
years  of  scientific,  educational  and  social  activity. 

In  the  years  of  his  youth,  Georgiy  Ivanovich  intimately  en¬ 
twined  his  scientific  interests  with  the  work  of  the  Progressive 
Soviet  School  of  Aerodynamics,  which  was  headed  at  that  time  by 
academician  Sergey  Alekseyevich  Chaplygin.  In  193^>  a  year  be¬ 
fore  graduating  from  the  Mechanical-Mathematical  Department  of 
Moscow  State  University,  Georgiy  Ivanovich  began  to  work  at  TsAGI. 

His  basic  works  of  this  period  were  devoted  to  questions  of  the 
propagation  of  vibrations  in  a  viscous  fluid  and  the  theory  of 
hydrodynamic  stability. 

During  the  years  of  World  War  II,  Georgiy  Ivanovich  partici¬ 
pated  in  work,  vitally  important  for  our  country,  in  perfecting 
aerodynamics  and  combat  characteristics  of  aircraft. 

In  1955,  already  being  a  well-known  scientist,  G.  I.  Petrov 
was  chosen  head  of  the  Department  of  Aeromechanics  and  Gasdynamics 
of  the  Mechanical-Mathematical  faculty  of  Moscow  State  University. 
Since  that  time,  a  significant  part  of  Georgiy  Ivanovich's  scientific 
and  educational  activity  has  taken  place  within  the  walls  of  his 
alma  mater.  In  addition  to  work  in  the  department,  G.  I.  Petrov 
is  the  scientific  supervisor  of  the  Department  of  Aeromechanics 
of  the  Institute  of  Mechanics  of  Moscow  State  University,  formed 
in  I960,  and  of  a  number  of  subunits  of  the  computer  center  of 
Moscow  State  University. 

Georgiy  Ivanovich  was  among  the  first  scientists  who  collec¬ 
tively  and  worthily  evaluated  the  role  of  computer  technology  in 
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the  development  of  science,  and  specifically,  in  theoretical  and 
applied  mechanics.  Already  many  years  ago,  a  seminar  led  by  him 
was  at  work  at  the  computer  center  at  Moscow  State  University. 

The  trends  associated  with  the  use  of  numerical  methods  in  differ¬ 
ent  branches  of  mechanics  are  being  successfully  developed  in  the 
works  of  his  students,  who  are  chiefly  graduates  of  the  Department 
of  Aeromechanics  and  Gasdynamics  of  Moscow  State  University. 

Many  young,  or  already  not  so  young,  students  of  Georgiy 
Ivanovich  consider  the  scientific  research  seminar  which  functions 
within  the  Institute  of  Mechanics  to  be  a  brilliant  school  of 
development.  The  variety  of  themes  of  reports  given  at  the  seminar 
reflects  the  breadth  of  scientific  interests  of  Georgiy  Ivanovich  - 
from  the  flight  mechanics  of  insects  and  aerodynamically  low 
velocities  to  the  hydrodynamics  of  processes  that  occur  on  the  sun 
and  the  machanics  of  formation  of  lunar  craters.  Unchanging  good¬ 
will,  a  burning  interest  in  anything  new,  no  matter  how  small  a 
result  with  respect  to  the  seminar  leader,  serves  as  good  payment 
for  the  risk  for  each  young  person  who  sets  out  on  the  thorny  path 
of  scientific  investigations. 

The  collective  of  the  Institute  of  Mechanics  of  Moscow  State 
University  congratulates  Georgiy  Ivanovich  from  the  heart  and 
wishes  him  strong  health  and  new,  great  creative  successes  in  his 
scientific,  educational,  and  social  activity  for  the  good  of  our 
great  motherland. 
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FOREWORD 


This  publication  of  scientific  transactions  is  the  second 
edition  of  a  collection  of  articles  on  the  aerodynamics  of  high 
velocities.  In  this  edition,  the  results  of  investigations  of 
a  number  of  two-dimensional  and  three-dimensional  problems  of 
exterior  aerodynamics  are  presented,  as  are  the  results  of  cer¬ 
tain  problems  associated  with  the  stability  and  the  excitation 
of  self-excited  vibrations  in  gas  flows. 

The  first  four  articles  examine  axisymmetrical,  supersonic 
streamlining  of  bodies  by  a  nonviscous,  perfect,  and  real  gas. 

The  article  by  M.  G.  Lebedev,  L.  V.  Pchelkina  and  K.  G.  Savinov 
presents  the  results  of  extensive  calculations  made  by  the  method 
of  establishing  supersonic  streamlining  of  bodies  of  different 
shape.  The  next  article  presents  a  differential  system  of  a  high 
order  of  accuracy  developed  by  V.  V.  Yeremin  and  Yu.  M.  Lipnitskiy 
This  system  makes  it  possible  to  use  the  establishment  method  to 
calculate  complex  two-dimensional  flows  without  the  explicit 
separation  not  only  of  shock  waves,  but  also  tangential  discon¬ 
tinuities,  and  examples  of  such  calculations  are  given  for  the 
first  time.  The  work  of  Ye.  G.  Shapiro  Is  devoted  to  the  system¬ 
atic  investigation  of  supersonic  streamlining  of  a  body  by  a  gas 
with  a  high  static  temperature;  the  basic  flow  characteristics  are 
presented  in  the  form  of  universal  relationships  with  certain 
parameters;  the  possibility  of  describing  an  equilibrium  gas  flow 
in  the  examined  case  with  the  aid  of  the  effective  adiabatic  ex¬ 
ponent  is  shown.  V.  P.  Stulov’s  article  evaluates  the  applicabil¬ 
ity  of  the  law  of  binary  similitude  during  streamlining  of  bodies 
and  obtains  a  boundary  equation  of  the  region  of  change  in  parame¬ 
ters  p^  ,  Rq,  where  this  law  is  valid. 

V,  N.  Mirskiy's  work  describes  the  algorithm  of  calculating 
conditions  on  the  incident  and  reflected  shock  waves  in  an  equili 
brium  gas  of  random  composition.  The  work  gives  an  example  of 
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calculation  for  a  mixture  of  1^0  +  Ar. 

The  next  two  articles  are  devoted  to  investigating  super¬ 
sonic,  three-dimensional  streamlining  of  blunt  bodies.  The 
article  by  Yu.  M.  Lipnitsfciy,  Yu.  A.  Mikhaylov  and  K.  G.  Savinov 
describes  the  differential  system  of  calculating  spatial,  super¬ 
sonic  streamlining  of  blunt  bodies  with  the  presence  of  angles  of 
attack  and  sideslip;  the  results  of  investigating  streamlining  of 
triaxial  ellipsoids  are  given.  The  work  of  A.  K.  Burdel'nyy, 

V.  B.  Minostsev  and  V.  P.  Shkadova  examines  nonequilibrium,  super¬ 
sonic,  spatial  streamlining  of  segmentally  conical  bodies. 

The  articles  by  G.  N.  Andreyev  and  V.  G.  Gromov  are  devoted 
to.  supersonic  streamlining  of  bodies  by  a  viscous  gas.  The  first 
of  these  gives  a  detailed  investigation  of  a  three-dimensional 
boundary  layer  on  the  flow  line  of  segmentally  conical  bodies  and 
spherically  blunted  cones.  Finally,  attention  is  given  to  study¬ 
ing  the  effect  of  the  character  of  pressure  distribution  on  the 
appearance  of  separation.  In  the  second  article,  a  solution  to 
the  Navier-Stokes  equations  in  the  vicinity  of  the  critical  line 
is  given  for  a  case  of  supersonic  streamlining  of  a  sphere  by  car¬ 
bon  dioxide  gas.  The  results  of  calculations  for  a  broad  range 
of  Reynolds  numbers  are  given. 

The  article  by  M.G.  Lebedev  and  G.  F.  Telenin  is  devoted  to 
an  Investigation,  suggested  by  the  authors  in  the  first  collection, 
of  the  mechanism  of  excitation  of  self-excited  vibrations  of  super¬ 
sonic  streams.  The  problem  Is  solved  in  the  complete  statement 
for  plane  and  axisyrametrical  streams.  A  comparison  is  made  with 
experimental  data.  The  article  by  S.  Ya.  Gertsenshteyn  and  A.  V. 
Kashko  cites  an  investigation  of  the  stability  of  an  axisymmetrical 
wake  In  a  supersonic  flow.  The  problem  is  s'olved  within  the  frame¬ 
work  of  an  ideal  fluid. 
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I.  PLOWS  OP  AN  IDEAL  GAS 


SOLVING  PROBLEMS  OP  GASDYNAMICS  BY  THE  ESTABLISHMENT 

METHOD 

M.  G.  Lebedev,  L.  V.  Pchelkina,  K.  G.  Savinov 


Beginning  in  1966,  M.  G.  Lebedev  and  K.  G.  Savinov,  at  the 
Institute  of  Mechanics  of  Moscow  State  University,  and  Yu.  N. 

D'yakonov  (deceased),  and  L.  V.  Pchelkina,  at  the  Computer  Center 
of  Moscow  State  University,  conducted  numerical  investigations  of 
various  gas  flows  by  the  establishment  method  [3] •  The  obtained 
results  have  been  partly  reflected  in  the  works  of  authors  [5-12] . 
In  this  article,  a  survey  of  the  basic  obtained  results  of  both 
gasdynamic  and  methodological  characters  is  given.  In  this  case, 
a  great  deal  of  attention  is  given  to  earlier  unpublished  results. 


1.  About  the  establishment  method 

In  1961,  K.  I.  Babenko  and  G«  P.  Voskresenskiy  suggested  a 
finite- difference  method  of  calculating  supersonic,  spatial, 
steady  state  gas  flows  [1].  A  monograph  [2]  contains  a  detailed 
presentation  and  the  foundation  of  the  method.  The  method  has 
been  extensively  used  for  calculating  the  supersonic,  spatial 
streamlining  of  sharp  and  blunt  cones  and  bodies  of  certain  other 
shapes.  The  results  of  the  conducted  investigations  are  in  works 


One  should  take  note  of  the  important  role  which  the  works 
of  Yu.  N.  D'yakonov  [5,6]  played  In  the  development  of  this  method. 
The  methodological  investigations  carried  out  by  Yu.  N.  D'yakonov 
made  it  possible  to  make  a  number  of  important  improvements  in  the 
method  [1]  and  to  apply  it  to  solving  a  number  of  gasdynamic  prob¬ 
lems  for  which  this  method  was  unsuited  in  its  initial  form 


The  experience  of  making  calculations  by  the  finite-difference 
method  [1]  makes  it  possible  to  obtain  numerical  solutions  with 
high  accuracy  to  various  complex  problems  of  gasdynamics,  without 
at  the  same  time  making  extreme  demands  on  the  speed  and  memory 
of  the  computer  used  for  making  the  calculations. 

In  1964,  V.  V.  Rusanov  [3,  21]  suggested  a  modification  of  a 
method  [1]  applicable  to  solving  two-dimensional,  unsteady  equa¬ 
tions  of  gasdynamics.  Inasmuch  as  particular  attention  attaches 
to  steady  state  solutions  which  can  be  obtained  from  unsteady 
equations  when  t  -*  «>  and  from  boundary  conditions  of  the  problem 
which  do  not  change  in  time,  the  method  1 3 1 >  like  other  methods 
of  the  solution  to  unsteady  problems,  is  frequently  called  the 
establishment  method. 

The  establishment  method  examined  by  the  authors  was  presented 
in  detail  in  I3,4,6jj  but  inasmuch  as  we  intend  to  dwell  on  certain 
results  of  a  methodological  character  in  this  article,  it  seems 
useful  to  briefly  describe  the  course  of  solving  the  gasdynamic 
equations  by  this  method.  We  shall  do  this  based  on  the  example 
of  a  problem  of  flow  arund  a  blunt  body. 

Let  CFig.  1)  a  blunt  rotational  body 
be  flowed  past  by  a  supersonic  gas  flow 
running  along  the  body  axis,  of  symmetry. 

Let  the  equation  of  contour  of  the  body 
in  the  polar  coordinate  system  in  the 
meridional  plane  be  * •'&<», t)  (we  assume 
that  the  body  contour  can  change  in  time 
according  to  a  certain  assigned  law),  and  let  r.rio.t)  be  the  equa- 

« 

We  note  that  there  are  also  other  modifications  of  the  examined 
method  today  Cfor  example,  see  the  article  by  Yu.  M.  Lipnitskiy, 

Yu.  Ya.  Mikhaylov  and  K.  G.  Savinov,  in  this  collection;. 


tion  of  the  regressed  shock  wave.  We  shall  insert  the  coordinate 
£  with  the  relationship: 

~  r(9.t)-6(0,t)  (1.1) 


such  that  £  =  0  on  the  body  surface,  and  5  =  1  on  the  shock  wave. 
The  gasdynamic  equations  for  a  perfect  gas  and  the  variables  t, 

£,  9  can  be  written  in  the  form: 


n  .it 
It  *a7S 


*  r 


» 0 


where 


(1.2) 
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pressure. 


—  velocity  components  in  the  direction 
P  —  density, Y  —  adiabatic  exponent 

t  ._*/**•  I(»r/)t-i6 Jit)  t  _  36/iOfi(»F/i0.»6/M) 

*  F*  TTa - 


x  .  9,  p  - 


The  precise  solution  to  problems  of  the  mathematical  state¬ 
ment  of  the  problem  of  steady  state  flow  past  a  blunt  body  with 
a  regressed  shock  wave,  i.e.,  questions  of  the  existence,  singular¬ 
ity  and  continuous  relationship  with  the  boundary  conditions,  poses 
great  difficulties  because  of  the  nonlinearity  of  the  gasdynamic 
equations.  Model  problems  which  conserve  the  basic  aspects  of  the 
studied  problem  are  usually  used  for  the  qualitative  examination 
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of  questions  of  uniqueness  and  the  continuous  relationship  with 
the  boundary  conditions.  One  can  use  plane,  potential  transonic 
currents,  described  in  the  plane  of  a  hodograph  by  an  equation 
of  the  Tricomi  type  CChaplygin  equation),  as  well  as  certain  linear 
equations  in  the  vertical  plane  of  such  models  [  1 9 , 4 ]  . 

Let  characteristics  of  the  second  and  first  families,  which 
intersect  at  a  point  c  CFig.  1),  correspondingly  be  drawn  from 
points  A  of  a  shock  wave  and  B  of  the  contour  of  a  body  in  which 
a  flow  is  supersonic.  Then  the  relationships  in  the  region  of  the 
shock  wave  AD,  the  boundary  condition  on  the  part  of  the  contour 
BE,  and  the  conditions  of  symmetry  on  the  part  of  the  axis  DE  will 
uniquely  determine  the  steady  state  solution  to  the  system  of  gas- 
dynamic  equations  in  the  region  of  the  shock  layer  DACBE  and  en¬ 
sure  a  continuous  relationship  with  the  boundary  conditions.  For 
the  existence  of  solutions  with  limited  derivatives  of  the  gas- 
dynamic  parameters,  it  is  sufficient  that  the  body  contour  not 
have  breaks  in  the  subsonic  flow  region,  and  that  the  distribution 
of  parameters  of  the  advancing  flow  and  the  body  contour  in  the 
supersonic  flow  region  ensure  the  absence  of  boundary  lines.  It 
should  be  noted  that  the  appearance  of  singularities  with  infinite 
derivatives  in  the  isolated  points  on  the  sonic  line  does  not  disrupt 
the  uniqueness  of  the  problem  Cfor  example,  when  the  sonic  line 
leaves  the  point  of  inflection  on  the  surface  of  the  body  or  when 
the  body  is  streamlined  by  a  stream  of  sufficiently  small  cross 
section) . 

We  finally  note  that  positions  closest  to  the  axis  of  all  the 
possible  positions  of  points  A  and  B  CFig.  1)  are  determined  from 
the  condition  such  that  the  characteristics  AC  and  BC  which  originate 
from  these  points  limit  the  minimum  region  of  effect  of  the  super¬ 
sonic  region  of  effect  on  the  subsonic  region. 

We  shall  now  assume  that  the  vector  x  and  the  position  F  of 
the  shock  wave  has  been  totally  assigned  at  a  moment  of  time  t  »  0 
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in  a  region  of  the  shock  layer  o«<*i,o«e*».  .  Then,  in  order  to 
determine  the  gas  parameters  in  this  region  when  t  >0,  it  is 
necessary  to  solve  a  mixed  Cauchy  problem  for  system  (2)  with  ini¬ 
tial  data  when  t  *  0  and  boundary  conditions  on  the  contour  of  the 
body,  the  shock  wave,  and  the  axis  of  symmetry.  In  order  that  the 
unsteady  problem  be  correctly  stated,  it  is  necessary  to  draw  a 

boundary  ray  9=0.  such  that  the  velocity  component  v  normal  to 

K  * 

the  ray  exceeds  the  local  velocity  of  sound  c  everywhere  on  it. 


We  insert  a  grid  with  nodes 

*  »  0„0««  .1,  £7*  0  i  t), 

in  the  investigated  region  v"'**-6/) 
We  designate  (Pig.  2) 


8  •  8  *10  4  ■(  < ( 
t  t->  •  '<  «  « 


t  4  •  4  *  »4 

•  • 

•  o 

0  9  • 

<-■'  .  V*.  • 


m- 0,t . n, 

e,*.f . l  . 


We  replace  derivatives  in  the  points  r».  »/*,«•♦//> ,t) 

with  the 

difference 

relationships  according  to  the  following 

formulas : 

9*  ( 

\  *  1  /  *  mot  m  . 

**m,e  ~*m,l  /* 

(1.3a) 

~ui f 

/it  t  r ,  *•*  •»• .  •  » 

l.'  f  ■.7^*-'.'  *  t  *  *mf  )]  , 

Cl. 3b) 

Here,  too,  we  note  that,  as  the  experience  of  calculations  shows, 
in  certain  cases  one  can  obtain  a  numerical  solution  to  the  problem 
without  significant  errors,  and  without  satisfying  the  conditions 
v  >c,  when  9=  9k,  and  even  without  drawing  the  boundary  ray  into 
the  region  of  supersonic  velocities.  It  is  usually  sufficient  for 
this  purpose  that  velocity  on  the  boundary  ray  be  a  magnitude  of 
precisely  the  same  order  as  the  speed  of  sound.  Certainly,  the 
question  of  the  permissibility  of  drawing  the  boundary  ray  in  the 
subsonic  region  requires  an  examination  in  each  specific  case.  One 
of  such  cases  is  presented  in  Sec.  5. 
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Cl. 3c) 


The  value  of  X  in  the  point  , ..  w  .  ~‘<n  .  e)  is  defined  as  the 
mean  arithmetical  value  of  X  in  the  four  nearby  nodes 


(*L,.e 


/  T  "**  1 

>  —  fx  ♦  x 

^  i  «•*»/// 


/ 

»  *  *ft  t  C  J 


(1.3d) 


In  formulas  (1.3) 


X 


•*»  ♦  t/f ,  c 


X 


) 


The  parameters  of  the  difference  system  a,  3  satisfy  condition 
«*<7>  j  in  most  of  the  calculations  that  were  made, 

it  was  assumed  that  -  -  o.aa,  a  -  0.4a. 


The  second  term  in  (1.3a)  is  artificial  viscosity  of  the  dif¬ 
ference  system;  the  coefficient  of  artificial  viscosity  is  « » 0  . 
In  certain  cases,  it  proves  that  one  cannot  obtain  a  numerical 
solution  to  the  problem  (see  Sec.  4)  without  inserting  artificial 
viscosity. 


By  substituting  Cl, 3)  and  Cl. 2)  for  the  point  *-.■/.  •  9t) 

we  proceed  to  the  system  of  difference  equations: 


m»t,t 


^ m*  tft ,  t  ^  m9t 


• ».  <ft 

c 


(1.4) 


where 


a  •  f  *  £«<  r,  4  f  I  *  £  - 1  •»  4  , 


c-2*tr»(£  '*/*,*)* m  * 


*  „  00  f  ,  *  *  . 

•  t  bU(X,  -X,  )*/(x.  -x,  )• 


(1.5) 
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,  t  *t 

*,  •  -Vi*  •  ,  E  —  unit  matrix. 

During  the  subsequent  presentation,  as  in  (.1.5),  the  obvious  indices 
shall  be  omitted. 


Thus,  the  boundary  value  problem  for  the  system  (1.2)  pertains 
to  a  boundary  value  problem  for  the  system  of  difference  equa¬ 
tions  (1.4)  with  initial  data  when  t  =  0  and  boundary  conditions 
on  the  surface  of  the  body  and  the  shock  wave. 

It  is  evident  from  the  structure  of  the  system  (1.3)  that  it 
is  situated  on  (L  +  1)  independent  subsystems.  In  this  connection, 
it  is  sufficient  to  examine  the  solution  to  difference  equations 
on  the  ray  9  =  const.  Inasmuch  as  the  coefficients  of  the  equation 
depend  upon  values  of  the  functions  for  ”nw  and  "n+1"  layers,  in 
the  process  of  solution  one  uses  a  method  of  iteration  according 
to  the  following  scheme:  At  the  first  iteration,  during  the  de¬ 
termination  of  matrix  elements  a,  b  and  column  c,  the  values  of 
all  parameters  are  assumed  to  equal  their  values  at  the  previous 
step.  The  obtained  vector  Xn+1  is  used  for  the  subsequent  itera¬ 
tion,  etc. 


The  trial  and  error  method  Is  employed  for  solving  the  system 
of  idifference  equations  of  each  ray.  In  addition  to  (1.4),  we 
shall  examine  the  system 


ix  •  e 

m*  t  ** 


Cl. 6) 


which  is  obtained  from  Cl. 2)  by  means  of  discarding  the  last  equa 
tion,  and  (the  stroke  signifies  transposition).  The 

boundary  condition  on  the  body  surface  is  presented  in  the  form 


*»  r 


m  n+t 

X 


(1.7) 


where 


.V •  (anj3 e',)',{o, 

’./’o',  w)'}'*  a  ■  u*/*t ) 
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One  can  obtain  the  recurrent  formulas  for  calculating  ^ 

and  (»•'.* in  the  following  relationship,  with  the  aid  of 

(1.6)  and  (1.7): 


'fm 


These  formulas  have  the  form 


Cl. 7a) 


*m  ( 3  )m.<n 


Cl. 7b) 


The  conditions  of  conservation  on  the  shock  wave  (.for 
)  can  5e  written  in  the  form: 

Cl. 8) 

f(v, •  7*  'f7'vi-jrr • 

V*  7  f 

■  "ll  *J?7  P~N  • 

Here  the  indices  »*t  ■,  t  are  omitted  everywhere;  values  with  the 
index  "°°"  are  determined  by  conditions  in  the  free-stream  flow 
and  are  considered  known.  *  .  •  v*  —  respectively,  velocity  of 

the  shock  wave  and  the  components  of  gas  velocity  that  are  normal 
and  tangential  to  the  shock  wave;  they  are  determined  through  the 
position  of  the  shock  wave  and  the  velocity  components  u,  v  accord¬ 
ing  to  the  formulas: 

Cl. 9) 
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If  the  values  F,  F0  are  considered  to  be  known  according 
to  the  previous  iteration  in  equations  Cl. 8),  then  (1.8)  is  a 
system  of  four  equations  relative  to  five  unknowns 
The  trial  and  error  relationship 

(1.10) 

when  m  *  M  closes  the  system  of  equations  on  the  shock  wave. 

One  can  solve  the  system  Cl. 8),  (1.10),  for  example,  by  the 

n+2 

Newton  method.  Having  thereby  determined  the  vector  ,  one 

■n  -4*1 

can  successively  find  x“  for  m  ■  M  -  1,  M  -  2,....  0,  l.e., 
carry  out  an  inverse  trial  and  error  run.  For  this  purpose,  one 
should  use  the  system  consisting  of  four  equations  obtained  by 
means  of  a  linear  combination  (.1.4),  (.1.7  a).  As  a  rule,  a  sy¬ 
stem  consisting  of  the  first,  second  and  fburth  equations  of  (1.4) 
and  the  trial  and  error  relationship  Cl. 7  a)  was  used. 


Finally,  the  new  position  of  the  shock  wave  can  be  calculated 
according  to  the  formula: 


>$(r'  * f"*' ) 

i  *  «  <  / 


Cl. 11) 


Moving  from  layer  to  layer  and  determining  the  gas  parameters 
for  t  •oa*  •  ••  .  L  for  each  of  them,  we  find  the  vector  X  in  the  en¬ 
tire  region  of  interest  to  us  t  »  a  ,  o* k * t  ,  a* 


The  conditions  of  stability  of  the  trial  and  error  method 
impose  limitations  on  the  use  of  the  method  in  the  described  form. 
One  can  derive  the  conditions  of  stability  of  the  trial  and  error 
method  and  the  establishment  method  similar  to  the  way  this  was 
done  in  [2],  They  have  the  form: 


4.  » lt  »»»*<.- 

u;  ♦(*•'<.)* 


(1.12) 

(1.13) 
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♦  c  >  o 


Conditions  Cl. 12)  and  (1.13)  can  be  rewritten  in  a  form  which 
makes  the  kinematic  meaning  explicit.  If  v  is  the  gas  velocity 
component  normal  to  the  line  £  =  const  (equation  1.1)),  and 
is  the  rate  of  displacement  of  this  line,  then  we  have  the  follow¬ 
ing  from  (1.12)  and  (1.13): 


vn*** 
VH  *VC 


(1.12a) 

(1.13a) 


In  the  case  of  a  steady  state  flow  *  0 ;  then  it  follows 
from  Cl. 12a)  that  angle  u  between  the  flow  line  and  the  line  £  * 
const  should  be  positive,  as  in  Pig.  3a;  it  follows  from  (1.13a) 
that  the  normal  velocity  V  in  this  case  should  be  less  than  the 
speed  of  sound  with  respect  to  its  absolute  value. 


Pig.  3 

0  -  line  £  *  0,  1  -  line  £  =  1 ,  2  -  intermediate  line 
£  *  const,  3  -  flow  line 

The  condition  of  stability  (1.13  a)  was  always  satisfied  for 
all  of  the  flows  we  calculated.  With  respect  to  condition  (,1.12a) , 
i  in  a  number  of  cases  it  was  violated.  For  example,  in  the 
case  of  flow  past  blunt  bodies,  it  follows  from  Cl.ll  that  the 
line  £  *  const,  drawn  near  the  contour  of  the  streamlined  body, 
will  trace  the  form  of  the  latter;  this  cannot  be  stated  about  the 
flow  line  drawn  near  the  contour,  if  the  shape  of  the  contour  is 
insufficiently  smooth.  In  these  cases;  condition  Cl. 12a)  will  be 
violated;  in  the  steady  state  flow  in  this  case,  there  will  be 
id  <  0  at  certain  points,  as  in  Pig.  3b, c. 


The  situation  depicted  in  Pig.  36  existed  when  making  the 
calculations  described  in  Secs.  2  and  U .  However,  it  proved 
that  If  condition  (1.12a)  is  violated  in  a  small  number  of  grid 
nodes,  then  this  does  not  have  an  effect  on  the  convergence  of 
the  process  of  establishment  and  the  accuracy  of  numerical  solu¬ 
tions.  However,  in  those  cases  when  nonsatisfaction  of  condition 
Cl. 13)  leads  to  the  impossibility  of  obtaining  a  solution,  one 
can  use  the  system  with  recalculation  developed  in  151  or  the 
system  of  two  directions,  which  ensure  stability  of  the  inverse 
trial  and  error  method. 

The  described  algorithm  of  the  numerical  solution  to  gas- 
dynamic  equations  was  programmed  on  the  computer.  The  program, 
devised  at  the  Institute  of  Mechanics  (subsequently,  the  IM  pro¬ 
gram)  was  written  for  a  differential  grid  with  11  x  21  points 
(11  in  the  £  direction,  21  in  the  9  direction);  the  program 
written  at  the  Computer  Center  (subsequently,  the  VTs  program) 
was  written  for  a  grid  with  9  x  17  points.  The  algorithms  realized 
by  the  IM  and  VTs  programs  differ  from  each  other  in  certain  de¬ 
tails  as  well  as  from  that  described  above.  For  example,  the  VTs 
program  selected  values  and  <•«'*:».  as  dependent  variables 

In  the  gasdynamic  equations,  while  the  IM  program  directly  used 
p  and  p.  The  singularity  of  the  equations  of  gasdynamics  on  the 
axis  of  symmetry  was  revealed  differently;  derivatives  according 
to  9  on  the  boundary  ray  9  *  9  etc.,  were  approximated  accord¬ 
ing  to  different  nonsymmetrical  schemes. 

We  note  that  the  VTs  program  made  it  possible  to  calculate 
the  gas  flow  with  the  presence  of  equilibrium  physico-chemical 
transformations,  while  the  IM  program  was  only  developed  for  cal¬ 
culating  flows  of  a  perfect  gas. 

A  check  of  the  satisfaction  of  certain  finite  relationships 
which  identically  follow  from  the  system  of  gasdynamic  equations 
was  provided  for  controlling  accuracy  of  the  obtained  numerical 
solutions  in  the  written  programs.  But  these  finite  relationships 


do  not  follow  from  the  approximating  system  of  difference  equa¬ 
tions,  namely  the  Bernoulli  equation,  conditions  of  constancy  of 
the  entropy  function  .*■/ :t  along  the  streamline  (specifically, 
along  the  streamlined  surface),  as  well  as  conditions  of  flow  con¬ 
servation  in  the  investigated  region  of  the  shock  layer.  The  accu¬ 
racy  of  the  numerical  solution  to  the  specific  problems  we  have 
examined  will  be  stated  in  the  appropriate  sections. 


Section  2.  Plow  past  a  sphere  and  cylinder  with 
a  flat  leading  part 

The  calculations  of  flow  past  a  sphere  were  made  with  num¬ 
bers  of  the  free-stream  flow  i  "'«*  •x>  .  The  numerical  solutions 

for  the  sphere  are  distinguished  by  high  accuracy.  For  example, 
when  a  4  in  a  point  of  the  outer  ray  lying  on  the  surface  of 
the  body  is  =?.  ,  the  Bernoulli  equation  is  satisfied  with 

an  accuracy  of  •  '.cr\  ,  and  the  condition  of  entropy  constancy 
is  satisfied  with  an  accuracy  of  » '•  ■>  o.v. ^  (with  calculations  made 

according  to  the  IM  program) .  One  notes  good  correspondence  of 
the  obtained  results  with  the  data  of  other  numerical  methods, 
specifically,  G.  P.  Telenin  II 8 ] . 

One  should,  however,  make  one  remark  on  the  subject  of  compar¬ 
ing  the  numerical  solutions  obtained  by  the  two  methods  with  small 
.  The  results  of  calculations  of  flow  past  a  sphere  when  1^J0<  2 
made  by  the  G.  P.  Telenin  method  were  presented  in  1161.  It  was 
concluded  on  the  basis  of  analyzing  the  numerical  data  when  M*  * 

1.2  in  II 6 ]  that  the  angled  between  the  shock  wave  and  the  sonic 
line  is  approximately  90°,  and  that,  consequently,  when  1^  $  1.2 
the  characteristic  bounding  the  minimum  region  of  effect  in  the 
3hock  layer  proceeds  from  the  sonic  point  on  the  shock  wave.  The 
positions  of  the  sonic  points  in  the  vicinity  of  the  shock  wave 
when  M,,,  *  1.2,  calculated  by  the  establishment  method,  coincide 
with  the  corresponding  positions  in  the  calculation  made  according 
to  the  G.  F.  Telenin  method  with  an  accuracy  of  several  percentages. 
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However,  differences  in  the  size  of  the  angle <S  proved  to  be 
greater,  inasmuch  as  the  latter  is  obtained  during  numerical  dif¬ 
ferentiation  of  equations  of  the  sonic  line  and  the  shock  wave. 

One  obtains  5  <90°  in  the  calculation  made  by  the  establishment 
method.  We  note  that  V.  V.  Rusanov  recently  demonstrated  analytic¬ 
ally  that  during  the  axisymmetrical  flow  past  bodies  when  Mot>  1, 
angle  6  is  always  less  than  90°  (with  certain  entirely  natural 
assumptions  made  about  the  shape  of  the  regressed  shock  wave). 

It  should  be  noted  that  differences  in  the  results  of  cal¬ 
culating  flow  past  a  sphere  when  *  1.2  by  the  two  methods  exist 
only  in  a  small  region  of  the  shock  layer  near  the  shock  wave  and 
at  a  great  distance  from  the  axis.  Outside  this  region,  when 

»  1.2,  as  well  as  throughout  the  entire  shock  layer  when  Mao>  1.2, 
the  numerical  solutions  obtained  by  the  establishment  and  G.  P. 
Telenin  methods  coincide  with  good  accuracy. 

Detailed  information  is  given  in  [6]  about  the  results  of 
calculations  of  flow  past  a  sphere  by  equilibrium  air  when  2*  m„«  ao. 
made  by  the  establishment  method.  We  note  here  that  the  equilibrium 
physico-chemical  transformations  have  a  weak  effect  on  the  distribu¬ 
tion  of  pressure  in  the  shock  layer,  in  comparison  with  the  case  of 
the  perfect  gas.  A  different  picture  is  observed  with  respect  to 
the  other  gas dynamic  values,  primarily  density.  With  hypersonic 
velocities  in  the  free-streara  flow,  the  magnitude  of  density  in  the 
shock  layer  can  be  more  than  two  times  greater  than  the  correspond¬ 
ing  density  value  during  flow  by  a  perfect  gas.  The  equili¬ 

brium  dissociation  and  Ionization  have  an  effect  as  well  on  the 
qualitative  picture  of  density  distribution  in  the  shock  layer. 

Thus,  when  a  perfect  gas  with  MB  *  00  flows  past  bodies,  density 
Pw  directly  behind  the  shock  wave  is  constant  and  is  ,  r  ♦  »)/</• >' 

The  distributions  of  density  through  the  shock  wave  during  the  flow 
of  equilibrium  air  around  a  sphere  depending  on  the  angles  with 
hypersonic  are  plotted  in 'Pig.  4.  The  characteristic  points  of 
these  curves  (the  extremum,  the  inflection)  correspond  to  the  cessa¬ 
tion  or  initiation  of  any  particular  equilibrium  processes  with  the 
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change  in  temperature  and  pressure  behind  the  shock  wave  (both  of 
these  values  diminish  monotonically  with  the  increase  in  angle  8 ) . 
Thus,  the  inflection  points  on  curves  for  *  20  and  25  correspond 
to  conditions  under  which  atmospheric  oxygen  is  practically  totally 
dissociated,  and  the  dissociation  of  nitrogen  has  still  not  started. 
The  maximum  point  on  the  curve  for  Moo  =  SO  when  9  “  90°  corresponds 
to  conditions  under  which  both  oxygen  and  nitrogen  are  totally 
dissociated,  but  the  ionization  of  the  atoms  has  still  not  begun; 
the  other  extremums  on  this  curve  are  due  to  the  occurrence  of 
ionization  of  nitrogen  and  oxygen  atoms. 


We  note  that  the  calculations  of  flow  past  a  sphere,  whose  re¬ 
sults  are  given  in  Fig.  4,  were  made  by  the  establishment  method 
according  to  the  VTs  program  up  to  8  =  66° ,  while  the  purely  super¬ 
sonic  flow  when  9  >  66°  was  calculated  by  the  method  of  grids  [1 ] . 


Calculations  of  flow  around  cylinders 
with  a  flat  leading  part  and  circular  edges 
were  made  according  to  the  IM  program 
(11  x  21  points)  with  ratios  R/d  *  0.4, 

0.3,  0.2  and  0,1  OR  -  rounding  radius,  d  - 
one-half  the  middle  diameter) .  The  geo¬ 
metrical  picture  of  their  flow  when  M^* 

3  is  shown  in  Fig.  5.  The  pole  of  the 
polar  coordinate  system  Or,  8)  was  located 
a  distance  d  from  the  critical  point.  When 
R/d  *  0.2  and  the  position  of  the  outer  ray 


*  0.835  at  a  point  on  the  surface  of  the 
body  5  »  0,  8  *  9^,  the  accuracy  of  satisfying  the  Bernoulli  integral 
is  •s-n.and  the  conditions  of  entropy  constancy  4*«a%,  .  in  this 


case,  the  normal  velocity  component  v  to  the  ray  is  greater  than 
the  local  speed  of  sound  c  throughout  the  entire  shock  layer  on 
the  ray  9  *  9^.  Calculation  was  also  made  with  the  outer  ray  pass¬ 
ing  through  a  point  c  of  conjugation  of  rounding  with  the  cylin¬ 
der  (8k  ■  0.8955).  In  this  calculation,  when  4*«u%, 

However,  in  this  case,  pressure  on  the  body  surface  is  calculated 
sufficiently  accurately,  while  the  high  magnitude  of  errors 


«  ao  v 


! 


and  •-»  is  chiefly  due  to  inaccuracy  in  deter¬ 
mining  density  and  velocity.  This  is  con¬ 
firmed  by  comparison  with  the  calculation  of 
supersonic  flow  on  the  curve  made  by  the 
method  [13,  which  was  made  by  V.  I.  Uskov 
from  our  initial  data  when  9  =  0.835. 

When  R/d  *  0.1,  the  accuracy  of  solution 
in  a  point  £  *  0,  9  *  9  ^  is  not  high:  when 
0k  =  0.8125  »s  •2%£a3  •  7  \  when  9^  =  0.8313 
(the  outer  ray  passes  through  c),  “  10  “  33  *•  .  We  note  that 

even  when  9fc  =  0.8313  in  certain  points  on  the  ray  0  *  0^  compris¬ 
ing  velocities  v<c.  One  can  obviously  obtain  a  solution  with 
v  >  c  when  9  »  9  with  significantly  higher  accuracy  by  shifting 
the  poles  of  the  coordinate  system  closer  to  the  critical  point 
and  by  drawing  the  outer  ray  through  the  point  on  the  body  where 
velocity  slightly  exceeds  sonic  velocity  (we  note  that  with  small 
R/d,  the  sonic  point  on  the  body  is  very  close  to  point  B  of  con¬ 
jugation  of  the  flat  leading  part  with  the  rounding). 


We  also  note  that  the  accuracy  of  the  numerical  solution  in 
the  point  Z  »  0,  6  *  9^  is  the  worst;  outside  its  small  vicinity 
Ca  radius  on  the  order  of  2-3  steps  of  the  differential  grid), 
error  was  no  more  than  1%  in  all  calculated  variations. 


The  distributions  of  pressure  p  re¬ 
lated  to  pressure  in  the  critical  point  pQ , 
along  the  surface  of  the  examined  bodies, 
are  plotted  in  Pig.  6  depending  on  the 
ordinate  y.  The  distribution  p/Pg  along 
the  surface  of  the  cylinder  with  a  flat 
leading  part  and  a  point  of  inflection 
(R/d  »  0)  calculated  by  the  G.  F.  Telenin  method  with  the  identi¬ 
fication  of  a  singularity  in  the  point  of  inflection  with  the  aid 
of  the  asymptotic  solution  (223  is  also  given  here.  The  points 
in  Pig.  6  plot  Yu.  Ya.  Karpeyskiy’s  experimental  data,  which 


coincide  well  with  the  calculation  re¬ 
sults  . 


We  note  that  such  values  that 
characterize  the  flow  as  a  whole  as 
the  regression  of  the  shock  wave  on 
the  axis  of  symmetry,  the  velocity 
gradient  in  the  critical  point,  and 
certain  others  change  with  the  change 
in  the  rounding  radius  from  R/d  =  1 
(sphere)  to  R/d  *  0  (disc)  approxi¬ 
mately  linearly.  At  precisely  the 
same  time,  the  local  flow  properties 
on  the  rounded  edge  can  change  ex¬ 
tremely  sharply,  particularly  with 
small  R/d.  Figure  7  plots  the  p/pQ  dis¬ 
tributions  along  the  rounded  surface  of  various  radii  depending 
on  the  angle  of  turn  v  of  the  flow  on  the  rounded  edge  (v  =  0  in 
the  sonic  point).  When  R/d  »  Q  in  the  inflection  point  (more 
accurately,  in  its  infinitely  small  vicinity),  a  Prandtl-Mayer 
flow  takes  place  which  begins  with  M  ■  1  122] .  As  is  evident  from 
the  data  of  Fig.  7,  the  supersonic  flow  on  the  surface  of  a  sphere 
(R/d  *  1)  is  similar  to  a  Prandtl-Meyer  flow.  However,  with  a  de¬ 
crease  in  the  size  of  rounding,  the  difference  between  flow  on  the 
rounded  edge  and  Prandtl-Mayer  flow  increases,  and  only  with  R/d  ^ 

Q.2  does  flow  on  the  rounded  edge  begin  to  tend  toward  the  limiting 
case . 

In  conclusion,  we  shall  illustrate  the  effect  of  the  Moo-number  on 
flow  around  a  cylinder  with  a  flat  leading  part  when  R/d  *  0.2 
(Fig.  8). 

Section  3.  Flow  around  blunt  bodies 
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Most  of  the  data  available  in  the  literature  on  flow  around 
blunt  bodies  relates  to  geometrical  shapes  in  which  curvature 


of  the  surface  is  either  constant  (a  sphere) 
or  increases  with  distance  from  the  critical 
point  (for  example,  ellipsoids  with  a  ratio 
of  the  vertical  axis  to  the  horizontal  axis 
of  5  >1,  cylinders  with  a  flat  nose  part, 
segmental  bodies,  etc.).  This  is  because 
it  is  precisely  such  shapes  which  are  favor¬ 
able  from  the  viewpoint  of  reducing  aero¬ 
dynamic  heating  of  bodies  streamlined  by  a 
supersonic  gas  flow.  Together  with  this, 
the  flow  around  bodies  with  reducing  (be¬ 
ginning  from  the  critical  point)  surface 
curvature  can  also  be  of  known  gasdynamic 
interest . 

Figure  9  gives  the  geometrical  pictures 
of  flow  past  a  rotational  body  when  =  2  and 
6.  The  contour  of  the  rotational  body  consists 
of  two  conjugated  arcs  of  circles  .  The 

radius  R^  of  the  circles  which  intersects  Fig.  9 

the  axis  of  symmetry  ("of  bluntness")'  is  less 

than  the  radius  of  the  second  circles  ("lateral  surface"); 

the  ratio  R£:R^  =  8.5.  The  angle  comprised  by  the  axis  of 
symmetry  and  the  line  which  connects  the  conjugation  point  with 
the  center  of  bluntness  is  greater  than  the  angle  6*  of  the  sonic 
point  on  the  surface  of  a  sphere  streamlined  by  a  flow  with  =  6. 
Therefore,  when  =  6,  the  lateral  surface  has  no  effect  on  flow 
around  the  blunt  part;  the  region  of  the  shock  wave  ahead  of  the 
blunt  point,  the  sonic  line  in  the  shock  layer  coincide  with  the 
corresponding  lines  during  flow  around  a  sphere.  The  pressure 
distributions  on  the  surface  of  a  sphere  with  a  radius  (dash 
line)  and  the  examined  body  are  plotted  in  Fig.  10.  It  is  evident 
from  this  graph  how  acceleration  of  the  gas  flow  decreases  on  the 
surface  of  a  complex  body  behind  the  point  of  conjugation. 

When  M„  *  2,  the  angle  9  <  *. ;  in  this  case,  flow  around  a 
blunt  body  is  affected  by  the  shape  of  the  lateral  surface.  As 


Pig.  10  Pig.  11  Fig.  12 


is  evident  from  Pig.  10,  pressure  differences  on  the  sphere  and 
on  blunting  of  a  complex  body  only  become  noticeable  in  the  vicin¬ 
ity  of  the  conjugate  point.  The  gas  slows  in  this  region  of  the 
body;  at  a  certain  distance  from  the  ponjugate  point,  the  accelera¬ 
tion  of  the  gas  again  begins,  but  occurs  significantly  more  slowly 
than  against  the  blunt  part  of  the  body.  In  accordance  with  this, 
the  sonic  line  is  displaced  significantly  higher  during  flow  around 
a  complex  body  than  in  the  case  of  the  sphere. 

Flow  around  a  blunt  cone  of  finite  length  with  a  large  expan¬ 
sion  angle  when  »  6  was  also  calculated  by  the  establishment 
method  according  to  the  VTs  program.  The  radius  of  spherical  blunt¬ 
ing  was  taken  to  be  R/d  *  0,45,  while  the  angle  of  inclination  of 
the  conical  surface  to  the  axis  of  symmetry  was  @c  *  51  .  The 
examined  body  had  a  slight  rounding  in  the  region  of  the  middle  in 
the  form  of  a  circumferential  arc  with  a  radius  r/d  *  0.2,  conju¬ 
gated  with  a  conical  surface.  We  note  that  the  angle  0^  *  51  is 
less  than  the  maximum  angle  for  flow  around  a  sharp  cone  by  a  flow 
with  *  6,  although  it  was  similar  to  it.  We  note  further  that 
the  angular  coordinate  of  the  sonic  point  on  a  sphere  streamlined 
by  a  flow  with  8,  e.>  x>’  -  or  -  ;  thus,  the  lateral  surface 

affects  flow  around  a  blunt  body. 

Figures  11  and  12  show  the  geometrical  picture  of  flow  around 
a  cone  and  the  distribution  of  pressure  on  the  surface.  We  note 
the  presence  of  a  point  of  inflection  on  the  shock  wave  and  the 
fact  that  the  sonic  line,  running  practically  from  the  conjugate 
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point  of  the  conical  surface  with  rounding  nearly  "spreads  out" 
along  the  body  over  its  entire  course.  The  distribution  of  pres¬ 
sure  in  the  vicinity  of  conjugation  of  the  cone  with  the  spherical 
blunting  is  "spoon  shaped";  similar  behavior  of  pressure  in  this 
region  is  also  characteristic  for  the  flow  around  cones  with  small 
expansion  angles,  when  flow  in  this  region  is  supersonic  16]. 
Pressure  on  the  surface  of  a  cone  with  0 , 55  _<  y  _<  0 . 75  is  nearly 
constant  and  is  similar  to  pressure  on  the  surface  of  a  sharp  cone 
with  precisely  the  same  angle  9^ .  However,  other  gasdynamic  parame¬ 
ters  significantly  differ  from  the  corresponding  values  on  the  sur¬ 
face  of  the  sharp  cone.  Thus,  when  a  sharp  cone  is  streamlined, 
the  local  Mach  number  in  the  examined  case  is  M  *  1.25;  in  the 
case  of  a  blunt  cone,  a  subsonic  value  with  M  *  0.7  corresponds  to 
precisely  the  same  surface  pressure  value. 

The  examples  examined  above  show  that  the  presence  of  regions 
where  a  flow  of  gas  with  negative  acceleration  occurs  is  possible 
on  the  surface  of  bodies  whose  curvature  increases,  beginning  with 
the  critical  point.  This  effect  appears  significantly  more  sharply 
when  there  is  flow  around  bodies  whose  curvature  changes  sign;  this 
question  is  examined  in  the  next  section. 

Section  *J.  Plow  around  convexo-concave  bodies 

We  shall  examine  a  family  of  bodies  with  the  equation  of  the 
contour  generatrix  Cin  the  polar  system  of  coordinates  in  the 
meridional  plane) 


When  m  *  0,  a  body  with  a  contour  (3.1)  is  a  sphere  with  a  unit 
radius;  when  m  >  0,  the  examined  body  has  a  depression  whose  maxi¬ 
mum  depth  is  achieved  when  9  =  45°  and  is  determined  by  the  parame¬ 
ter  m  .  Parameter  K  characterizes  curvature  of  the  depression. 

The  establishment  method  was  used  to  make  calculations  of 
flow  around  bodies  (3.1)  by  a  perfect  gas  (y  *  1.4)  and  equilibrium 
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Pig.  13 

1 .  Sphere  /m  “  °  » a. «  -  oxx. 


Fig.  14 

m  =  0 . 1 ,  k  =  8 


air  in  the  range  of  change  of  parameters  which  determine  the  body 
shape:  0  _<  m_<  0.2,  2  _<  k  _<  8 .  Calculations  were  made  according 
to  the  IM  program  for  3  _<  M,,,  <  6  with  position  of  the  outer  ray 
■  90°,  and  according  to  the  VTs  program  for  a«M*ao  when  6,  * 

O  O  **  iC 

90  and  70  .  When  *  6 ,  a  comparison  of  results  obtained  accord¬ 
ing  to  the  two  different  programs  was  made;  these  results  coincided. 


Figure  13  compares  the  geometrical  pictures  of  flow  past  a 
sphere  and  a  body  with  a  slight  depression  (m  *  0.0  5)  when  M,*  *  3. 

As  is  evident,  even  a  slight  change  in  the  shape  of  the  body  causes 
a  noticeable  change  in  the  picture  of  flow  in  the  shock  layer.  Gas 
accelerates  more  rapidly  on  the  lower  part  of  a  body  with  a  depres¬ 
sion,  where  the  inclination  of  the  contour  to  the  axis  decreases 
more  rapidly  than  in  a  sphere  and  the  gas  accelerates  more  rapidly  than  on  a  sphere. 
This  is  due  to  lowering  of  the  sonic  line  in  the  shock  layer  ahead 
of  the  body  with  a  depression,  in  comparison  with  the  sonic  line 
position  in  the  sphere.  The  change  in  the  sign  of  surface  curva¬ 
ture  on  the  body  depression  (4.1)  leads  to  gas  deceleration  in  this 
region  and  to  the  appearance  of  a  local  zone  of  subsonic  flow  which 
"lies"  in  the  depression,  with  a  repeated  change  in  the  sign  of 
curvature  of  the  surface,  the  velocity  of  the  gas  again  reaches 
supersonic  values  and  subsequent  gas  expansion  occurs  approximately 
as  on  a  sphere.  The  described  flow  characteristics  due  to  the 
double  change  in  the  sign  of  surface  curvature,  extend  through  the 
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shock  layer  and  lead  to  the  appearance  of  two  inflection  points 
on  the  regressed  shock  wave,  which  are  distinguishable  in  Fig. 
13.  We  also  note  that  the  thickness  of  the  shock  layer  during 
the  transition  from  the  sphere  to  the  body  with  the  depression 
decreases . 


Figures  14  and  15  show  more  detailed 
pictures  of  flow  in  bodies  with  depression 
m  =  0.1  and  m  =  0.2  (when  =  3)-  Shock 
waves,  sonic  lines  and  (dash  line)  stream¬ 
lines  are  plotted  here.  The  deepening  of 
the  depression  leads  (Figs.  13-15)  to  a 
certain  decrease  in  the  thickness  of  the 
shock  layer  in  the  region  of  subsonic  and 
transonic  flow  and  to  a  sharper  curvature 
of  the  shock  wave  (the  inflection  points 
on  the  shock  wave  are  the  more  clearly  pronounced,  the  deeper  the 
depression),  and  to  an  increase  in  the  dimensions  of  the  local 
zone  of  subsonic  flow  in  the  depression  both  in  absolute  units , 
and  relative  to  the  thickness  of  the  shock  layer.  When  m  *  0.2, 
the  inflection  of  the  shock  wave  is  so  strong  that  a  local  sub¬ 
sonic  zone  forms  on  the  wave  which  connects  with  the  region  of  sub 
sonic  flow  lying  in  the  depression  (.Fig.  15). 

Deceleration  of  the  gas  to  subsonic  velocities  during  flow 
around  the  depression  can  lead  to  the  appearance  of  suspended 
shock  waves  in  the  shock  layer  and  the  triple  points  on  the  shock 
wave.  The  programs  that  were  written  do  not  provide  for  the  ex¬ 
plicit  identification  of  characteristics  of  this  type,  and  there¬ 
fore  the  suspended  shock  waves  which  arose  will  be  "blurred"  by 
the  difference  scheme. 


Figure  16  compares  the  geometrical  pictures  of  flow  around 
a  body  with  a  depression  k  »  8,  ra  ■  Q.l  when  *  3,  4,  6  and  20. 
We  note  that  the  position  of  the  shock  wave  changes  strongly  with 
the  change  in  M»  from  3  to  6,  and  slightly  with  the  change  in 
from  6  to  2Q .  When  Mo,  increases,  the  sonic  line  is  displaced  to 
the  axis  of  symmetry  —  slightly  on  the  surface  of  the  body,  and 


quite  strongly  on  the  shock  waves.  The 
dimensions  of  the  local  zone  of  subsonic 
flow  in  the  depression  decrease  with  the 
increase  in  M^. 

Pigure  17  gives  a  comparison  of  the 
picture  of  flow  around  bodies  with  m  =  0 
(sphere),  m  =  0.1,  and  0.2  when  k  »  8, 
when  M,,,  =>  6 .  The  picture  of  flow  in  the 
shock  layer  with  deepening  of  the  depres¬ 
sion  changes  when  *  6  in  precisely  the 
same  way  as  when  =  3*  We  note  that 

unlike  the  case  of  M  ■  3,  the  local  zone 

00 

of  subsonic  flow  in  the  depression  occupies 
only  a  part  of  the  shock  layer,  without 
reaching  the  shock  wave. 

The  regressed  shock  waves  in  bodies 
which  have  an  identical  depression  depth 
(m  »  0.2),  but  different  depression  curva¬ 
ture  (for  values  of  the  parameter  K  *  8, 

4,  2)  are  plotted  in  Fig.  18  (for  *  6). 

The  shock  wave  approaches  the  body  with 
the  decrease  in  depression  curvature  and 
the  shock  waves  which  correspond  to  bodies 
with  different  values  of  the  parameter  K 
can  be  considered  approximately,  geometric¬ 
ally  similar.  We  also  note  that  the  bound¬ 
aries  of  the  local  subsonic  zone  in  the  de¬ 
pression  depend  weakly  on  curvature  K,  while 
the  transverse  dimension  of  this  zone  signifi¬ 
cantly  increases  with  the  increase  in  K. 


Fig.  16 
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Fig.  17 
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Pig.  18 


We  shall  dwell  on  the  distribution  characteristics  of  pressure 
through  bodies  with  a  contour  (4.1).  The  distributions  of  pressure 
p/pQ  along  the  surface,  depending  on  the  angular  coordinate  6,  are 
plotted  in  Pigs.  19  and  20.  The  characteristic  feature  of  these 
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Fig.  20 
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distributions  is  their  nonmonotonic  character.  The  increase  in 
pressure  in  the  region  corresponds  to  deceleration  of 

gas  on  part  of  the  surface  having  negative  curvature.  The  value 
of  the  minimum  in  pressure  distribution  decreases  with  the  in¬ 
crease  in  depth  or  curvature  of  the  depression,  while  the  value 
of  the  maximum  increases.  The  maximum  value  of  pressure  in  the 
depression  differs  from  pressure  in  the  critical  point  by  only 
about  1Q%  for  a  body  with  the  deepest  and  steepest  of  the  exam¬ 
ined  depressions  On  »  0.2,  K  »  8),  We  note  the  presence  of  a 
small  "gap”  in  the  distribution  of  pressure  immediately  before 
pressure  reaches  the  local  maximum  for  this  same  body.  A  similar 
effect  was  obtained  during  the  calculations  of  flow  around  bodies 
with  a  complex  shape  by  the  S,  K.  Godunov  method  in  a  work  [23]. 

The  effect  of  the  M  number  on  the  pressure  distribution 
along  the  surface  is  manifested  in  a  decrease  in  the  value  of 
the  local  pressure  minimum.  With  respect  to  the  local  maximum, 
then  for  the  case  m  ■  0.1,  its  value  noticeably  decreases,  but 
for  the  deeper  depression  m  *  Q.2,  it  remains  approximately  con¬ 
stant,  The  positions  of  both  the  minimum  and  the  maximum  of 
pressure  for  both  depressions  shift  toward  the  middle  with  an  in 
crease  in 
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The  calculated  distributions  of 
pressure  along  the  body  are  compared 
with  pressure  according  to  the  modi¬ 
fied  Newton  formula  in  Pigs.  19  and 
20. 


The  Newton  formula  describes  the 
actual  distribution  of  pressure  along 
the  body  at  hypersonic  M*,,  particular' 
ly  when  m  -  0.1,  well. 


*n  -  0,2  ,  <  '  A, 

Pig.  21 

m  *  0 . 2 ;  K  *  8 


-  flow-around  by  a  per- 

The  effect  of  equilibrium,  physico-  feet  gas  /j  -  t.«  i 

chemical  transformations  in  the  shock  _  flow-around  by  eauili- 

layer  on  the  picture  of  flow-past,  is  brium  0,1 

shown  by  Pig.  21.  It  appears  precisely 

the  same  as  in  the  case  of  flow  around  convex  bodies,  and  no  new 
qualitative  characteristics  arise  here.  We  note  that  the  size  of 
the  local  zone  of  subsonic  flow  in  the  depression  also  decreases 
with  the  decrease  in  thickness  of  the  shock  layer  caused  by  the 
effect  of  the  real  gas.  A  comparison  of  the  results  of  calculat¬ 
ing  flow  around  bodies  C4.1)  by  equilibrium  air  and  a  perfect  gas 
with  a  specially  selected  adiabatic  exponent  shows  that  the  method 
of  approximately  calculating  the  real  properties  of  air  during  the 
flow  around  the  facing  surfaces  of  bodies,  developed  In  [171,  is 
also  Justified  for  the  case  of  complexly  shaped  bodies  (4.1). 


We  shall  now  dwell  on  the  methodological  characteristics  of 
the  numerical  solution  to  the  problem  of  flow  around  convexo- 
concave  bodies.  A  significant  aspect  here  is  that  one  cannot  ob¬ 
tain  a  numerical  solution  for  flow  around  bodies  with  a  signifi¬ 
cant  depression  depth  (namely,  when  m  >  0.05)  without  introducing 
artificial  viscosity  into  the  difference  scheme  (1.3). 


When  making  calculations  of  the  flow  past  bodies  whose  curva¬ 
ture  does  not  change  sign  (Secs.  2  and  3),  it  was  unnecessary  to 


28 


resort  to  artificial  viscosity  and  the  coefficient  a  in  (l.3a) 
was  assumed  to  equal  zero.  On  the  other  hand,  during  calculations 
of  the  flow  around  bodies  with  depressions  when  a  *  0 ,  the  estab¬ 
lishment  process  diverges.  With  a  certain  minimum  value 
one  can  obtain  a  steady  state  solution,  but  the  distributions  of 
the  gasdynamic  parameters  in  the  region  of  the  depression  will 
have  a  sawtoothed  character.  With  a  certain  s.  >6~  ,  one  can 

smooth  these  distributions.  The  results  of  calculations  made  in 
this  section  ('Pigs.  13-21)  were  obtained  with  values  of  the  artifi¬ 
cial  viscosity  coefficient  •j-*..  The  value  of  a  increases  with 
the  increase  in  M  :  when  3 , s.*o.2a,  when  «.»o.s. 

OO^ 

Most  of  the  calculations  were  made  with  several  values  of 
the  parameter  a  in  order  to  identify  the  relationship  between  the 
numerical  solutions  and  this  value.  It  proved  that  the  value  of 
a  exerts  a  noticeable  effect  on  the  values  of  the  gasdynamic  values, 
which  determine  flow  in  the  depression.  Thus,  the  local  minimum 
in  the  pressure  distribution  along  the  body  (Figs.  19  and  20)  de¬ 
creases  with  the  increase  in  cr,  while  the  local  maximum  drops. 

The  accuracy  of  making  these  checks  in  the  region  which  corresponds 
to  a  depression  in  the  body  deteriorates  with  the  increase  in  a. 
However,  the  change  in  a  within  quite  broad  limits  (for  example, 
from  0.25  to  1  when  M  *  3)  has  no  noticeable  effect  on  the  numerical 
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solution  outside  the  depression  of  the  body  in  precisely  the  same 
way  as  it  does  on  accuracy  of  making  the  checks. 

With  a  value  «•<.  ,  the  accuracy  with  which  the  Bernoulli 

integral  is  satisfied  in  the  region  of  the  body  depression  and  the 
condition  of  entropy  constancy  comprises  approximately  2-M  %  when 
m  ■  Q.l,  and  approximately  5-8  %  when  m  *  0.2.  In  the  remainder 
of  the  shock  layer,  the  accuracy  of  making  checks  is  higher  and 
comprises  approximately  1-2  %.  These  data  pertain  to  calculations 
made  on  a  grid  with  11  x  21  points  when  9 k  *  90°;  the  accuracy  of 
calculations  made  on  a  9  x  17  grid  when  9 k  ■  70°  is  approximately 
the  same.  The  accuracy  of  calculations  made  on  the  9  x  17  grid 
when  ^  ■  9Q°  is  somewhat  lower. 
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In  conclusion,  one  should  note  that  the  effect  of  viscosity, 
as  a  rule,  can  be  taken  into  account  well  within  the  frameworks 
of  boundary  layer  theory  for  convex  bodies  whose  contour  curvature 
does  not  change  sign.  Therefore,  the  calculations  made  on  the 
basis  of  the  model  of  a  nonviscous  gas  give  the  correct  pressure 
distributions  on  the  surface  of  the  streamlined  body,  as  well  as 
the  correct  values  of  the  other  gasdynamic  parameters  in  the  entire 
field  behind  the  closed  boundary  layer  zone.  In  the  examined  case 
of  convexo-concave  bodies,  high  positive  pressure  gradients  on  the 
surface,  as  in  Pigs.  19  and  20,  can  lead  to  boundary  layer  separa¬ 
tion  and  gas  flow  in  the  depression  will  be  significantly  viscous. 
This  is  specifically  indicated  by  comparing  the  calculation  and 
experimental  data  on  the  flow  around  convexo-concave  bodies  in 
[23].  Therefore,  the  results  of  solving  the  problem  of  flow  around 
bodies  of  such  shape,  evidently,  should  be  viewed  as  the  first  ap¬ 
proximation  to  the  real  case. 

Section  5.  Nonequilibrium  gas  flows  around  surfaces 

In  works  [8-11],  the  establishment  method  was  applied  to  cal¬ 
culating  the  flow  around  spherical  surfaces  of  varying  radii  and, 
specifically,  around  a  plane  surface  by  supersonic,  nonequilibrium 
gas  flows.  Calculations  were  made  according  to  the  IM  program 
Cll  x  21  points). 

In  [8-10],  the  investigation  was  made  for  free-stream  flows 
which  are  nonuniform  with  respect  to  value  and  direction;  cases  of 
flow  engendered  by  the  supersonic  spatial  source,  and  a  stream 
flowing  In  a  vacuum  from  an  axisymmetrical  nozzle,  are  examined. 

The  flow  of  a  compressed  gas  from  a  spatial  source  exists  out¬ 
side  a  sphere  with  a  radius  R*,  in  which  the  velocity  of  the  flow 
is  equal  to  the  velocity  of  sound.  The  Mach  number  M  at  a  distance 
r  from  the  center  of  the  source  can  be  determined  from  the  equation 


(r 


(5.1) 


30 


With  hypersonic  M  ,  i.e.,  at  great  distances  from  the  source, 
one  can  use  the  approximate  formula 


The  velocity,  density  and  pressure  are  determined  according  to 
the  known  M  with  the  aid  of  the  formula  of  the  isentropic  flow. 

The  parameters  of  the  gas  in  the  supersonic,  freely  expand¬ 
ing  stream  were  taken  from  calculations  made  by  the  grid  method 
[15]  and  the  method  of  characteristics  [20].  We  note  that  the  flow 
in  such  a  stream  at  great  distances  from  the  nozzle  has  the  follow¬ 
ing  properties  [20]:  the  flow  lines  are  similar  to  straight  lines, 
while  the  distributions  of  parameters  along  them  satisfy  the  source 
law  C5.1)  or  C5.2)  with  good  accuracy;  the  center  of  the  source  is 
located  in  the  nozzle  cross  section.  However,  the  radius  of  the 
equivalent  source  changes  with  the  transition  from  one  stream 

line  to  the  other.  In  this  case,  the  radius  *•»  changes  weakly 
in  the  region  which  encompasses  the  axis  of  symmetry  and  is  bounded 
by  a  certain  characteristic  surface  y  =  y*Cx)  C.Fig.  22).  Intensive 
gas  expansion  occurs  in  the  indicated  region  in  the  centered  rarefac¬ 
tion  wave,  proceeding  from  the  edge  of  the  nozzle;  a  sharp  decrease 
in  the  value  of  *•»  corresponds  to  it. 

Figure  22  plots  the  geometrical  picture  of  flow  around  a  flat 

wall  by  a  stream  flowing  from  a  nozzle  with  a  uniform  flow  at  the 

outlet  and  a  Mach  number  M  =2.5;  the  wall  is  located  a  distance 

ot 

of  50  r^  from  the  cross  section  of  the  nozzle  CrQ  is  the  nozzle 
radius).  In  Fig.  22,  the  dash  line  plots  the  picture  of  flow  around 
a  plane  by  a  flow  from  a  source  whose  radius  R#  is  equal  to  the 
radius  of  the  equivalent  source  *„  for  flow  on  the  stream  axis 

.  In  the  region  between  the  axis  and  the  characteristic 
y*(x),  the  shock  waves  for  both  flows  around  the  object  differ  little 
from  each  other.  Deviations  in  the  field  of  the  gasdynamic  values 
in  the  shock  layer  are  somewhat  greater,  which  can  be  seen  by  compar¬ 
ing  the  positions  of  the  sonic  lines;  to  the  point,  dimensions  of 
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the  subsonic  regions  in  both  places  are  ap¬ 
proximately  equal. 

Differences  between  the  two  flows 
around  objects  become  significant  about  the 
characteristic  y*Cx).  In  the  case  of  the 
source,  the  velocity  of  the  gas  behind  the 
shock  wave  monotonically  increases  with 
distance  from  the  axis.  In  the  case  of  the 
stream,  a  sharp  decrease  in  pressure  and 
density  in  the  rarefaction  wave  above  y#(x) 
leads  to  the  intensive  deviation  of  the 
shock  wave  toward  the  nozzle  and  to  an  in¬ 
crease  in  the  angle  between  the  shock  wave 
and  the  direction  of  the  free-stream  flow. 
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Because  of  this,  velocity  behind  the 
shock  wave  once  again  begins  to  decrease  and 
reaches  the  subsonic  values;  in  this  case, 
a  second  subsonic  zone  adjacent  to  the  shock 
wave  (region  ABC  in  Fig.  22)  forms  in  the  shock  layer 


Fig.  22 

streamline  flow 


-  stream  flow  from 

source 


The  presence  of  this  region  has  no  effect  on  the  solution  in 
the  subsonic  and  transonic  regions  of  the  shock  layer  which  ‘surround 
the  axis  of  symmetry.  However,  when  the  Ma  number  increases  Cor 
with  the  increase  in  y) ,  the  second  sonic  point  B  rapidly  descends 
along  the  shock  wave  and  with  a  certain  value  of  Ma,  one  has  a 
closure  of  the  two  subsonic  regions.  A  narrow  zone  of  subsonic 
flow  forms  In  the  shock  layer,  extending  along  the  shock  wave. 

It  was  stated  in  Sec.  1  that  the  correct  statement  of  the 
problem  of  flow  around  a  body  with  a  regressed  shock  wave  on  a 
line  which  limits  the  calculated  region  of  flow  from  above  requires 
that  the  normal  component  of  gas  velocity  to  this  line  should  exceed 
the  velocity  of  sound.  Unsuccessful  efforts  were  made  to  obtain  a 
numerical  solution  with  such  a  property  to  the  problem  of  flow  around 
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a  flat  wall  by  a  stream  in  the  examined 
case.  The  velocity  was  subsonic  in  the 
calculations  that  were  made  in  a  certain 
part  of  the  boundary  line.  Naturally, 
the  question  of  the  validity  of  such  cal¬ 
culations  arises. 

In  order  to  prove  that  the  solutions 
obtained  in  this  case  correspond  to  reality, 
the  following  numerical  experiment  was  per¬ 
formed.  Plow  around  a  plane  by  a  flow  from 
a  source  whose  radius  changed  depending  on 
angle  0  according  to  an  assigned  law  was  cal¬ 
culated.  In  this  case,  the  flow  field  was 

similar  to  flow  in  a  stream  when  M  =2.65. 

a 

The  calculation  was  made  in^the  spherical  co¬ 
ordinate  system  with  its  center  in  the  center 
of  the  source  with  positions  of  the  outer  ray 


Fig.  23 

**»•  O.Mj - #,  -  O.T4  | 

- O.tM, 

<C.»  U.74J-  -  -  0.*»» 

•>...»«  •  0.7». 


0,  =  0.64  and  0.74.  The  fields  of  flow  in  the 

A 

boundary  layer  coincided  with  high  accuracy  for  both  cases  (Fig.  23). 


Attempts  to  obtain  a  numerical  solution  when  0.  >  0.74  were  unsuccess¬ 
ful. 


Calculation  was  further  made  for  a  case  when  the  radius  of  the 
source  remained  constant,  beginning  with  a  certain  value  0  =  .  . 

In  the  calculations,  it  was  assumed  that  ®const  =  0.74  and  0.79  when 
9^  =  0.84.  in  the  free-stream  flow  when  9  >  9const,  the  intensive 
expansion  of  gas  ceases  and  the  factors  which  cause  the  shock  wave 
to  deviate  sharply  toward  the  nozzle  disappear.  Velocity  behind  the 
shock  wave  rapidly  reaches  supersonic  values  and  the  outer  ray  0^ 
wholly  lies  within  the  supersonic  region. 


By  comparing  the  geometrical  pictures  of  flow  around  a  plane 
in  Fig,  23  for  all  three  cases  <  •.  .  0.74  4  *.  -  qj*.  ax  i 
one  can  see  that  up  to  9  =  9cons.t*  shock  waves  and  the  sonic  line 
do  not  practically  differ  from  each  other,  and  they  diverge  signifi¬ 
cantly  only  when  0  >  6const* 


Consequently,  the  perturbations  introduced  above  into  the 
flow  9*  ®const  are  Poetically  rot  transmitted  to  the  subsonic 
region.  This  fact  can  be  explained  because,  although  the  veloci¬ 
ties  are  subsonic  in  a  certain  region  of  the  boundary  line,  they 
are  still  comparable  in  value  to  the  speed  of  sound.  Because  of 
this,  information  through  the  subsonic  corridor  adjacent  to  the 
shock  wave  is  not  practically  transmitted  upstream.  Thus,  one  can 
obtain  a  solution  which  corresponds  to  reality  without  introducing 
the  outer  ray  into  the  region  of  supersonic  velocities. 

However,  such  resistance  to  perturbations  does  not  exist 
throughout  the  entire  subsonic  "corridor."  Velocity  in  the  shock 
wave,  having  reached  a  maximum,  decreases;  with  distance  from  the 
region  of  highest  velocities,  the  subsonic  flow  in  the  "corridor” 
becomes  more  and  more  sensitive  to  perturbations  arriving  from  the 
stream  boundary. 

Thus,  a  certain  region  of  high,  subsonic  velocities  exists 
in  the  "corridor,"  weakly  transmitting  perturbations  upstream.  It 
is  a  unique  analog  of  the  supersonic  region  between  two  sonic  points 
on  the  shock  wave  in  the  case  of  flow  with  two  subsonic  regions . 

We  note  here  that  a  suspended  shock  wave  forms  in  the  stream  flowing 
from  the  axisymmetrical  nozzle  into  space  with  counterpressure;  flow 
in  the  stream  is  precisely  the  same  as  with  flow  in  a  vacuum  in  the 
region  between  this  wave  and  the  axis.  If  a  flat  obstacle  is  placed 
in  such  a  stream,  then  the  suspended  shock  wave  can  intersect  the 
regressed  shock  wave  before  the  obstacle  at  the  point  where  velocity 
behind  the  regressed  wave  is  subsonic.  Nevertheless,  on  the  basis 
of  the  conducted  Investigations,  one  can  assert  that  the  perturba¬ 
tions  introduced  by  the  suspended  shock  wave  have  practically  no 
effect  on  flow  in  a  significant  part  of  the  shock  layer  ahead  of 
the  obstacle. 

Figures  24  and  25  compare  the  results  of  calculating  the  plane 
of  the  stream  with  the  data  of  experiment  [24].  These  data  pertain 
to  the  case  of  a  stream  of  air  C y  *  1.4),  flowing  from  a  nozzle  with 


Fig.  24 

1.  Calculation  when  <J>  =  0 ;  1. 

2.  calculation  when  4>  =  15°;  2. 

-  experiment  [24]  o 

a  Mach  number  on  the  section  =  5  for  a  pressure  drop  in  the  nozzle 
receiver  and  the  surrounding  space  o;pf  .  )t>*  (thence  io*  ).  The 

position  of  the  regressed  shock  wave  ahead  of  the  flat  wall  located 
a  distance  D  *  32.6  ra  from  the  nozzle  cross  section  obtained  in  the 
experiment,  is  plotted  by  the  dash  line  in  Fig.  24.  Here,  too,  the 
x's  mark  points  of  intersection  of  the  regressed  wave  before  the 
wall  with  a  suspended  wave  in  the  stream  and  a  stream  boundary.  In 
[24],  there  are  no  indications  relative  to  the  character  of  flow  at 
the  nozzle  outlet  (it  is  only  stated  that  the  stream  flows  from  a 
long,  isentropic  nozzle}.  Therefore,  the  experimental  data  are  com¬ 
pared  with  the  calculation  data  for  the  case  of  a  stream  flowing 
from  a  nozzle  with  a  uniform  flow  at  the  outlet  (.4)  =  0)  and  from  a 
conical  nozzle  with  a  coning  angle  4>  =  15°. 

The  experimental  and  calculated  positions  of  the  regressed 
wave  coincide  well  at  distances  from  the  axis  of  symmetry  of  r*  o,4Pr 
When  in  the  experiment,  one  observes  an  extremely  sharp  devia¬ 

tion  of  the  regressed  shock  wave  toward  the  nozzle;  the  numerical 
solution  when  y  >0.4  D  for  a  stream  flowing  from  a  nozzle  with 


Fig.  25 

Calculation  when  <j>  *  0  ; 
calculation  when  4>  =  15°; 
experiment  [24] 


*  0  could  not  be  obtained. 


The  coincidence  of  the  calculation  and  experimental  data  on 
determining  pressure  on  the  flat  wall  (Fig.  25)  is  satisfactory. 
Checks  of  accuracy  of  the  obtained  solutions  showed  that  the  error 
in  satisfying  the  Bernoulli  equation  •  ojs  *  2%  ;  the  scatter  of 
values  of  the  entropy  function  on  different  flow  lines  comprised 
0.5  -  3*.  The  results  pertained  to  the  case  of  a  stream;  in  the 
case  of  the  source,  errors  «*  and  *3  did  not  exceed  0.5%. 

A  work  [11]  examined  flow  around  a  sphere  parallel  but  not 
equal  in  value  to  a  supersonic  flow,  which  can  be  viewed  as  two  co¬ 
axial  streams  which  have  a  transition  zone  between  them.  Calcula¬ 
tions  were  made  both  for  a  case  when  velocity  Voo^  of  the  internal 
stream  is  greater  than  the  velocity  of  the  outside  stream,  and 
for  the  case  when  •  The  transitional  zone  occupied  a  regio 

*,  .*«*,.»  .  The  velocity  profile  V*,  of  the  transition  zone  was  as¬ 

signed  by  the  formula 


f  /y-*, ,w/'  - 

*  '  j  Vm,9Vmt 

(~  )  1  when  •«  *  V-’ 


when 


(5.3) 


(5.4) 


Static  pressure  and  total  enthalpy  throughout  the  entire  flow  were 
assumed  to  be  constant. 


The  radius  of  the  outside  stream  is  considered  to  be  suffi¬ 
ciently  great  such  that  the  outer  boundary  of  the  stream  has  no  ef¬ 
fect  on  the  calculated  region  of  flow.  Most  calculations  are  made 
when  *  0.4,  h  »  0.2. 

The  geometrical  pictures  of  flow  around  a  sphere  by  coaxial 
streams  are  plotted  in  Figs.  27  and  28.  The  Mach  number  of  the 
internal  stream  in  the  cases  Mw^  *  3  is  presented  in  Figs.  26  and 
27.  In  Fig.  26,  velocity  of  the  outside  stream  is  less  than  that 
of  the  inside  stream  (M^  *  2.5  and  2),  while  Fig.  27  shows  the 


reverse  case  CMoo2  *  3.5  and  4).  The  examined  flows  are  character¬ 
ized  by  a  more  complex  shape  of  the  shock  wave  and  the  sonic  line 
in  comparison  with  the  flow  of  a  uniform  flow  around  a  sphere.  We 
note  that  when  M092  <  ,  one  has  a  significant  increase  in  part  of 

the  region  of  effect  adjacent  to  the  body,  but  when  Mao2  >  M^,  one 
has  an  increase  in  the  part  of  the  region  of  effect  adjacent  to  the 
shock  wave.  A  region  of  the  shock  wave  which  contains  the  inflection 
points  corresponds  to  the  transition  zone  in  the  free-stream  flow 
C5.3)  • 


The  gas  parameters  behind  the  shock  wave  change  very  sharply 
here.  Calculations  show  that  when  M^  is  small  and  Nk2  >  one 

observes  the  appearance  of  a  local  supersonic  zone  in  this  region 
adjacent  to  the  shock  wave  (.for  example,  when  =  1.3,  N^2  *  1.5) 
When  decreases,  this  zone  blends  with  the  basic  supersonic 
region  (for  example,  when  *  1.1,  r^2  =  1.5). 

We  note  that  regression  of  the  shock  wave  from  the  sphere  on 
the  axis  of  symmetry  during  flow  passed  by  coaxial  streams  with 
Mach  numbers  and  increases  in  comparison  with  regression 
ahead  of  the  sphere  streamlined  by  a  uniform  flow  with  »  Mm1 
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both  with  a  decrease  and  an  increase  in  by  comparison  with 
(Figs.  26  and  27).  The  increase  in  in  comparison  with  leads 

to  a  decrease  in  pressure  on  the  sphere  and,  on  the  contrary,  a  de¬ 
crease  in  1^2  leads  to  an  increase  in  pressure  on  the  body. 

In  the  calculations  that  were  performed,  the  accuracy  of 
satisfying  the  Bernoulli  equation  was,  as  a  rule,  *6-1-2%  •  The 

scatter  of  values  of  the  entropy  function  on  different  flow  lines 
was  within  limits  of  0.5  -  b%.  The  accuracy  of  the  calculations 
increases  with  the  decrease  in  the  difference  of  velocities  of  the 
internal  and  external  streams  with  a  fixed  width  of  the  transition 
zone.  Accuracy  also  increases  with  the  increase  in  the  width  of 
this  zone  with  fixed  stream  velocities. 

Section  6.  Solving  the  direct  problem  for  a  Laval  nozzle 

The  problem  of  flow  in  a  Laval  nozzle  of  assigned  shape  was 
solved  in  a  work  [71  with  the  aid  of  the  establishment  method. 

Let  y  *  G (x )’  be  the  contour  equation  of  the  nozzle,  and  F  be 
the  length  of  the  nozzle.  We  insert  the  normed  variables 

t.s/f  .  :•  */*<*'  ■  (6.1) 

The  outlet  cross  section  C.K  -  1)  is  arranged  such  that  condi¬ 
tion  u  >c  is  satisfied  everywhere  on  it;  then  when  5  »  1,  one  need 
not  impose  boundary  conditions;  relative  to  nozzle  intake  (£  *  0), 
more  shall  be  stated  below.  Let  the  initial  condition  of  a  gas  be 
assigned  at  a  moment  t  ■  0.  t*  signifies  a  finite  moment  of  time 
and,  consequently,  the  region  of  interest  to  us  is  such: 

0«  t  «  tm  ,  O  »  4  *  1,  O  «  t  «  I. 

We  shall  briefly  describe  the  algorithm  of  calculations  with 
direct  and  inverse  acceleration  when  solving  a  system  of  differen¬ 
tial  equations  similar  to  (l.M).  For  the  sake  of  simplicity,  we 
shall  assume  that  only  one  point  exists  such  that  for  all 
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4-  « «.(<"-  o,i  1)  inequality  u  <  c  is  satisfied.  This  ensures 

stability  of  trial  run  (Sec.  1),  but  when  <m  »«.  ,  u  >  c.  With 
the  aid  of  an  auxiliary  difference  scheme  described  in  [5],  we 

shall  determine  the  values  of  the  gasdynamic  parameter  u  and  *.e»„ 
when  C  s  on  a  time  layer  t  =  tn  .  We  shall  write  the  trial 
and  error  relationship  in  the  following  form  when  £  *  £ ^ 


(6.2) 


(the  index  ”1”  pertains  to  line  r  =  r^).  We  shall  assume  *.•<>  } 

>*.  •  -oa  *o.8  ;  when  /  ■■<««  +om*  . 


Further,  with  the  aid  of  recurrent  relationship  similar  to 
(1.7),  we  shall  determine  the  trial  and  error  coefficients  when 
m  =  k  -1,  k  -2,...  0,  i.e.,  we  shall  perfect  the  direct  trial  and 
error  run  in  the  direction  toward  the  nozzle  intake.  Having  de¬ 
termined  the  gasdynamic  parameters  when  m  *  0  (see  below),  we  can 
successively  find  their  values  when  m  =  1,2,...,  k,  i.e.,  complete 
the  inverse  trial  run.  When  m  *  k  +  1,  k  +2,...  M,  the  inverse 
trial  run  is  accomplished  using  only  the  system  of  difference 
equations  similar  to  (1.4),  without  using  the  trial  and  error  rela¬ 
tionship  . 


Calculations  are  made  for  all  ♦,  (  (  -o,i . o  .  During  cal¬ 

culations  on  the  nozzle  contour  Cl  *  L),  naturally,  the  boundary 
condition  of  non-flow  is  used.  By  transferring  from  one  time  layer 
to  another,  we  find  the  flow  field  in  the  entire  region. 

Most  of  the  calculations  were  carried  out  with  a  difference 
grid  having  7  x  64  points  (7  in  the  ♦  direction,  and  64  in  the  5 
direction) . 


Figure  28  gives  certain  results  of  calculations  which  pertain 
to  the  steady  state  flow  regime.  The  nozzle  1,  shown  in  Fig.  28, 
was  obtained  as  the  result  of  solving  the  inverse  steady  state 
problem  in  the  work  of  U.  G.  Pirumov  [191.  All  linear  dimensions 
are  related  to  the  radius  of  the  nozzle  throat.  The  figure  shows 
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I  -  nozzle  from  a  work  [191; 

II  -  nozzle  with  a  contour  (6.3); 

_  distributions  of  parameters 

on  the  side  of  the  nozzle 

and  the  shape  of  the  sonic  line; 

-  distributions  of  parameters 

on  the  axis  of  the  nozzle; 

0  -  according  to  the  results  of  a 
work  [19  1 . 


Fig.  28 


the  nozzle  region  close  to  the  outlet,  and  for  convenience,  the 
origin  of  calculation  is  shifted  along  the  x  axis  in  comparison 
with  the  equation  presented  in  the  calculation  (for  the  given 
nozzle,  the  distance  from  the  intake  cross  section  to  the  throat 
X  *  9.8971. 


Figure  28  gives  the  position  of  sonic  lines  and  the  distribu¬ 
tions  of  parameters  along  the  axis  and  the  contour.  Pressure  is  re¬ 
lated  to  the  value  pQ,  velocity  to  the  value  fpTTt?  ,  where  P0»Pq 
total  pressure  and  density;  y  *  1.4  »  const. 

Of  the  flow  characteristics,  we  note  a  slight,  approximately 
0.1  -  0.4  of  pQ  increase  in  pressure  on  the  side  (in  comparison  with 
its  value  when  x  »  0),  which  exists  before  the  beginning  of  the  in¬ 
tensive  acceleration  of  gas.  It  is  almost  unnoticeable  in  Fig.  28 
due  to  the  smallness  of  scale. 

The  transition  surface  through  the  velocity  of  sound  signifi¬ 
cantly  differs  from  the  rectilinear  one  in  the  regime  of  steady  state 
flow.  The  sonic  line  extremely  rapidly  (in  comparison  with  flow  on 
the  whole)  acquires  a  shape  similar  to  the  stationary  one. 

The  difference  between  the  results  obtained  during  solution 
of  the  direct,  unsteady  problem  and  the  data  [19)  in  the  vicinity 
of  the  throat  comprises  a  value  on  the  order  of  1% .  This  difference 
is  basically  due  to  certain  inaccuracies  in  assigning  the  nozzle 
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contour.  We  still  note  that  already  when  *»o.e  ,  flow  noticeably 
differs  from  one-dimensional  flow  Cparticularly  in  the  vicinity  of 
the  critical  cross  section) . 

When  solving  the  direct,  unsteady  problem  for  the  nozzle  de¬ 
picted  in  Pig.  28,  the  initial  condition  of  the  gas  was  determined 
from  a  calculation  made  by  a  method  [191.  In  the  process  of  solu¬ 
tion,  steady  state  values  of  all  gasdynamic  parameters  were  main¬ 
tained  at  the  intake  ( Z  *  0).  Inasmuch  as  the  basic  goal  of  the 
calculation  of  flow  in  the  given  contour  was  a  comparison  with  the 
data  of  [191,  the  indicated  method  seems  the  natural  one. 

The  results  of  a  calculation  of  flow  in  a  nozzle  whose  con¬ 
tour  was  assigned  by  the  equation 


f  2.32 

[1.06  *  0.08  »*^{*-7i.a 


when  o  «  x  •  njoa 

51  when  7>.fi33‘  **  74.825. 


(6.3) 


were  also  cited  in  [7] .  Figure  28  shows  this  contour.  The  initial 
state  of  the  gas  in  this  case  was  determined  from  the  one-dimensional 
theory  with  the  additional  assumption  of  a  linear  law  of  change  of 
the  angle  of  inclination  of  the  velocity  vector  between  the  contour 
and  the  axis  in  cross  section  x  *  const.  Steady  state  parameter 
values  were  maintained  in  the  intake  to  the  nozzle  (£  =  0).  Such 
conditions  when  £  *  0  are  essentially  equivalent  to  conditions  at 
infinity,  inasmuch  as  the  length  of  the  subsonic  part  was  selected 
such  that  perturbations  over  time  t  =  tB  did  not  reach  from  the  in¬ 
take  to  the  region  of  the  contour  shown  in  Pig.  28,  which  is  of 
interest  to  us.  In  a  known  sense,  the  obtained  solution  simulates 
(with  an  accuracy  up  to  the  initial  data  when  t  =  Q)  the  operation 
of  a  unique  wind  tunnel. 

We  note  still  another  statement  of  the  examined  problem,  based 
on  the  model  of  a  tank-type  installation.  A  quite  large  reservoir 
filled  with  gas  is  attached  in  one  way  or  another  to  the  investigated 
contour.  If  G  >>  g,  where  G  is  the  mass  of  the  gas  contained  in  the 


container*  and  g  is  the  mass  of  the  gas  which  passes  through  the 
nozzle  over  a  time  t,*,  necessary  for  entry  into  the  steady  state 
regime,  then  as  the  result  of  the  calculation  when  t  =  t^,  we  ob¬ 
tain  the  steady  state  solution.  In  this  case,  the  problem  is 
completely  closed  and  no  additional  conditions  except  the  boundary 
conditions  of  non-flows  are  required. 


We  note  in  conclusion  that  the  method  of  establishment  was 


used  in  a  work  [12]  for  calculating  unsteady  flow  in  a  plane,  super¬ 
sonic  stream  interacting  with  external,  acoustical  waves  of  assigned 
frequency  and  amplitude.  This  problem  was  examined  in  connection 
with,  the  investigation  of  the  phenomenon  of  sonic  radiation  and  a 
discrete  frequency  by  supersonic  streams  flowing  into  a  flooded 
space.  The  investigation  of  the  relationship  between  solutions 
which,  describe  the  perturbed  flow  in  the  stream  and  the  amplitudes 
of  the  external  acoustical  waves  given  in  [12]  made  it  possible  to 
draw  certain  conclusions  about  the  mechanism  of  limiting  pulsations 
of  the  stream  which  emits  sound  at  a  discrete  frequency,  leading 
to  the  establishment  of  self-excited  vibrations. 
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A  DIFFERENCE  SCHEME  OF  THE  THIRD  ORDER  OF  ACCURACY 

FOR  CALCULATING  TWO-DIMENSIONAL  FLOWS  WITH  CONTACT 

SEPARATION 

V.  V.  Yeremin,  Yu.  M.  Lipnitskiy 

The  presently  employed  difference  methods  of  through  cal¬ 
culation  basically  have  the  first  or  second  order  of  approximation. 
Difficulties  of  increasing  the  order  of  the  difference  schemes 
consist  in  the  necessity  of  using  a  significantly  larger  number 
of  points  of  the  initial  field  than  in  schemes  of  the  first  or 
second  order,  which  has  a  significant  effect  on  calculation  time. 
This  has  an  effect  particularly  when  the  order  of  approximation  is 
increased  for  two-dimensional  schemes  .  The  increased  requirements 
of  the  accuracy  of  determining  the  gasdynamics  of  flows,  on  the  one 
hand,  and  the  increase  in  the  speed  and  memory  of  the  computer,  on 
the  other,  advance  the  problem  of  developing  and  adopting  methods 
of  increased  accuracy  for  solving  complex  gasdynamic  problems. 

Works  [1,2,3]  have  been  devoted  to  the  development  of  such  schemes. 
These  works  have  suggested  difference  schemes  of  increased  accu¬ 
racy  for  solving  one-dimensional  gasdynamic  problems.  The  fact 
that  the  schemes  of  a  higher  order  of  accuracy  produce  significant¬ 
ly  less  oscillation  near  the  separation  and  significantly  reduce 
the  zone  of  its  "blurring”  is  a  significant  result  of  these  works. 
The  numerical  results  for  two-dimensional  problems  of  gasdynamics, 
obtained  with  the  aid  of  third  order  devices,  have  not  yet  been 
published,  although  such  difference  schemes  based  on  a  25  point 

pattern  have  been  constructed  in  works  [2,  3J • 

This  article  suggests  a  quite  compact  device  of  the  third 
order  of  approximation  for  solving  two-dimensional  problems  of  gas- 
dynamics  based  on  a  1 6  point  pattern  and  calculation  examples  are 
given  in  it . 

We  shall  write  a  system  of  gasdynamic  equations  for  which 
the  difference  scheme  is  constructed  in  the  divergent  form 


(1) 


where 
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The  symbols  for  the  gasdynamic  functions  are  the  standard  ones. 
System  (1)  is  closed  by  the  equation  for  specific  energy,  which  has 
the  following  form  for  a  perfect  gas: 

t  i 

'-’‘-i'T- ~  ■ 


We  insert  a  difference  grid  with  steps  Ax,  Ay,  At  in  a  space 
Cx,  y,  t).  In  addition  to  the  main  •  difference  grid,  we  insert 
an  auxiliary  one  (Fig.  1),  where  -*  •  '"•-*/  •  »•*•■»< 

We  shall  use  both  integral  and  half-integral  values  of  the  spatial 
indices  k  and  ,  having  designated  values  of  functions  at  nodes 
of  the  basic  grid  ,  and  at  nodes  of  the  auxiliary  grid 

The  solution  to  the  system  Cl)  on  the  step  shall  be  sought  as  the 
result  of  the  following  iteration  process: 


/  /  *  I  /«*»  .*•*  ftl  ,  tv  f  ...»  I 


In  order  to  calculate  each  /$*  ,  we  shall  use  values 


/  *./ 


We  Insert  the  symbols  for  the  averaging  operators  and  for  taking 
the  difference  according  to  x  and  y  : 


,/•¥»  ,(•!/!  f ’•  ft,  /-  <ft)]  • 


(2) 


We  shall  require  that  the  sought  difference  scheme  be  a  natural 


continuation  of  schemes  of  the  first  and  second  orders,  i.e.,  tha 
the  first  step 
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and  the  second  step 

.  cu  ,  1,1  t,  6  ™  Cl 

i'j'J™  shall  be  sought  in  the  form: 
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The  unknown  coefficients  of  formula  (5)  need  to  be  selected 
such  that  the  expansion 


/../  i  it*  i 


where  the  time  derivatives  are  expressed  with  the  aid  of  equation 
Cl)  through  functions,  and  their  derivatives  according  to  the 


tW 

spatial  coordinates  coincide  with  the  expansion  '».t  in 


powers  of  At  up  to  terms  of  the  order  Atf,  inclusively. 


Making  the  expansions  is  an 
extremely  laborious  task  because 
of  the  nonlinear  character  of  the 
relationship  F (f),  G(f)  and  H(f). 
Omitting  these  cumbersome  inter¬ 
mediate  calculations,  we  give  the 
final  result: 


i.e. , 
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There  can  be  infinitely  many  schemes  of  the  third  order  based  on  the 

16  point  pattern,  inasmuch  as  the  equation  system  for  de¬ 
termining  the  scheme  coefficients  has  ten  equations  with  16  un¬ 
knowns.  Of  all  the  solutions  to  this  system,  of  the  most  practical 
value  are  schemes  which  satisfy  (3)  and  (4)  and  have 

the  least  number  of  components  in  the  right-hand  side.  The  form 
(5).  was  chosen  in  accordance  with  these  conditions. 

We  shall  limit  ourselves  to  an  examination  of  a  linear, 
differential  equation  of  the  following  type  during  the  investiga¬ 
tion  of  stability  of  the  suggested  difference  scheme  : 


The  difference  scheme 
ing  way : 


m  *S  *Jfm°  •  /•/««. y.t) 

for  this  equation  is  written  in  the  follow- 
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Having  expressed  values  /.?  by  the  values  Js,  ,  and  by 

examining  a  solution  of  the  type  fit  •  «•'*••<*>  for  determining 

the  region  of  stability  of  the  difference  scheme  >  we  obtain  an 
inequality:  V*  i  ,  where 


f  n-v  '  ■ 
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The  analytical  investigation  of  this  inequality  is  extremely 
difficult,  and  therefore  the  region  of  its  solution  was  determined 
numerically .  With  fixed  oo^  and  max  IM  was  sought  according 

to  a  and  3  . 

The  relationship  between  max  |X|  and  cj^  and  was  constructed 
as  the  result  of  this.  The  investigation  of  this  relationship  showed 
that  the  scheme  is  stable  when  «»*♦  «/<**,  where  R  *  1/3. 

As  an  example,  Pigs.  2-10  give  the  results  of  calculations 
of  flow  around  a  flat  tip  by  uniform  and  nonuniform  supersonic  flows 
of  a  perfect  gas.  The  general  picture  of  flow  in  the  case  of  a 
homogeneous  free-stream  flow  with  a  number  m_-  3, .  T .  M  is  shown  in 
Pig.  2.  It  gives  the  horizontal  lines  which  are  designated  by  the 
indices  1,  2,  3,  4,  and  the  distributions  of  the  gasdynamic  func¬ 
tions  alone  which  lines  are  drawn  in  Pigs.  3-6.  The  presented  values 
of  pressure  p,  density  p,  and  the  projections  of  the  velocity  vectors 
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Fig.  2 


r 

!-i 

u,  v  are  respectively  related  to  the  doubled  velocity  head  ?-v-  ,  , 

the  density  head  p^,  and  the  velocity  head  of  the  free-stream 

flow.  • 

M 

i 

Figures  3  and  4  clearly  show  regions  of  unperturbed  flow  \ 

and  the  flow  between  the  shock  wave  and  the  body,  and  the  jump  is 
"blurred"  over  2-3  points  of  the  difference  scheme.  Errors  in  p 

satisfying  the  Hugoniot  conditions  did  not  exceed  1%.  j 

I 

A  characteristic  region  of  rarefaction  appears  in  the  dis¬ 
tribution  of  the  gasdynamic  parameters  in  Figs.  5  and  6,  inasmuch 
as  the  third  and  fourth  rays  (Fig.  2)  intersect  the  fan  of  the 
rarefaction  waves  proceeding  from  the  body  contour  break  point. 

From  the  viewpoint  of  investigating  the  possibilities  of  the  ■ 

scheme,  the  examples  of  calculating  flow  around  a  body  by  a  non-  [ 

homogeneous,  supersonic  flow  with  contact  separation  is  most  inter- 

!■ 

esting.  The  corresponding  picture  of  flow  is  shown  in  Fig.  7,  where  i: 

the  dash  line  depicts  contact  separation,  while  pressure  and  density 
in  the  free-stream  are  constant  and  the  velocity  value  changes  upon 
transferring  through  it  two  times  (v  r  •  m>  .  Contact 

separation  in  the  examined  case  falls  on  the  shock  wave  below  the 
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sonic  line,  which  corresponds  to  flow  with  a  number  M  *  3  CFigs . 

2,  7).  This  leads  to  a  significant  restructuring  of  the  flow  pic¬ 
ture  and  to  the  formation  of  a  trinary  point  on  the  shock  wave  (the 
sonic  line  and  contact  separation  on  the  shock  wave  have  a  common 
point).  In  connection  with  the  decrease  in  the  vertical  size  of 
the  region  of  effect,  shock  wave  regression  increases  and  the  sonic 
line  deforms  significantly.  The  distributions  of  pressure  and  the 
longitudinal  velocity  component  along  the  cross  sections  1,  2,  3 
and  4  (Fig.  7)  are  shown  in  Figs.  8-9.  Pressure  remains  constant 
in  cross  section  1,  being  slightly  perturbed  in  the  contact  separa¬ 
tion.  The  second  cross  section  intersects  the  shock  wave  at  a  small 
angle,  which  leads  to  ''blurring”  of  the  shock  wave  over  a  large  num- 
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ber  of  grid  points. 

The  flow  is  still  not  perturbed 
in  the  vicinity  of  the  contact  discon¬ 
tinuity,  and  its  pressure  change  is 
slight.  The  contact  discontinuity  is 
situated  in  the  region  behind  the  shock 
wave  on  rays  3  and  4  and  the  character¬ 
istic  pressure  buildup  behind  it,  associ¬ 
ated  with  the  increase  in  pressure  after 

the  shock  wave,  proportional  to  the  square  of  the  M  number  is  notice¬ 
able.  The  change  in  the  longitudinal  velocity  component  in  cross 
sections  1,  2,  3  and  4  makes  it  possible  to  note  precisely  the  same 
characteristic  flow  regions. 


The  suggested  difference  scheme  makes  it  possible  to  make 

calculations  for  flows  with  contact  discontinuity  even  in  the  case 
when  the  velocity  value  during  passage  through  it  changes  from  super¬ 
sonic  to  subsonic.  Figure  11  shows  an  example  of  such  calculation. 

In  Fig.  11,  the  dot-dash  line  shows  contact  discontinuity  and  the 
velocity  value  during  transit  through  it  changes  six  times 

(M„*  3,  M_7  0&,  r  m  M> 

The  solid  line  depicts  the  shock  wave  and  the  dash  line  shows  the 
sonic  line. 
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Calculation  accuracy  is  controlled  by  checking  the  Bernoulli 
relationship,  which  is  an  integral  of  the  differential,  but  not 
the  difference  equations.  The  value  of  the  Bernoulli  integral 

for  a  homogeneous  flow  differed  from  the  value  calculated  accord¬ 
ing  to  the  parameters  of  the  free-stream  flow  by  no  more  than  0.5$. 
On  the  contact  discontinuity,  the  value  of  Bernoulli  integral  con¬ 
stants  changes  in  jump  fashion,  according  to  whose  '’blurring"  value 
one  can  judge  the  quality  of  function  of  the  difference  scheme. 

When  constructing  the  difference  system  in  the  work,  contact 
separations  in  all  regions  of  flow  (Pigs.  1,  2,  3>  *0  were  "blurred” 
on  no  more  than  4-5  points,  which  makes  it  possible  to  determine  the 
discontinuity  line,  not  only  by  plotting  the  streamline,  but  also 
according  to  values  of  the  Bernoulli  integral  (Pig.  10). 


The  authors  consider  it  their  duty  to  thank  V.  3.  Balakin  for 
his  useful  discussions  and  the  attention  he  gave  to  the  work. 
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INVESTIGATION  OF  SUPERSONIC  FLOW  AROUND  A  SPHERE  BY 
A  FLOW  OF  AIR  WITH  A  HIGH  STATIC  TEMPERATURE 


Ye.  G.  Shapiro 


The  Investigation  of  supersonic  flows  near  blunt  bodies  has 
been  made  in  many  works  [1],  Approximate  laws  of  similarity  have 
been  established  for  such  flows,  and  the  necessity  of  taking  the 
physico-chemical  processes  in  the  boundary  layer  into  account  has 
been  shown  for  hypersonic  flows. 


It  has  been  assumed  in  most  of  these  works  that  the  gas  ahead 
of  the  wave  has  a  temperature  T,,,  ■  200  -  300°  K,  and  primary  empha¬ 
sis  has  been  placed  on  investigating  the  hypersonic  flows.  However, 
in  certain  cases  associated  with  experiments,  one  has  had  to  deal 
with  a  flow  of  gas  with  a  high  static  temperature  and  a  low  Mach 
number  [Machao<  6  ]  even  at  high  velocities.  This  work  is  devoted 
to  the  numerical  Investigation  of  such  flows. 


Statement  of  the  problem.  Supersonic  flow  of  air  past  a 
sphere  is  examined.  The  equations  of  motion 

,(9r)w  *f"wfO  , 

,  •  Cl) 

f**(**2  **)•»  ,  *-*WJ , 

are  written  in  the  spherical  system  of  coordinates  whose  center 
coincides  with  the  center  of  the  sphere  and  an  axis  running  along 
the  flow  axis.  As  sought  functions,  we  choose  the  radial  function 
u  and  the  circular  function  v  of  the  velocity  component,  as  well 
as  pressure  and  temperature  p,  T. 


las 


The  non-dimensional  variables  are  Inserted  according  to  formu- 

i 

■*•««.«  .  »*••»«,*  ,  r  . 

r'*VT.T’*V  •  C2) 
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As  usual,  the  non-flow  condition  on 
the  body  and  the  laws  of  conservation  on 
the  shock  wave  were  taken  as  the  boundary 
conditions.  Functions  h  and  p  were  cal¬ 
culated  according  to  the  A.  N.  Krayko 
approximation  formulas  [2]  not  only  in 
the  boundary  layer,  but  also  ahead  of  the 
shock  wave  front,  inasmuch  as  T*,  -  2000°  K 
in  most  calculations,  and  at  such  tempera¬ 
tures  as  is  known,  physico-chemical  proc¬ 
esses  already  occur  in  the  air. 


Fig.  1 

1  -  sonic  line 


The  problem  was  solved  by  the  G.  F.  Telenin  method, 
following  was  taken  for  the  independent  variable 


where  f;  is  the  shock  wave  contour,  -4  —  contour  of  the  body. 

The  gasdynamic  parameters  were  approximated  by  Lagrangian  poly¬ 
nomials  according  to  a  five-ray  method  (-■>.  -o;  t.  •  ■,  . 

■  O.Tua  i  '  •  0,300).  . 

2.  Relationship  between  the  position  of  the  shock  wave  and 
the  Mach  number,  temperature,  and  pressure  of  the  free-stream  flow. 

The  calculations  that  were  made  /«.-  i  an*.  ■,  Tm -  x*x>‘ k  ( 
showed  that  th.e  shock  wave  monotonically  approaches  the  body  in 
precisely  the  same  way  as  in  the  case  of  standard  external  condi¬ 
tions  when  the  Mach  number  increases,  while  the  sonic  line  is  dis¬ 
placed  toward  the  axis  of  symmetry  (Figs.  1  and  2). 


Already  in  the  first  works  devoted  to  hypersonic  flows  around 
blunt  bodies,  it  was  shown  that  the  gas  behind  the  shock  wave  cannot 
be  considered  a  perfect  gas.  At  high  Mach  numbers,  temperature  be¬ 
hind  the  shock  wave  becomes  sufficient  for  physico-chemical  processes 
to  develop  in  the  shock  layer,  in  connection  with  which  regression 


I 


of  the  shock  wave  e  from  M,,,  is  more  complex  in  comparison  with  the 
case  of  a  perfect  gas.  The  position  of  the  shock  wave  in  the  per¬ 
fect  gas  is  determined  by  and  by  the  heat  capacity  relationship 
y  and  does  not  depend  on  pressure  and  temperature  ahead  of  its  front. 
Inasmuch  as  the  increase  in  temperature  in  the  equilibrium  gas  is 
accompanied  by  physico-chemical  processes,  the  increase  in  tempera¬ 
ture  ahead  of  ths  shock  wave  front  with  a  fixed  M*.  causes  physico¬ 
chemical  processes  in  the  shock  layer,  and  the  position  of  the  shock 
wave  will  depend  upon  To,,,  monotonically  (Fig.  3).  For  an 

illustration,  we  shall  examine  the  process  of  dissociation  of 
and  O2  molecules  in  the  shock  layer.  Figure  4,  a,b  shows  concen¬ 
trations  with  molecules  of  Og  (curve  1)  and  Ng  (curve  2)  related 
to  standard  atmospheric  conditions  taken  from  tables  [3]  according 
to  conditions  behind  the  direct  shock  wave  for  =  2  and  4,  de¬ 
pending  upon  T,,,.  Here,  too,  the  dash  line  plots  the  corresponding 
values  eQ  (T<*>).  The  initial  segment  of  the  curve  Eg  (Tc)  when 
*  2  corresponds  to  the  beginning  of  the  process  of  oxygen  dis¬ 
sociation;  for  M.  *  4  with  precisely  the  same  T*,  we  have  an  al¬ 
ready  developed  process  of  02  dissociation.  When  the  02  dissocia¬ 
tion  process  culminates,  regression  of  the  shock  wave  on  the  £g 
axis  acquires  its  maximum  value  in  both  cases.  The  process  of 


dissociation  is  still  not  complete  in  both  cases,  and  therefore 
there  is  no  second  maximum  on  the  e^CT*)  curves.  Higher  tem¬ 
peratures  in  the  shock  layer  correspond 
to  high  Mach  numbers  ahead  of  the  shock 
wave  with  identical  external  conditions, 
and  therefore  the  higher  degrees  of 
molecular  dissociation  and  the  associ¬ 
ated  maxima  and  minima  of  curves  ^(Ta,) 
are  displaced  toward  the  low  tempera¬ 
tures  when  the  Mach  number  increases 
(Pig.  3).  On  the  other  hand,  the  in¬ 
crease  in  Moo  is  accompanied  by  an  in¬ 
crease  in  pressure  behind  the  wave, 
while  the  pressure  increase  reduces  the 
degree  of  dissociation  at  precisely  the 
same  temperatures.  Therefore,  the  curves 
in  Pig.  3  become  smoother  at  the  high 
Mach  numbers.  Pig*  ** 

Pressure  ahead  of  the  shock  wave  front  has  a  significant  effect 
on  the  position  of  the  shock  wave  in  the  equilibrium  gas.  Figure  5 
shows  the  relationship  between  regression  of  the  shock  wave  along 
the  five  rays  t  Q.287,  0.538,  0.725,  0.8  when  Moo  *  3,  Too  = 

60QQ°  k  with  pressure  of  the  free-stream  flow.  The  relationship  be¬ 
tween  and  poo  has  precisely  the  same  complex,  nonmonotonic  charac¬ 
ter  as  with  T^.  Here,  the  nonmonotonic  nature  of  the  behavior  of 
e  is  already  not  associated  with  changes  of  temperature  in  the  shock 
layer,  which  are  relatively  small,  but  with  changes  in  pressure  in 
it  and  the  maximum  £ C p^)  is  due  to  culmination  of  the  nitrogen  dis¬ 
sociation  process. 

The  shape  of  the  shock  wave  changes  weakly  with  a  change  in 
p^  and  with  the  change  in  T*.  Only  its  position  relative  to  the 
sphere  changes.  When  the  Mach  number  of  the  free-stream  flow  de¬ 
creases,  the  shock  wave  not  only  regresses  from  the  body,  but 
straightens  (the  radius  of  curvature  of  the  shock  wave  increases  on 


the  axis ) . 

3.  Certain  properties  of  similarity  of 
equilibrium  flows  near  the  sphere.  Plow 
around  a  Body  of  assigned  form  By  a  perfect 
gas  is  totally  determined  by  two  parame¬ 
ters  —  Mo,,  and  y  —  and  equilibrium  flow  is 
determined  by  the  Mach  number,  the  equation 
of  state,  temperature  T*,  and  pressure  pM. 
However,  analysis  of  calculation  results 
showed  [4,5]  that  in  the  case  of  an  equili¬ 
brium  gas,  the  number  of  parameters  can  be 
reduced.  An  approximate  formula  was  given 
in  a  work  I1*] 

t.’O.SM  , 


Pig.  5 

M  =  3,  T  *  6000°  K 
J  *  00 


CM 


reflecting  the  relationship  between  regression  of  the  shock  wave  on 
the  axis  with  compression  in  the  shock  wave  .  Here,  too, 

a  formula  is  given  for  the  dimensionless  velocity  of  flow  in  the 
critical  point  .  These  results  are  extended  to  ellipsoids 

in  a  work  [6],  which  showed  that  the  value  of  the  parameter  •  uni¬ 
formly  determines  flow  rate. 


Approximating  relationships  for  e  .are  suggested  in  a  work  [7], 
as  are  flow  rates  for  blunt  bodies  of  any  shape.  In  all  of  these 
works,  the  number  of  determinative  parameters  is  reduced  to  one. 
However,  as  was  indicated  in  15],  one  can  achieve  greater  accuracy, 
particularly  if  one  is  dealing  not  with  too  high  Mach  numbers,  if 
one  leaves  out  two  parameters  —  the  compression  coefficient  and 
the  Mach  number  of  the  free-stream  flow.  Such  a  refinement  was 
made  in  a  work  [8],  which  gives  the  results  of  calculating  flow 
around  a  sphere  by  a  perfect  gas  in  a  broad  range  of  numbers  Mo,  * 

1.5  -  1000,  y  *  1*05  -  1-6.  On  the  basis  of  analyzing  calculations 
in  [8],  similarity  criteria  were  suggested  as  follows: 
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H'  •  *  .  0,07 / * 


(5) 


v,  •■}}*(*  ‘ip,)  ■  t/*m 
A/j  -  *  •  0,r/rtm 

ft't  ‘( t  •  m)tui/a  {<•  ct tut  yjT ]  J 


(6) 

(7) 

(8) 


which  make  it  possible  to  present  the  results  in  the  form  of  analy¬ 
tical  relationships: 


f.  *  \  ' 

—  -  —  » J.il  V,  f  •  0.  '5  V.  • 

J.  ■  JS5«  •  KjC.t-0,  iii'f,  ) 


(9) 

CIO) 

Cii) 

C12) 


Here,  r^  and  ^  are  distances  from  the  axis  of  symmetry  to  the  sonic 
points  on  the  surface  of  the  body  and  the  shock  wave.  The  compres¬ 
sion  coefficient  *  enters  into  the  similarity  criterion  as  the 
basic  parameter.  Furthermore,  values  M,,,  and  *■*,/.*-  ulm,  are  used. 

The  authors  of  a  work  [ 8 J  pointed  out  that  the  error  of  the 
formula  for  eQ  did  not  exceed  4.5%  for  r^  *  0.8%,  and  10%  for 
(when  M,,,  2  2  -  6.6%).  Formula  (10)  was  least  successful  for  the 
rate  of  flow  —  the  error  of  the  formula  reached  24%. 

Formula  (9)  at  high  Me  is  quite  similar  to  relationship  (4). 
However,  at  small  Mach  numbers  (9)  gives  more  accurate  results. 

Thus,  for  M®  *  2,  y  -  1.4,  formula  (9)  gives  1.5%  error,  while 
formula  (4)  gives  15%  error. 

A  check  of  the  suitability  of  the  formulas  (5-12)  was  made 
in  this  work  for  the  examined  case  of  flow  around  a  sphere  by 
equilibrium  air  with  a  high  static  temperature.  Figures  6,  7  and 
8  give  the  curves  calculated  according  to  the  formulas  for  eQ ,  r^ 
and  r^i  In  Figures  6-9,  the  following  symbols  are  used:  light 
circles  and  light  triangles  —  M*  *  2,  4,  respectively,  dark  circles 
and  triangles  —  M,,  *  3,  5  (various  values  of  the  similarity  cri- 
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teria  for  identical  were  obtained  as  the 
result  of  different  T^,),  the  rectangles  note 
data  which  correspond  to  the  fractional 
values  of  the  number  1.5  <  M<»<  6  when  pM  = 

1  atm,  Ta,  =  2000°  K. .  It  is  evident  that 
formulas  (9),  (11)  and  (12)  determine  the 
position  of  the  shock  wave  and  the  sonic 
points  in  the  examined  class  of  flows  quite 
well . 


Pig.  5 


Formulas  (6)  and  (10)  produced  the 
greatest  error.  In  addition  to  the  basic 
determinative  parameter  •  ,  they  also  in¬ 

clude  the  Mach  number  of  the  free-stream 
flow  and  the  dimensionless  pressure  behind 
the  shock  wave  .  In  the  litera¬ 

ture  devoted  to  hypersonic  flow  around  blunt 
bodies,  there  are  indications  that  the  veloc 


Fig.  7 


ity  gradient  in  the  critical  point  depends  solely  on  compression  * 
for  an  equilibrium  gas  [** ,  6].  For  the  case  of  flow  past  a  sphere 


by  an  equilibrium  gas,  the  following  relationship  is  suggested  in 

ttl 


K  t  dr 

-  —  *  I  l.b  ■  I,  5  * 


(.13) 


In  both  of  these  works,  velocities  were  related  to  the  maximum  veloc¬ 
ity  umax>  Furthermore,  it  was  noted  that  when  a  perfect  gas  flows 
past  bodies,  the  parameter  «  must  be  replaced  by  the  value  (r-r 


In  a  work  [7],  velocities  were  related  to  u,*  and  a  universal 
formulas  was  suggested  for  dimensionless  flow  rate  which  is  suitable 
for  both  the  equilibrium  and  the  perfect  gas 


*i  <», 


(.1*0 


where  weakly  depends  on  the  shape  of  the  body.  For  the  Heiz 

solution  [9],  y  *  4/3. 


_ • 
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The  necessity  of  replacing 
*  with  the  value  (r~'Vff 
with  the  examined  flow  of  a  per¬ 
fect  gas  is  explained  in  the  fol¬ 
lowing  way.  We  shall  assume  that 
the  velocity  of  flow  is  determined 
by  compression  and  the  relationship 
has  the  form  g  d 

•  (15) 


Figure  8 

where  .  Considering  that  this  relationship  is  conserved  for 

the  perfect  gas,  we  proceed  to  a  different  velocity  scale  u 

max  • 

Inasmuch  as, 

7T  l  ’ 

and  compression  in  the  direct  discontinuity  is 

t (r- 1 ) 

** 

then  for  non-dimensional  flow  rate  we  obtain 


Ail-.a/Zl 

“«««  *  ’  T" 


Thus,  the  selection  of  umax  as  the  velocity  scale  led  to  relationships 
of  the  type  (16)  and  precisely  this  kind  of  relationship  was  obtained 
on  the  basis  of  analyzing  the  calculation  of  flow  around  blunt  bodies 
by  a  perfect  gas.  The  calculation  of  equilibrium  flow  detected  rela¬ 
tionships  of  the  type  (13)  although  u  was  taken  as  the  velocity 

ITlaX 

scale.  The  fact  of  the  matter  is  that  such  flows  with  high  velocities 
corresponded  to  equilibrium  flows  when  the  kinetic  energy  of  the  gas 
was  much  greater  than  the  thermal  energy,  and,  therefore, 

the  relationship  of  the  type  (13)  practically  coincides  with  relation¬ 
ships  (15).  When  the  number  is  not  too  great  or  when  the  free- 

stream  flow  has  a  high  temperature,  the  difference  between  u,  and 

max 

uB  becomes  significant,  and  the  correct  choice  of  the  velocity  scales 
becomes  important.  Great  deviations  of  the  empirical  formula  (10) 
from  the  results  of  calculation,  and  the  presence  of  additional  para¬ 
meters  and  */•.•*—■  in  it,  are  evidently  associated  with  the  unsuc¬ 
cessful  selection  of  the  velocity  scale. 
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The  calculations  made  in  this  work  confirmed  the  assumption 
that  flow  rate  in  the  critical  point  is  a  uniform  function  of  com¬ 
pression.  The  data  whLch  correspond  to  the  entire  investigated 
range  of  Mach  ,  numbers,  temperatures  and  pressure  were  approximated 
well  by  the  relationship  *.  4*  .  rir 

£2  — i — i — I — 1 — |  i  K  n 

r  j  *  %  -  We  note  that  all  formulas  (10),  (13), 

I  (14)  and  (17)  produce  similar  results  with  high 

compression  values  when  *«'  and  t  .  They 

/  |  are  all  described  by  a  relationship  of  the 

j  !  ”  •  type  (15)  where  the  coefficient  a  changes  with- 

«  i  1  3*  *  in  limits  from  1.6  to  1.76. 


The  velocity  gradient  in  the  critical  point  can  be  approximated 
by  the  following  relationship _ 

^4 — 1  •£(%■ 

* - U .j~  (x8! 


4/  44 


Figure  10 

for  an  ellipsoid  of  rotation  with  a  semi-axis  ratio  of  »*-£, 

Figure  10  shows  curves  calculated  according  to  formula  (18),  which 
plot  the  results  of  calculating  the  flow  around  ellipsoids  according 
to  the  data  of  works  [6,10,11]. 

4,  Modeling  equilibrium  flow  of  a  perfect  gas.  A  work  [5]  sug¬ 
gested  modeling  equilibrium  flow  at  the  forward  part  of  a  blunt 
body  by  flow  of  a  perfect  gas  with  an  adiabatic  exponent 


1  •/  «  *  ff  .  *  > 


(19) 


The  possibility  of  such  modeling  for  the  flew  around  a  sphere  by  a  gas 
with  a  high  static  temperature  was  tested  in  this  work. 
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Pip'S  .  il  -  13.  •  -T.  -  23UO*K 
J  .  -  T_  -  4000'  K 

•  -  ^  -  0.3  •  lO'*  «TM 


Figures  11-13  give  the  distributions  of  the  velocities  and 
pressures  for  a  perfect  gas  Ma  »  3,  y  *  1.18  according  to  the  £ 
coordinate  on  different  ray  .  The  symbols  on  these  curves 

plot  points  which  correspond  to  values  of  precisely  the  same  para¬ 
meters  for  flow  around  the  sphere  m„*  3,c-  i  «tv.  t„»  23oo  k,  m.-  3,  p.  -  i  .™,  rm.  -tooo' 
and  p-  “  °-3  •  aooo"  k. 


t 


l 


Compression  is  identical  in  all  the  examined  cases:  «  = 

0.1845.  Differences  in  the  fields  of  velocities,  pressures  and  in 
values  '*■('*)  do  not  exceed  3%»  The  temperature  field  for  the  exam¬ 
ined  cases  (curves  o  -  j  -  ua,  i  -  t.  -  hoo'k,  2  ~Tm  -  -kxo'k,  3  -  0.3  io  1  »tm> 

are  presented  for  the  zero  ray  in  Figure  14  for  the  last  ray  in  Figure 
15.  As  is  evident,  the  temperature  difference  is  more  significant 
and  reaches  14$.  The  fact  of  the  matter  is  that  in  a  real  gas  enthalpy 
will  already  not  be  proportional  to  temperature  and,  therefore,  even 
in  a  case  wher  enthalpy  in  the  real  and  simulated  flows  is  identical, 
the  distribution  of  temperature  in  the  shock  layer  will  generally  be 
different.  In  the  case  examined  here,  which  is  flow  of  high  tempera¬ 
ture  air  around  a  sphere,  gas  ahead  of  the  shock  wave  cannot  be  con¬ 
sidered  a  perfect  gas.  Totalenthalpy  of  the  real  flow 
need  not  be  equal  to  total  entropy  of  the  simulated  flow 
with  identical  Mach  numbers  M^  and  compression  *  but 

Enthalpies  behind  the  shock  wave  in  these  flows 


different  p^ 


and 


«STr* 


T7T 
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will  not  be  equal.  Thus  arises 


'■j  iSM  >&«  * 


the  impossibility  of  directly  simulating  temperature  distribution. 
All  the  same,  despite  this,  one  can  obtain  a  temperature  distribu¬ 
tion  in  the  shock  layer  for  the  real  gas  from  the  model  solution. 


At  the  same  time,  by  knowing  T 


and  Mw,  it  is  easy  to  find  the 


total  enthalpy  of  the  free-:  stream  flow  of  the  real  gas  H*.  Having 
taken  the  velocity  distribution  from  tbs  model  solution,  we  find 
enthalpy  in  the  shock  layer  *<*>r  i  .  Now,  inasmuch 

as  the  distribution  of  pressures  is  also  known  to  us,  it  is  not  diffi¬ 
cult  to  find  temperature. 


Figure  14 
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BOUNDARIES  OP  APPLICABILITY  OF  THE  LAW  OF  BINARY  SIMI¬ 
LARITY  DURING  THE  FLOW  OF  A  DISSOCIATING  GAS  AROUND  BODIES 

by 

T.  P.  Stulov 

When  Investigating  non-equilibrium  flow  around  blunt  bodies, 
the  law  of  binary  similarity  plays  an  important  role  [1].  Accord¬ 
ing  to  the  law  of  binary  similarity,  non-equilibrium  flows  of  a 
given  gas  assigned  condition  v  ,  T^  near  geometrically  similar, 
bodies  are  similar,  if  the  condition,  pm  R„  *  const  is  satisfied 
for  these  flows,  where  pM  is  pressure  in  the  free-stream  flow, 

Rq  is-  the  characteristic  size  of  the  body.  The  applicability  of 
the  binary  similarity  in  the  dissociating  gas  is  limited  due  to 
taking  recombination,  which  proceeds  in  the  trinary  collisions,  into 
account  [1],  Inasmuch  as  the  importance  of  the  recombination  terms 
increases  proportional  to  the  approach  to  equilibrium,  then  it  is 
clear  that  one  should  in  fact  anticipate  great  deviations  from  the 
rule  of  binary  similarity  in  precisely  this  region.  In  the  practi¬ 
cal  regard.  It  is  extremely  important  to  determine  the  region  in  the 
plane  of  variables  (p*,  Rq)  in  which  the  binary  similarity  law  is 
valid  with  the  assigned  accuracy. 

Violation  of  binary  similarity  during  flow  around  a  blunt  body 
primarily  has  an  effect  on  the  distribution  parameters  in  the 
vicinity  of  the  critical  point,  for  in  the  critical  point  Itself 
the  composition  of  gas  always  corresponds  to  local  thermodynamic 
equilibrium.  With  the  increase  in  the  parameters  p_  R_ ,  the  region 
of  quasiuniform  in  the  shock  layer  increases  and  finally  encompasses 
the  entire  shock  layer  with  the  exception  of  a  narrow  relaxation  zone 
behind  the  curvilinear  shock  wave.  One  can  take  the  following  con¬ 
dition  as  the  criterion  for  describing  the  region  of  validity  of  the 
binary  similarity  law  for  the  shock  layer  in  general;  the  law  of 
binary  similarity  is  valid  when  the  length  of  the  relaxation  zone 
behind  the  direct  shock  wave  exceeds  a  certain  proportion  of  the 
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length  of  regression  of  the  shock  wave  ahead  of  the  body  under 
precisely  the  same  conditions,  i.e., 

t.kt.  (1) 

Choice  of  the  arbitrary  parameter  k  insures  satisfying  the  law  of 
binary  similarity  in  the  sought  region  with  the  accuracy  assigned 
in  advance.  Obviously,  one  can  only  find  the  value  of  k  empirically, 
on  the  basis  of  numerical  solutions  to  the  problem. 


Equality  1  ■  ke  gives  the  boundary  of  the  region  in  the  plane 
V*  exP1iGit  formula  will  be  obtained  below  for  this  bound¬ 
ary  and  numerical  estimates  will  be  given  for  specific  cases. 


In  order  to  establish  the  relationship  between  regression  of 
the  shockwave  e  and  the  length  of  the  relaxation  zone  behind  the 
direct  shock  1,  we  shall  use  the  approximate  law  of  flow  similarity 
at  the  leading  surface  of  the  body  examined  [2],  Formula  (5)  of 


this  work 


7  ■ 7  , 


(2) 


establishes  the  relationship  between  e  and  the  distribution  of  den¬ 
sity  behind  the  direct  shock  wave  p(x)  under  precisely  the  same 
conditions.  Formula  (2)  was  obtained  with  the  following  assumptions: 
1)  the  flow  of  the  mass  of  gas  changes  along  the  critical  stream 
line  according  to  a  linear  law;  2)  the  ratio  of  regression  of  the 
shock  wave  to  the  radius  of  curvature  of  the  body  RQ  is  proportional 
(with  a  coefficient  A)  to  the  integral  compression  of  gas  on  the 
flow  axis. 


In  order  to  determine  the  density  profile  p(x),  we  obtain  an 
approximate  solution  to  the  problem  of  flow  behind  the  front  of  the 
direct  shock  wave  in  the  dissociating  gas.  We  shall  examine  the 
relaxation  zone  behind  the  direct  shock  wave  in  a  mixture  of  gases 
M2  and  X  formed  by  the  dissociation  of  component  M2  into  atoms. 

The  equation  for  the  change  of  the  atomic  component  is  written  as 
follows : 


f»7Ti2?  *• 


1  ' 


(3) 


66 


where 


*4'*. (T)lkx (*t)r<nt) ♦  *,  W/ (*) * *, Wjffll .  ■ 


By  excluding  the  value  yCMj)  from  the  brackets  in  the  right  hand 
part  of  equation  (3),  with  the  aid  of  the  material  balance  equation 


(4) 


and  density  p  with  the  aid  of  the  equation  of.  state,  we  obtain  the 
equation  for  determining  y(M) 

V  ■  W  -6W  . 


(5) 


where 


I  -m(*)r(*)_  Kjr).ftTK(T)  . 

2m(M)  0  *  ' 

Thus,  if  one  ignores  the  relationship  between  the  multiplier  before 
the  parenthesis  in  the  right  hand  part  of  (5)  and  y(M),  then  the 
rate  of  dissociation  becomes  a  quadratic  function  of  the  concentra¬ 
tion  of  the  atomic  component.  We  designate 

*  zp**4  * 

«.*T iu)  •  * Tjn  (g 

In  order  to  solve  equation  (5)  with  the  initial  condition 


when 


(7) 


we  shall  consider  that  the  coefficient  t,  and  a^  in  the  relaxation 
zone  are  constant.  Below,  when  making  the  calculations,  we  shall  use 
certain  average  values  of  these  components.  The  solution  to  equation 
( 5  J  has  the  form 


(8) 

Solution  (8)  takes  a  particularly  simple  form  in  the  case  of  disso¬ 
ciation  of  a  pure  gas,  i.e,,  when  y(X)  ■  0 


It  follows  from  the  solution  (8)  that  one  should  view  the  value 
i-C,/ «  as  the  length  of  relaxation  in  the  dissociating  gas,  and  not 

simply  t.  .  In  order  to  obtain  the  explicit  relationship  between 


t-  and  p^,  we  shall  assume  that  the  following  approximate  equalities 
are  valid: 


£ 


Cio) 


by  substituting  (10)  in  the  expression  for  we  obtain 


Thus,  the  length  of  relaxation  in  the  dissociating  gas  is  inversely 
proportional  not  to  paJ  as  with  the  purely  binary  processes,  but 
to  the  value  *-^*«*-  .  The  second  factor  arises  due  to  the  rela¬ 
tionship  between  the  equilibrium  degree  of  dissociation  and  the 
characteristic  pressure.  We  recall  that  temperature  should  be  taken 
in  a  certain  average  point  of  the  relaxation  zone  for  calculating 
a  and  b. 


The  distribution  of  p(x)  necessary  for  substitution  in  equation 
(2)  should  be  found  from  the  algebraic  equations  of  mass,  momentum, 
and  energy  conservation  using  solution  (8).  In  this  case,  the  rela¬ 
tionship  between  and  )  would  be  non-linear, 

as  in  formula  (8).  However,  in  order  to  obtain  the  criterion  curve 
(1),  it  is  more  simple  to  use  the  linear  relationship  between  p(x) 
•*r(- */H  .  Such  a  replacement  somewhat  changes  the  value  of  the 
integral  (2),  but  it  significantly  simplifies  the  calculations. 


Thus,  assuming 


(12) 


and  by  substituting  this  formula  in  (2)  we  obtain  the  relationship 


between  the  length  of  relaxations  1  and  shock  wave  regression  e 


m 


fte  ■ {frf,)jrt  *  A*.J »- 


(13) 


Here  ft  •  ft  -  density  of  gas  at  the  beginning  and  end  of  the  relax¬ 
ation  zone,  respectively.  Values  *,-** and  */ • 
are  equal  to  regression  of  the  shock  wave  in  the  equilibrium  and 
the  frozen  flows.  By  using  these  symbols,  we  obtain 

t 


4 


*/*4 


(14) 


By  excluding  e  from  ( 1 4 )  with  the  aid  of  the  criterion  equality 
we  obtain 

tfit,  ■  t 


f.ki  , 


(15) 


We  shall  insert  the  dimensionless  values  of  shock  wave  regression 

and  **  •*//*«  .  By  substituting  these  formulas  and 

expressions  [11)  for  (  in  formula  (15)  we  obtain  the  boundary 
equation  of  the  region  of  validity  of  the  binary  similarity  law  in 
the  following  form: 


urn j_  ■  _  r  i 


(16) 


where  it  is  taken  into  account  that  the  position  of  the  shock  wave 
in  the  equilibrium  flow  depends  on  pressure  in  the  free-stream  flow, 
although  quite  weakly  [3]. 


If  the  curve  (16)  is  plotted  in  the  Cartesian  coordinate  system 
,  then  the  region  of  applicability  of  the  binary  similarity 
law  will  lie  on  the  left  side  of  the  curve,  i.e.,  in  the  direction 
of  decrease  of  the  derivative  a.*.  during  motion  from  the  curve,  inasmuch 
as  the  condition  is  satisfied  in  this  region.  Inasmuch  as  the 
relationship  is  weak,  then  with  sufficiently  small  p^,  when 

is  small  in  comparison  with  unity,  curve  (16)  is  similar  to 
a  hyperbola  .  With  the  increase  in  p^,,  curve  (16)  inter¬ 

sects  the  hyperbola  with  a  decreasing  constant  c  (the  decrease  in  c 
occurs  according  to  the  law  )  and  when  p«— •  we  have 

i.e.,  Ro  stripes  to  0  more  rapidly  than  for  the 
hyperbola.  Thence  follows  the  important  conclusion  that  with  any 
value  p^  Rq,  the  region  of  validity  of  binary  similarity  law  is 
limited  with  respect  to  the  high  pressure  values  in  the  free-stream 
flow. 
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: _ '..misjar:  ■i'V- 


r  .C* ,  -  ■  ’ 


Examples  of  curves  (16)  for  the  case  of  air 
are  shown  In  Figure  1  in  the  logarthmic  coordi¬ 
nate  system  (RQ  in  CM  ,  p,,  in  atm).  Curves  1 
and  2  relate  to  conditions  in  the  free-stream  flow 
“  13,  *  250°K,  when  the  dissociation  reac¬ 

tion  of  02  predominates.  For  the  curve  1  ,  k  =  1, 
k  *  0.5  for  curve  2.  It  is  evident  that  the  change  in  k  leads  to 
a  small,  nearly  parallel  shift  of  the  curve.  Curve  3  was  plotted 
for  a  case  nm-  a.  t.  -  ao’K  (k  *  1),  when  nitrogen  dissociation  is  the 
predominant  reaction.  Inasmuch  as  the  characteristic  length  of 
relaxation  in  air  monotonically  decreases  with  the  increase  in  tem¬ 
perature,  then  the  boundary  curves  (16)  when  i3*m.*25  should 

be  situated  between  curves  1  and  3  {i) . 


Figure  1 


Curve  (16)  can  be  viewed  in  precisely  the  same  way  as  the  arbi¬ 
trary  boundary  between  the  equilibrium  and  nonequilibrium  flows,  if 
one  terms  a  flow  to  be  equilibrium  flow  in  which  the  length  of  the 
relaxation  zone  does  not  exceed  a  certain  small  proportion  of  shock 
wave  regression.  Inasmuch  as  ,  then  with  small  k  equation  (16) 

can  be  somewhat  simplified  and  rewritten  in  the  following  form: 


t(r > 


(17) 
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CONDITIONS  ON  THE  INCIDENT  AND  REFLECTED 
SHOCK  WAVES  IN  AN  EQUILIBRIUM  GAS 


V.  N.  Mirskiy 


Basic  Equations 

An  equilibrium  flow  of  a  gas  mixture  in  shock  tubes  behind 
the  incident  and  reflected  shock  waves  is  examined.  The  flow  is 
described  by  a  system  of  equations  which  express  the  laws  of  con¬ 
servation  of  the  mass  flow,  momentum,  total  energy,  and  the  equa¬ 
tions  of  state.  We  shall  write  the  equations  form  a  plane,  incident 
wave  in  the  system  of  coordinates  associated  with  the  wave. 

2  •  ■  P.  I” 

r*  « 

A——,  (1) 

a". 

4.  _r 

V.1  *  *  '•*  *  S*  *  2 

.»,  •  r.: 

where  -  gas  parameters  ahead  of  the  wave,  t 

-  parameters  behind  the  incident  wave  subject  to  determination. 

In  order  to  establish  gas  parameters  behind  the  reflected  shock 
wave  and  to  calculate  its  velocity  D2,  we  proceed  to  a  coordinate 
system  which  moves  with  the  reflected  wave  and  shill  write  the  laws 
of  conservation  for  it,  taking  into  account  that  the  gas  behind  the 
reflected  wave  should  be  at  rest  relative  to  the  tube  end.  We  immediately 
note  that  gas  velocity  ahead  of  the  reflected  wave  will  be  an  unknown 
value  relative  to  the  incident  wave,  unlike  the  equation  of  the 
incident  wave.  The  equations  have  the  following  form: 
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h  •  *  *  p  : 

<  S  ■  Kl- 


»  •  <  '  ^  ' 
.  I  •  ‘  t  *  «' 


(2) 


where  D2  is  the  velocity  of  gas  behind  the  reflected  shock  wave, 

•  v  —velocity  ahead  of  the  reflected  shock  wave  relative  to 
the  sides  of  the  tube,  a  value  known  from  solving  problem  (1), 

(V^  +  D2) — gas  velocity  ahead  of  the  reflected  wave  in  the  labor¬ 
atory  coordinate  system. 

We  obtain  the  following  equation  system  after  simple  transfor¬ 
mations  of  equation  (2): 


uVa.»* 


f  / 


^  .  t  (?»/?.)' 
v,1  i  (?>/•■’)' 


0,  m  t 


(2M 


Equations  (1)  and  ( 2* )  can  be  written  in  the  universal  form 
-with  a  certain  parameter  i: 

fLu_  *  JL^ZL-i .  -£i_  . 

A. . .  ?  ?  \  ./A)’ 


P~  "  u*  i 

n  .A  (P  ,T  ) 

_ ' 


(3) 
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When  i  *  0,  we  obtain  equations  of  the  incident  shock  wave  (1), 
when  i  *  1 — reflected  wave  (2*). 


For  a  perfect  gas  f""**  from  (2')  or  (3)  when  i  =  1,  one 
can  obtain  a  gas  state  behind  the  reflected  shock  wave  if  the 
pressure  ratio  p. on  the  front  of  the  incident  wave  is  known. 
This  can  be  determined  from  the  system  of  equations  (1),  using  the 
measured  value  of  the  Mach  number  MQ  in  the  form  (1): 


*  jiij •  -  (3r-o}pr')M!  ■  :2l 


(4) 

(5) 


The  ratio  of  the  velocity  of  the  reflected  wave  front  to  the  velocity 
of  the  free-stream  flow  can  be  determined  in  the  form  [1] 


D,  t 

•7  '  cl  7 

•  r5  ~K 


or  by  replacing 


2  '? 


with  its  value  via  and  y,  we  obtain: 


h  >hK'frV 

V.  *  (f* 


(6) 


CALCULATION  OF  THE  EQUILIBRIUM  GAS  PARAMETERS 

A  program  of  calculating  chemical  equilibrium  of  an  arbitrary, 
multicomponent  gas  system  was  used  to  find  h  and  p  in  the  system  (3). 
In  order  to  calculate  the  integral  and  differential  thermodynamic 
characteristics  c,  . e,.  j.  af  *  it  is  necessary  to  know 

the  equilibrium  gas  composition.  Finding  the  equilibrium  composition 
of  a  chemically  reacting  gas  at  assigned  temperature  and  pressure 
pertains  to  solvine  a  system  of  nonlinear,  algebraic  equations  in 
which  the  following  are  mathematically  written:  The  law  of  conser¬ 
vation  of  matter,  Dalton's  law  and  the  law  of  moving  masses.  An  equa 
tion  of  neutrality  of  the  gas  mixture  is  inserted  into  the  general 
system,  if  ionized  components  are  considered  in  the  reaction  products 
The  ionized  reaction  products  are  included  in  the  material  balance 


equation,  the  Dalton  equation,  and  increase  the  number  of  equations 
of  the  law  of  acting  masses.  The  nonlinear,  algebraic  equation  system 
is  solved  by  the  Newton  method.  The  linear  equation  system  obtained 
in  this  case  for  the  logarithmic  correction  is  solved  by  the  Gaussian 
method  with  selection  of  the  main  component. 

Practical  work  with  the  program  includes  the  following  stages. 

The  number  K  is  determined  on  the  basis  of  the  initial  chemical  compo¬ 
sition  of  the  gas  system.  One  then  determines  the  number  i,  where  K 
is  the  number  of  elements,  including  the  electrons,  while  i  is  the 
number  of  components.  A  strict  order  of  gas  system  components  is  pre¬ 
liminarily  established,  and  then  all  atomic  components  and  electrons 
are  arranged  at  the  beginning  of  the  list.  A  matrix  of  stoichio¬ 
metric  coefficients  m**  is  constructed  according  to  the  selected 
i  and  K  numbers.  Coefficients  of  polynomials  +T  for  the  thermo¬ 
dynamic  potentials,  the  level  of  determining  enthalpy  •*<  ,  the 

material  constants  c.-Zlkum  >  and  the  initial  approximation  &*.*.  , 

. are  written  in  the  point  from  which  the  calcula¬ 
tion  begins.  «.*  -«*•  —total  number  of  moles,  while  < — number 

of  moles  of  the  base  components. 

The  program  provides  for  orientation  to  any  set  of  the  K  base 
components  of  all  sets  i.  The  transition  in  the  calculation  process 
to  orientation  according  to  the  predominant  components  provides  a 
gain  during  calculations  and  expands  the  region  of  convergence  of 
the  Newton  method. 

When  making  the  series  of  calculations  in  a  certain  range  of 
temperature,  it  is  more  convenient  to  begin  the  calculation  with 
the  higher  temperatures,  inasmuch  as  one  can  take  the  composition 
of  the  fully  dissociated  mixture  as  the  initial  approximation  in  this 
case. 


The  final  results  subject  to  print  out  are  arranged,  beginning 
with  a  cell  ,  where  *  •  dooo),  .  The  following 

values  are  printed  out  r .  *  *  . 

c, ,  r • 
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where  T  and  P  are  temperature  and  pressure  at  which  the  thermodynamic 
values  and  the  mixture's  composition  are  calculated,  Nq — total  num¬ 
ber  of  moles,  N. — number  of  moles  of  the  components  i,  s — entropy, 
of  mixture,  h — enthalpy,  u — internal  energy,  and  Cv — heat  capa¬ 
cities  at  constant  pressure  and  volume,  respectively,  r-c>  c. 
adiabatic  exponent,  y — molecular  weight  of  mixture,  p — density. 

The  last  occupied  MOZU  cell  (p)  has  the  address  * 


This  program  can  be  used  as  a  subprogram,  for  example,  for 
solving  the  gas  dynamic  problem. 

The  maximum  possible  values  of  i  and  K  can  be  determined  from 
the  approximate  relationship: 


la 

J  A 


H 

Vi 

■"si 

li 

r- 


or 


5! 2  .  it  K‘  ti  -I  *  S<  *  '«  •  iO?: 


lx  *  :a  *  sk  «  2S4j 


Assuming  that  K  =*  10,  we  obtain  *  *  km  ,  etc. 


A 


One  can  familiarize  one's  self  in  detail  with  the  mathematical 
process  and  the  method  of  solution  [2], 

Calculation  time  on  the  triple  address  computer  with  an  oper¬ 
ational  memory  volume  of  ( 4096 ) 10  cells  of  1  point  (T,  P)  ranges 
from  a  few  seconds  to  1-1/2  minutes,  depending  on  K  and  i,  as  well 
as  on  the  method  of  assigning  the  initial  approximations  during 
transition  from  one  point  to  another.  The  simplest  method  was 
chosen  in  the  given  version:  The  solution  in  the  preceding  point 
is  taken  as  the  initial  approximation  in  the  next  point. 

THE  NUMERICAL  SOLUTION 


Numerical  calculations  of  the  thermodynamic  and  gas  dynamic 
parameters  behind  the  plane  incident  and  reflected  shock  waves  in  a 


t 
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mixture  of  gases  consisting  of  Ar  and  K^O  were  made  as  an  example. 

The  system  of  equations  (3)  was  solved  by  the  Newton  method  with 
the  following  initial  conditions:  Pq  *  0.0065  atm,  Tq  *  293.15°K; 
a  range  of  the  velocities  of  the  incident  wave  from  1500  m/sec  to 
5000  m/sec  with  a  step  AUq  *  250  m/sec  was  examined.  The  program 
took  into  account  the  following  components:  h  ,o  .  At .o, .  .  e*  . 

Calculations  were  made  with  various  initial  molar  concentrations 
of  water  in  the  mixture:  -  io'*;  o,ooi  ;  o,oi  ;  o.os  ;  o.i  ;  0,2  *  0,3  ;  0,5. 

Figures  1-8  give  the  results  of  calculating  the  gas  dynamic  para¬ 
meters  behind  the  incident  and  reflected  shock  wave  depending  on  the 
velocity  of  free-stream  flow  and  the  molar  concentrations  of  the 
gas  in  the  cold  state.  The  dash  curve  in  all  the  graphs  is  a  cal¬ 
culation  of  the  mixture  of  gases  consisting  of  tixo-io-'  and 

1  -  io"’  in  the  cold  state,  i.e.,  an  "almost"  perfect  gas. 

It  follows  from  an  examination  of  the  thermodynamic  potential  of 
argon 

<p‘(At)  •  49,4668  *  4.968Zh  T 

that  the  enthalpy  of  argon  acquires  the  form: 

where  Cp  -  4.9682.  It  follows  from  the  relationship  that 

y  *  1.66. 

According  to  the  theory  of  the  reflected  wave,  gas  pressure  and 
temperature  increase  approximately  as  the  square  of  velocity,  while 
density  tends  toward  a  finite  value  in  a  perfect  gas,  with  the 
increase  in  shock  wave  velocities  (formulas  (4)  and  (5)). 

The  velocity  of  the  reflected  shock  wave  was  calculated  accord¬ 
ing  to  formula  (6)  for  certain  values  of  velocity  of  the  free-stream 
flow  and  a  comparison  was  made  with  the  numerical  results.  The 
deviation  does  not  exceed  4  m/sec. 

Figures  1  and  2  give  the  relationship  between  the  velocity  of 
gas  behind  the  incident  and  reflected  shock  waves  in  the  laboratory 
coordinate  system  and  the  velocity  of  the  free-stream  flow. 
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Figure  1 

Key:  (1)  1  M/sec 


Figure  2 

Key:  (1)  1  M/sec 


/ 

V 


On  each  curve,  the  molar  concentration  of  water  in  the  cold  state 
which  corresponds  to  this  version  is  indicated.  The  great  differ¬ 
ence  between  the  dash  curve  and  the  solid  lines  is  due  to  the  trans¬ 
port  of  part  of  the  kinetic  energy  of  the  gas  to  vibratory  energy 
of  the  molecules  in  the  real  gas,  by  energy  expenditures  on  the 
dissociation  of  molecules  into  atoms  and  electron  excitation  of  atoms. 
The  given  calculation  did  not  take  into  account  the  ionization  of 
atoms  and  molecules,  but  this  effect  does  not  introduce  additional 
difficulties  into  the  function  of  the  program  and  can  easily  be 
taken  into  account  if  desired. 

Velocity  behind  the  inciden*  and  reflected  shock  waves  begins 

to  diminish  with  the  increase  in  the  velocity  of  the  free-stream 

flow  V  from  a  value  V„  =  1500  M/sec.  This  decrease  is  associated 
o  o 

with  excitation  of  the  internal  degrees  of  freedom  of  the  molecules 
•and  dissociation  of  molecules  into  atoms.  It  is  interesting  to  note 
one  characteristic  of  this  problem.  In  those  points  in  which  velocity 
behind  the  incident  and  the  reflected  shock  .waves  acquires  its  least 
value,  the  concentrations  of  0,  and  OH  reach  their  maximum  value. 
With  the  subsequent  increase  in  velocity  of  the  free-stream  flow 
-hese  components  begin  to  dissociate.  This  is  because  the  equilibrium 
” stints  of  #2*  °2  are  very  similar  at  these  temperatures. 

'i  -'.at Ion  practically  terminates  with  the  increase  in  velocity  of 
*  lent  wave,  and  the  gas  behaves  as  a  perfect  gas.  At  high 
-«  f  the  free-stream  flow,  the  electron  excitation  of  atoms 
s  -r. tails  a  loss  of  kinetic  energy  of  the  gas,  and  as  a 
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consequence,  a  decrease  in  the  growth  of  gas  velocity  behind  the 
incident  wave,  and  even  a  slight  drop  in  velocity  of  the  reflected 
shock  wave. 


Figures  3  and  4  give  the  distribution  graphs  of  density  behind 
the  incident  and  reflected  waves  depending  on  the  velocity  of  the 
free-stream  flow.  As  should  have  been  expected,  the  ?. !? .  ratio 
for  the  "nearly"  perfect  gas  tends  toward  4  with  the  increase  in 
velocity  of  the  free-stream  flow,  inasmuch  as,  according  to  the  for¬ 
mula  j> »  if— —  when  m, — -  and  —  UlL. .  * 

■■•ir.tV .  fr-n 


Figure  3  ,  G 

Key:  (1)  1  m/sec2^ 

(2)  (M)sec 


Figures  5  and  6  show  the  relationship  between  temperature  and 
the  velocity  of  the  free-stream  flow.  At  low  velocities  of  the  free- 
stream  flow,  only  the  advancing  and  rotational  degrees  of  freedom  of 
the  gas  molecules  are  excited.  In  this  case,  as  in  the  case  of  the 
ideal  gas,  temperature  is  the  subordinate  to  an  approximately 
parabolic  relationship  with  shock  wave  velocity.  Heat  capacity 
of  the  gas  increases  with  the  Increase  in  the  velocity  of  the  incident 
wave  as  the  result  of  the  excitation  of  the  vibratory  degrees  of  the 
molecules.  Therefore,  its  temperature  decreases  more  slowly  than  the 
temperature  of  an  ideal  gas.  This  slowing  becomes  particularly 
noticeable  when  thermal  energy  is  expended  on  dissociation. 
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Figure  7 
Key :  1_M_  2 

sec 


atm 


Figure  8 
Key:  1  M/sec 

2  atm 


Figures  7  and  8  depict  graphs  of  pressure  behind  the  Incident 
and  reflected  shock  wave  depending  on  velocity  of  the  free-stream 
flow. 
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INVESTIGATION  OF  A  SPATIAL,  SUPERSONIC  FLOW  AROUND  THE 
LEADING  PART  OF  BLUNT  BODIES  BY  THE  ESTABLISHMENT  METHOD 


by 

Yu.  M.  Lipnitskiy,  Yu.  Ya.  Mikhaylov  and  K.  G.  Savinov 

In  order  to  solve  the  problem  of  supersonic  flow  around  blunt 
bodies,  the  establishment  method  was  used.  This  method  is  based  upon 
a  difference  device  [1,2]  in  combination  with  the  trigonometric 
approximation.  The  results  of  calculations  of  triaxial  ellipsoids 
are  given  as  an  example.  Primary  attention  was  paid  to  investigating 
flows  without  a  symmetry  plane.  The  described  algorithm  is  a  modifi¬ 
cation  of  methods  presented  in  works  [3,5] . 


We  shall  take  a  system  of  gas  dynamic  equations  in  the  spherical 
coordinate  system  »•  *  ,  which  describes  the  unsteady  three-dimen¬ 
sional  flows,  in  the  following  characteristic  form: 


A 


0 

j/c  •  1  0  0  0 

t 

j/C  1  0  0  0 

V 

A- 

0  0  1  0  0 

w 

0  0  0  r  0 

* 

0  0  0  0  1 

Here  t'  is  time;  £ « l (t. *.•.+)  — a  function  which  satisfies  con¬ 
ditions  *.+)  and 

equations  of  the  surfaces  of  the  body  and  the  shock  wave,  y — 
relationship  of  specific  heat  capacities,  c— speed  of  sound, 

is  the  projection  of 

velocity  vector  on  the  normal  to  the  surface 

and  V  and  W  are  projections  on  2  linearly  independent  vectors  of 
the  orthogonals  <7. <*«,-.<,/»)/:■«  .  w*f* * 

-components  of  the  velocity  vector  along  axes  (*,«,?>.  A— 
diagonal  matrix  with  elements  of  the  type 
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(2) 


namely,  a,  XO.Ah-a.-^-aW  . 

In  addition  to  the  sought  functions  and  their  derivatives  according 
to  0  and  <p ,  and  *«#  also  enter  into  the  composition 

of  the  vector  of  the  right  hand  part  of  /  . 

The  solution  is  sought  in  the  region  hounded  by  the  surface  of 
the  body,  the  shock  wave,  and  the  conical  surface  .  This  sur¬ 

face  is  selected  such  that  the  velocity  component  normal  to  it  is 
greater  than  the  speed  of  sound.  One  should  impose  the  following 
boundary  conditions  for  unity  of  solution  to  the  examined  problem. 

On  the  body  condition, 

(3) 

should  be  satisfied  and  ($•”  should  be  satisfied  on  the  shock  wave- 
conditions  of  conservation  of  flows  of  mass,  energy,  and  momentum 


(4) 

(5) 

?'  ■  3)' 

(V- J))‘*  2H-(um~D)\ZH  f  H  *  H( p,  f) 

(6) 

V  •  V  ,  W  •  W' 

w  40  • 

(7) 

Here  the  subscript  °°  designates  parameters  ahead  of  the  shock 
wave  while  v  designates  the  relative  velocity  of  displacement  of 

the  shock  wave  in  the  normal  direction. 

It  is  assumed  that  t  ■  0  assigns  the  position  of  the  shock  wave 
and  the  distribution  of  all  gas  dynamic  parameters  in  the  region 
where  the  solution  is  sought  at  an  initial  moment  of  time.  The  steady 
state  solution  is  obtained  as  the  unsteady  limit,  when  t  -*■  «  with 
boundary  conditions  of  the  problem'  that  do  not  change  in  time  (para¬ 
meters  of  the  free-stream  flow  and  the  conditions  of  non-flow  on  the 
surface  of  the  body)*  It  is  subsequently  considered  that  *  0. 
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We  take  the  following  grid  during  construction  of  the  numeri¬ 
cal  method 

*„•'»/*  ,  B,  •  tB'/L  ,  f'-UK/K  ;  n.o.t,..., 
m  •  0,1 . M  ;  t*OJ . L  ;  k- 0,t,...,  K-t  . 

Following  work  [6],  we  used  the  trigonometric  approximation  accord¬ 
ing  to  ip  for  smooth  bodies,  and  for  the  nonsymmetrical  bodies,  in 
order  to  improve  accuracy,  we  introduce  the  radius  of  the  body  into 
the  approximation  formulas.  With  an  even  K,  such  an  approximation 
has  the  form: 


X*G(9,f)[s  Jjf *1  (.tint*] 
1  •  •  1 


The  inclusion  of  the  trigonometric  approximation  makes  it  possible 
to  insure  sufficient  accuracy  without  a  significant  increase  in  the 
grid  nodes  . 


Equations  of  the  system  which  is  obtained  following  the  trigono¬ 
metric  approximation  are  replaced  by  difference  equations  in  the 
following  way:  All  equations  which  correspond  to  the  first  two 
equations  in  system  (1)  are  approximated  in  the  nodes  (»'\ .  «•  j 


,1  ..t 


2*  1  1 |  *  "  . 
jt'7 'K.'yK.,',)- 

.,y"  „x"  -x"  -2x’ 

Jf  '  m.t'f  t  rn.t.k  m»l,rn,k 

■  TpfCt -  CJ 

H  ♦-* 

*’?{*.**  ‘X-u*  *  ' 
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The  remaining  equations  are  approximated  in  the  nodes 


.b(£lf 

sriS/CL 

C’.Mu. -C.J/ 

f/  *  MT  I 

'•r(v.-rv..  ) 


(9) 


In  formulas  (8)  and  (9)  «,-*A  .  vV*,/*,* »>.•£)*  *'/»•'  .->/»*.  »  -  4*. 

And  when  f.r  •*,»«,#,  >«  and  when 

Such  an  approximation  insures  the  second  order  of  accuracy  for  the 
for  the  steady-state  solution. 


-I 


Boundary  equations  are  written  in  the  following  way  in 
difference  form: 
on  the  body 

on  the  shock  wave 


°w  ■« 


the 

(10) 


fm  1  -0) 

(11) 

*  *  •>« 

•7  ^  ,  ”~T  . 

%  ’t.  •*) 

(12) 

t  **  r  *,ir  « 

(Um-D)  >ZHm  .(Vm 

(13) 

... 

m  *  —r  **  • 

(14) 

. w.  ■  k. 

♦  .  «»t  ♦  .  t  ♦ .  /  r 

*■■/•«  T-(Vi  y  "zrtM  % 

v*-[\  '(y)  %  //w. 

/•  •• I//’* *r 

B'F  /l«\ 
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Here  the  coordinates  of  the  shock  waves  and  their  derivitives 
according  to  0  and  ip  are  taken  from  the  previous  iteration. 


Unilateral  differences  are  used  on  the  boundary  of  the  region 
situated  in  the  supersonic  part  of  the  flow  (i.e.,  when  1  *  L)  in 
place  of  the  symmetrical  differences  for  approximating  the  der¬ 
ivations  according  to  0.  Another  deviation  from  the  presented 
approximation  is  the  0*0  line,  where  the  existing  characteristics 
are  preliminarily  revealed  according  to  theL'Hopitairule.  Here,  too, 
the  unilateral  differences  are  used  for  approximating  all  expressions 
which  contain  differentiation  according  to  0. 


After  replacing  the  differential  equations  with  differences, 
the  determination  of  all  sought  values  of  breaks  down  into 

a  series  of  independent  boundary  problems.  For  and 

with  identical  k,  1  indices,  this  is  a  system  of  the  type 


t  o 
0  c/ 


frA  » 

V 

rft1  j 

tn*  0,1 . M 

jointly  with  the  boundary  conditions  (10)-(13). 
are  obtained  for'  the  other  values 


(15) 

The  following  systems 


The  indices  ..j  , 


y.v  ,  or  W,  or  6 

t-ldW,*(»0  in  t*td*,l(-0  (U  «  ,/-« 

w  ■  —  ,  a  ■  -  — . ,  — ■  ,  « 

r-  id *,*(•')  * 

k  and  1  are  omitted  in  systems  (15)  and  (16). 


(16) 


Such  separation  of  the  difference  equations  is  achieved  thanks 
to  the  transition  from  variables  to  variables  u,v,w  . 

Otherwise,  one  must  jointly  solve  the  double  flow  difference  equations 
as,  for  example,  in  a  work  [1]. 

The  system  (15)  is  solved  by  the  trial  and  error  method,  when 
the  direct  course  consists  in  the  successive  construction  of 


relationships 


m.  0,t 


(17) 


When  m  *  0,  the  boundary  condition  (10)  is  taken  as  (17).  All  sub¬ 
sequent  trial  and  error  relationships  are  obtained  by  excluding  am 
from  the  second  equation  (15)  with  the  aid  of  the  preceding  trial  and 
error  relationship  and  the  first  equation  (15).  Relationship  (17) 
when  m  *  M,  Jointly  with  the  boundary  conditions  on  the  shock  wave 
C 11 ) — ( 1 3 )  make  it  possible  to  determine  all  sought  functions  on  the 
shock  wave.  The  determination  of  in  the  other  nodes  takes 

place  on  the  basis  of  the  first  equations  of  the  system  (15)  and 
relationships  (17). 


Inasmuch  as  is  proportional  to  X  (0),  then  because  of  the 

boundary  condition  (10)  we  have  .  Thus,  the  difference 

equation  (16)  when  m  *  0  can  be  interpreted  as  a  boundary  condition. 
Consequently , finding  V,  W,  and  5  after  solving  the  problem  (15) 
pertains  to  solving  boundary  value  problems  for  three-point  scalar 
equations,  whose  solution  methods  are  widely  known. 


The  satisfaction  of  conditions  of  uniform  correctness  of  the 
boundary  value  problems  is  necessary  and  sufficient  for  stability 
of  the  trial  and  error  method.  According  to  a  work  [7],  this  is 
equivalent  to  satisfying  the  inequalities 

(18) 


i(.  >  o 


(19) 


we  note  that  in  the  case  of  using  the  approximation  formulas  (8)  for 
all  five  equations  (1),  the  conditions  of  uniform  correctness  are 

*<"**  (20) 
>(♦»)>« 

(21) 


Conditions  (18)  significantly  expand  the  region  of  stability  of  the 
trial  and  error  method  in  comparison  with  (20)  which  makes  it 
possible  to  select  the  transformation  proceeding  from  the 

condition  of  the  best  approximation  of  these  sought  functions. 
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The  choice  of  the  coefficient  u ^  is  made  taking  into  account 
the  inequality 

(M*ef  jf,  *  >**  •  (22) 

This  condition  is  a  sufficient  condition  of  stability  of  the  employed 
difference  method  for  the  scalar  equation.  When  conditions 

(22)  are  also  necessary. 


Figure  1 

Calculations  were  made  in  the  spherical  coordinate  system 
».  • .  *  ,  whose  center  is  situated  on  the  axis  x  of  the  streamlined 
ellipsoid  (Figure  1).  The  distance  from  the  center  of  the  coordinate 
system  to  the  center  of  the  ellipsoid  is  symbolized  by  xQ.  The 
selected  difference  grid  is  characterized  by  the  following  parameters 
M  =  10,  L  ■  20,  K  *  6.  The  transformation  of  was  taken 

in  the  form 

« e(e,+)]  . 

The  velocity  components  a  *,  relate  to  the  maximum  velocity  V  , 

max 

pressure  ,  density  -  *  p-  .  The  linear  dimensions 

pertain  to  the  value  of  the  semi-axis  of  the  ellipsoid  in  the 
direction  of  the  axis  »►  and  the  ellipsoid  was  characterized  by 
the  values  j.-S.  ,  it.^L  («,  »,•}—  dimensions  of  the  ellipsoid  semi¬ 
axis  in  the  direction  of  axes  «.y,  »;  .  The  Mach  number  of  the  free- 
stream  flow  M^  «  3  for  all  calculations,  and  the  adiabatic  exponent 
Y  -  1.4. 


With  the  goal  of  making  a  comparison  with  the  results  obtained 
by  other  methods,  calculation  of  flow  around  an  ellipsoid  with 
i,  -  a,  /,  -  i  at  an  angle  of  attack  >  -  m>. 

was  made.  Figure  2  shows  the  picture  of  intersection  of  the  shock 
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Figure  2 

Key  1:  work  [8] 


wave  and  the  sonic  surface  by  the  plane  . 

In  this  same  figure,  the  circles  plot  the  values 
for  the  corresponding  calculation  obtained  by  the 
G.  F.  Telenin  method  [8].  As  is  evident,  the  coin¬ 
cidence  of  the  geometrical  pictures  is  very  good. 

The  comparison  of  the  gas  dynamic  values  in  the  shock 
layer  that  was  made  showed  that  the  maximum  difference 
that  exists  on  the  edge  rays  does  not  exceed  4%. 


In  order  to  investigate  the  effect  of  the  non-axisymmetric 
nature  of  the  body  on  the  character  of ‘flow,  flow  was  examined  at  a 
zero  angle  of  attack  on  ellipsoids  with  fixed  6^  *  2  and  6 2  =  1; 

Sj  58  1.25;  <5  =  1.5  (in  these  calculations,  0R  *  1.125;  xq  =  1.5). 

It  Is  interesting  to  trace  how  the  change  in  flow  occurs  in  the 
cross  section  ip  =  0,  inasmuch  as  the  contour  of  the  body  here  is 
identical  for  all  three  ellipsoids.  The  position  of  the  shock  waves 
and  the  sonic  lines  for  this  cross  section  is  shown  in  Figure  3. 
Pressure  along  the  body  according  to  the  coordinate  •*•/•»  is  shown 
In  Figure  4.  Figure  5  and  Figure  6  give  an  idea  of  flow  In  the  plane 
*•{  .  The  distributions  of  velocity,  pressure  and  density  on  the 
axes  plotted  according  to  the  coordinate  £  are  identical  for  all  three 
cases  (figures  7  and  8).  Consequently,  the  conclusion  of  a  work(9) 
that  flow  on  the  axis  Is  solely  determined  by  parameters  of  the  free- 
stream  flow,  which  conclusion  was  drawn  for  axisymmetric  bodies, 
can  be  generalized  to  a  broader  class  of  bodies. 

The  study  of  properties  of  spatial  flow  with  the  absence  of  a 
plane  of  symmetry  in  the  shock  layer  is  of  particular  interest.  Such 
an  investigation  was  made  for  an  ellipsoid  with  parameters  6^  *  3 
and  52  *  2.  With  a  fixed  angle  of  attack  a  =  15°,  the  slip  angle 
changed  from  0°  to  15°.  Figures  9-11  give  the  geometrical  pictures 
of  flow  in  the  plane  f  .  o-j  .  +.  $  •  4* ,  $  r 

for  angles  of  0°,  5°,  10°  and  15°  (respectively,  curves  1,  2,  3,  4). 
The  distribution  of  pressure  along  the  body  according  to  the  coordi¬ 
nate  0  #'{*,«  14)  in  the  indicated  planes  is  shown  in  Figures  12-14. 

When  the  slip  angle  Increases  from  0°  to  15°,  the  greatest  changes 
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in  pressure  distribution  occur  in  the  half-  planes 
and  .  For  example,  when  0  =  38°  in  the  half  plane 

pressure  increases  by  27%,  and  pressure  decreases  by  26%  in  the  half 
plane  .  I  /VI.  iH  I" 
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Analyzing  the  flow  field  in  the  shock  layer,  one  can  conclude 
that  in  the  examined  range  of  change  in  the  angle  x  with  sufficiently 
large  angles  0,  the  relationship  between  the  gas  dynamic  parameters 
with  x  111  the  nodfiS  C,  0,  ^  is  nearly  linear  or  is  nearly  linear. 
The  linear  relationship  is  particularly  accurately  satisfied  in  the 
planes  and  .  Figure  15  shows  the  change  of 

pressure  on  the  body  when  0  *  38°,  as  an  example. 

One  can  Judge  how  significantly  the  character  of  flow  changes 
in  the  shock  layer  with  the  presence  of  a  slip  angle  according  to 
Figures  16  and  17,  which  show  distribution  over  the  body  when  0  *  38°, 
pressure  and  the  circular  velocity  component  w  according  to  the 
coordinate  ij;. 

Angle  *  •  o*.  30*.  uo‘,  so*,  in  figure  18  for  an  ellipsoid  with  c  4  -  1.  J,*  i  • 
when  the  vector  of  the  free-stream  flow  forms  an  angle  <*' -  >d* 
with  the  x  axis  and  is  parallel  to  the  plane  which  passes  through 
the  axis  x  and  stops  with  the  axis  z.  In  order  to  illustrate  dis¬ 
placement  of  the  critical  point,  isobaric  projections  on  the  0  *  90° 
plane  are  given  (the  dash  line  corresponds  to  a  sonic  line,  pressure 

relates  to  the  double  velocity  head). 
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INVESTIGATION  OP  SPATIAL,  NON-EQUILIBRIUM 
VOLUME  AROUND  SEGMENTAL  BODIES 


A.  K.  Burdel'nyy,  V.  B.  Minostsev  and  V.  P. 
Shkadova 


Symbols:  -  velocity  component  vKj 

-  maximum  velocity; 


4  +  • 

mSm 

r  \ 

UtTm 

*/ 

fi/‘ 


-  density; 

-  pressure; 

-  temperature; 

-  enthalpy  of  a  mass  unit ; 

-  ratio  of  frozen  heat  capacity; 

"i/fmt  -  concentration  of  the  small  i-th  component 


-  number  of  moles  of  the  i-th  component  per  unit  of 
volume; 

-  molecular  weight  of  mixture; 

-  internal  energy  of  vibratory  degrees  of  freedom 
of  the  i-th  component  per  1  mole; 


93 


R  -  gas  constant ; 

0  -  energy  function  of  the  effect  of  physico-chemical 

processes ; 

M  -  Mach  number;  12a 

R^  -  bluntness  radius; 

o 

*.«.*  -  generalized  curvilinear  coordinate  system; 

#.*.♦  *  sphereical  coordinate  system; 

*.*.+  -  cylindrical  coordinate  system; 

-  generalized  normed  coordinate  system; 


*<•  *»• 


V** 


(or 


Lame  parameters 


) 


c  -  *  ....  -  —  i 

a 

*  « V t », » \‘j 
5  3  — - £— 

•  Ht '«C'*  ♦  *$,♦ 

W«*  .  \ • 


All  geometrical  values  are  related  to  RQ. 

Indices:  •  -  parameters  in  the  free-stream  flow,  s  -  on  the 
shock  wave  ,  S  -  .on  the  surface  of  the  body. 

The  state  of  the  dense  high  temperature  layer  which  forms  ahead 
of  the  blunt  body  flying  at  a  high  supersonic  velocity  in  the  atmos¬ 
phere  of  the  Earth  depends  on  the  physical  transformations  of  air 
particles  on  the  whole  (dissociation,  ionization),  and  upon  the  excita¬ 
tion  of  internal  states  of  the  molecules.  The  necessity  of  strictly 
taking  these  phenomena  into  account  when  determining  the  gas-dynamic 
characteristics  of  the  flow  significantly  complicates  the  investi¬ 
gation  of  the  problem  of  flow  around  blunt  bodies. 
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The  basic  system  of  equations  of  the  non-equilibrium  flow  con¬ 
sists  of  the  gas  dynamic  and  relaxation,  equations 


*j»l/  « o 

V(vp-a't  v?).  mfp<P  .  Vvj  .  ^ 


(1) 


The  relaxation  equations  in  system  (1)  which  describe  behavior 
of  concentrations  of  the  gas  components  and  the  oscillatory  energies, 
are  significantly  nonlinear.  A  singularity  arises  in  these  with  the 
approximation  of  relaxation  processes  to  equilibrium.  One  should  add 
the  equation  of  state  and  the  caloric  equation  to  the  system  (1) 


p* 


f! 


*  *£  (Tlr, 


(2) 


Processes  of  dissociation,  associative  ionization,  exchange  processes 
and  the  excitation  of  the  oscillatory  degrees  of  freedom  of  the  02 
and  N2  molecules  are  distinguished  by  the  basic  relaxation  processes 
in  air.  Here  one  examines  the  components  o, ,  no ,  o .  v.  «;•  ,n,\'no\  r 

between  which  13  chemical  reactions  take  place,  and  the  energy  of  the 
vibratory  modes  is  1 and  l».. 


The  boundary  conditions  of  the  problem  are  set  on  the  shock  wave 
and  the  surface  of  the  body.  In  a  nonviscous  flow,  the  shock  wave  is 
a  front  of  zero  thickness  on  which  the  one  kind-  Rankine-Huboniot  relation¬ 
ship  is  satisfied.  The  chemical  components  do  not  undergo  a  dis¬ 
continuity,  and  the  internal  fluctuations  of  the  degree  of  molecular 
freedom  are  not  excited,  i.e., 

r„-rim  .  • 0  (3) 

The  position  of  the  shock  wave  in  space  is  unknown  and  is  determined 
in  the  solution  process. 

The  condition  of  non-  flow  is  satisfied  on  the  surface  of  a 

body  whose  shape  is  assigned  by  the  analytical  expression,  i.e., 
the  normal  velocity  component  vanishes. 

As  is  known,  the  system  of  gas  dynamic  equations  in  the  sub¬ 
sonic  region  of  flow  is  a  system  of  the  elliptical  type,  that  is 
hyperbolic  in  the  supersonic  region. 
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Therefore,  in  order  to  calculate  flow  in  each  of  these  regions,  one 
needs  a  corresponding  numerical  schemes  which  will  take  into  account 
the  differing  character  of  propagation  of  small  perturbations  in  the 
flow.  In  this  article,  the  solution  to  the  problem  of  flow  around  the 
body  of  rotation  by  a  supersonic  air  flow  which  is  chemically  active 
behind  the  shock  wave  was  obtained  with  the  use  of  two  such  methods. 


The  calculations  were  made  for  segmental  bodies  for 
whose  spatial  streamlining  a  number  of  results  already  exists. 

A  work  [1]  contains  a  survey  of  investigations  that  have  been  con¬ 
ducted  of  the  spatial  flow  of  a  perfect  gas  near  segmented  bodies. 

In  work  [2],  an  investigation  was  made  of  the  problem  of  flow  around 
the  leading  part  of  these  bodies  by  a  nonequilibrium  gas  at  small 
angles  of  attack.  The  problem  of  non-equal  flow  around  a  segmented 
body  in  general  was  solved  here.  In  this  case,  the  data  on  the 
kinetics  of  the  physico-chemical  processes  were. taken  from  a  work  [3]. 

1.  PLOW  AROUND  A  BLUNT  BODY 


The  solution  to  the  problem  of  nonequilibrium  flow  in  the  sub¬ 
sonic-transonic  region  was  obtained  by  the  G.  P.  Telenin  method. 
Generalization  to  the  case  of  a  nonequilibrium  gas  and  a  detailed 
presentation  of  the  characteristics  of  employing  this  method  relative 
to  spatial  flow  is  given  in  a  work  [2].  The  very  same  work  contains 
the  result  of  flow  around  the  leading  part  at  angles  of  attack  a  from 
0°  to  10°.  In  this  case,  these  results  are  expanded  for  an  angle  of 
attack  from  10°  to  25°  according  to  nonequilibrium  flow. 


A  segmented  body  which  is  circular  in  the  middle  region  is 
examined.  In  the  Cartesian  coordinate  system,  the  contour  of  the  body 
in  the  region  of  rounding  has  the  following  analytical  expression. 


7 * {(* •  fS*(i *(j. "**1  }  *° . 


(JO 


where  f..t,  -  change  in  the  origin  of  the  coordinates  with  respect 

to  the  apex  of  the  segment 

. .  V; 
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Q*  -  half-angle  of  the  segment,  0*  -  angle  of  the  inverse  cone. 

The  position  of  the  origin  of  the  coordinates  with  different  a  is 
selected  depending  upon  the  size  and  position  of  the  subsonic  region 
of  flow  ahead  of  the  body.  For  example,  for  an  angle  of  attack 
-  -  is*  t,  •  o, M83,  t,  -  -  0,022,  ,  and  for  —  7,  e,  -  -o^*io. 

In  this  case,  0^  ■  1,116.  In  this  case,  the  flow  is  supersonic  on 
all  outer  rays  and  throughout  all  half  planes. 

The  geometrical  picture  of  flow  of  the  leading  surface  in  the 
symmetry  plane  when  -  -  is*  and  «-«*  is  shown  in  Figure  1. 
The  free-stream  flow  has  the  following  parameters  ■  20,  T^  »  250°K, 

— u 

P  *  10  atm.  the  bottom  radius.  R*  1  cm.  With  an  increase  in 
00  w  O 

the  angle  of  attack,  the  shock  wave  approaches  the  body  from  the 
windward  side  and  straightens  in  the  region  *  0.  The  dimen¬ 

sions  of  the  subsonic  region  significantly  decrease.  The  body's 
critical  point  approaches  the  edge  of  the  body.  Figure  1  shows  the 
position  of  stagnation  points  with  different  angles  of  attack,  while 
the  functional  relationship  (y  stagnation) * with  a  is  given  in  Figure  2. 


We  shall  examine  the  question  of  the  effect  of  nonequilibrium 
processes  in  the  flow  on  the  flow  gas  dynamic  characteristics  in 
more  detail.  It  was  shown  in  work  [2]  that  the  spatial  flow  effects 
directly  ahead  of  the  body  smooth  as  the  result  of  a  change  in  the 
chemical  composition  of  the  mixture.  A  sublayer  of  gas  of  increased 


density  which  is  homogeneous  with  respect  to  temperature  and  com¬ 
position  forms  ahead  of  the  body.  The  flow  rate  of  the  gas  along 
the  shock  layer  decreases  in  comparison  with  the  case  of  flow  of  a 
perfect  gas 

Figure  2  gives  the  velocity  gradients  and.  £  in  the 

stagnation  point;  the  dot-dash  line  was  obtained  from  a  solution  to 
the  problem  for  a  perfect  gas. 

Pressure  in  the  entire  shock  layer 
at  the  leading  surface  of  the  body  in 
the  case  of  the  nonequilibrium  flow  has 
a  higher  value  than  is  the  case  for  a 
perfect  gas.  A  symmetry  in  its  distri¬ 
bution  along  the  surface  is  basically 
determined  by  the  spatial  character  of 
flow  (Figure  3). 

For  angles  of  attack  ^10°,  spatial  effects  had  practically  no 
effect  on  the  distribution  of  gas  composition  in  space.  The  distri¬ 
bution  of  concentrations  y^  was  nearly  axe symmetrical  throughout 
the  entire  flow  region  (for  example,  y^^  in  figure  5  for  -.3*  ). 

At  large  angles  of  attack,  one  observes  a  noticeable  effect  of 
the  three  dimensionality  of  flow  on  the  thermo-chemical  structure 
of  the  gas  in  the  shock  layer.  All  functions  which  describe  this 
state,  p,  T,  yj^  at  angles  of  attack  a  >  10°  have  clearly  pronounced 
maxima  (minima)  in  the  vicinity  of  the  stagnation  point  and  change 
sharply  along  the  body’s  surface. 

In  Figure  4,  one  can  trace  how  gas  temperature  changes  along 
the  shock  layer  when  a  +  25°.  The  solid  lines  T(0)  pertain  to  the 
plane  CO,  II),  while  the  dash  lines  to  the  plane  /f .  y  J  .  The  asymm¬ 
etry  of  temperature  distribution  decreases  toward  the  body  (£  *  0) 
as  the  result  of  the  contribution  of  energy  of  the  noneauillbrium 
processes  to  the  total  flow  energy,  but  temperature  along  the  surface 
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already  clearly  is  not  constant  for  the  case  a  *  5°  [4]  as  well 


It  is  interesting  to  trace  the  character  of  the  change  in  the 
electron  concentration  (Figure  5)  along  the  body  on  the  symmetry 
line.  A  sharp  drop  in  y,  at  the  high  angles  of  attack  is  a  con¬ 
sequence  of  the  more  rapid  expansion  of  flow  from  the  leeward  side 
and  its  sharper  turn  from  the  winward  side. 

2.  FLOW  AROUND  THE  LATERAL  SURFACE  OF  A  BODY 

The  previous  section  gave  the  results  of  calculating  non- 
equilibrium  flow  in  a  shock  layer  near  the  leading  surface  of  a 
segmented  body.  If  the  surface  0  »  0^  is  a  surface  of  the  spatial 
type,  i.e..  If  the  normal  flow  velocity  on  it  is  greater  than  the 
local  frozen  velocity  of  sound  a^.  everywhere,  then  the  data  on  the 
surface  can  be  used  as  initial  conditions  for  calculating  the  super¬ 
sonic  flow  region. 

For  the  case  of  near-equilibrium  flow  of  the  leading  surface, 
a  work  [5]  suggests  an  approximate  method  of  calculating  the  effect 
of  inperfection  of  the  gas  with  hypersonic  flow  around  segmented  bodies 
which,  as  the  conducted  investigations  showed,  proved  to  be  effective 
for  investigating  hypersonic  flow  around  segmented  bodies  of  short 
length  with  a  large  angle  of  inclination  of  the  trailing  cone 
generatrix  and  at  comparatively  small  angles  of  attack. 

In  this  method,  flow  around  the  leading  surface  is  considered  to  be 
completely  equilibrium  and  is  simulated  by  the  flow  of  a  perfect  gas 
with  an  adiabatic  exponent  determined  from  the  condition  of  equality 


of  the  relationship  of  densities  in  the  direct  shock  wave  in  the 
equilibrium  and  simulated  flows.  The  adiabatic  exponent  is  considered 
to  be  precisely  the  same  for  the  entire  supersonic  flow  region,  except 
for  points  on  the  shock  wave,  and  is,  in  fact,  determined  from  the 
condition  of  freezing  of  the  equilibrium  gas  composition  in  the  region 
of  the  body's  rounding.  A  work  [6]  gives  calculations  of  flow  around 
the  lateral  surface  for  conditions  of  developed  non-equilibrium  in 
the  shock  layer  near  the  lateral  surface  of  a  segmented  body  under 
the  condition  of  constancy  of  gas  composition  along  the  stream  line 
in  the  supersonic  flow  region.  The  composition  of  the  gas  slows  in 
each  surface  point  0*0^,  and  the  frozen  adiabatic  exponent  «,  was 
determined  according  to  the  frozen  composition.  Precisely  the  same 
initial  data  were  used  as  in  this  article.  The  calculation  of  flow 
in  the  supersonic  flow  region  was  made  according  to  the  method  of  a 
work  [7].  In  this  case,  the  adiabatic  exponent  */  is  determined 
on  the  layers  by  the  Lagrange  polynominal  according  to  £  and 
The  numerical  investigations  of  the  spatial  flow  around  bodies  with 
inclusion  of  the  chemical  kinetics  equations  in  the  calculation 
requires  a  significant  increase  in  calculation  time  and  complication 
of  the  program,  in  the  investigations  that  were  made.  Only  one  work  has' 
been  published  up  to  now  [8]  which  is  devoted  to  calculating  the  three- 
dimensional,  nonequilibrium  flow  of  a  gas  in  the  supersonic  region 
of  flow  around  bodies.  This  work  investigates  the  effect  of  non- 
equilibrium  dissociation  on  flow  near  blunt,  inverse  cones  situated 
.in  the  supersonic  oxygen  flow.  The  initial  data  for  calculating  the 
supersonic  region  was  taken  from  the  corresponding  numerical  solution 
for  an  oxygen  flow  in  the  shock  layer  near  a  sphere.  There  are 
presently  no  published  works  on  calculating  the  spatial  flow  of  non¬ 
equilibrium  air  in  the  supersonic  region  of  flow. 

As  was  already  noted,  the  effect  of  the  physico-chemical  processes 
on  the  gas  dynamic  parameters  in  the  supersonic  flow  region  was 
estimated  quite  well  with  the  aid  of  simulating  non-equilibrium  flow 
in  the  shock  layer  near  the  lateral  surface  of  a  body  by  a  flow  of 
gas  with  an  adiabatic  exponent  **  ,  This  property  of  the  total 

effect  of  relaxation  of  the  physico-chemical  processes  on  the  gas 
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dynamic  flow  parameters  can  be  used  when  developing  numerical  schemes 
for  calculating  nonequilibrium  flows  in  the  supersonic  flow  region. 

The  results  of  calculations  cited  in  this  section  were  obtained 
according  to  a  method  in  which  the  gas  dynamic  parameters  and  rates 
of  the  physico-chemical  processes  differ  sequentially  according  to  the 
different  schemes.  In  this  case,  the  precise  calculation  of  the  rates 
of  physico-chemical  processes  is  only  made  in  the  symmetry  plane  of 
flow. 


We  shall  examine  the  structure  of  a  system  of  equations  (1)  which  , 
describe  nonequilibrium  gas  flow.  The  first  three  relationships 
describe  the  behavior  of  gas  dynamic  parameters  of  flow.  The  next 
relationship  describes  the  change  in  the  chemical  composition  of 
gas  and  relaxation  of  the  oscillatory  degrees  of  freedom  of  the  molecular 
components.  In  order  to  calculate  flow  in  the  supersonic  flow  region, 
the  initial  conditions  are  assigned  on  the  surface  of  a  spatial  type 
0  ■  0^  from  a  solution  to  the  problem  of  nonequilibrium  flow  of  the 
leading  surface  of  a  body.  The  boundary  conditions  are  precisely 
the  same  as  for  flow  in  the  subsonic  flow  region. 

As  is  evident  from  the  system  (1),  the  effect  of  a  change  in 
gas  composition  and  the  excitation  of  oscillations  in  the  gas  dynamic 
equations  is  only  manifested  cumulatively  via  the  relationship 
between  the  frozen,  specific  heat  capacities  mt  (or  the  frozen 
velocity  of  sound  af)  and  the  energy  function  of  the  effect  of 
nonequibrial  processes  4. 


In  order  to  solve  the  gas  dynamic  equations  in  the  supersonic 
region  of  flow  with  a  known  distribution  of  ,*/.  and  4,  one  uses  a 
scheme  suggested  by  Yu.  Ya.  Mikhaylov  for  a  perfect  gas.  According 
to  this  scheme,  the  gas  dynamic  equations  are  written  in  in  the 
characteristic  form 
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Equations  (5)  are  written  in  the  common,  curvilinear  coordinate 
system  and  are  obtained  from  the  system  (1)  by  excluding  the  deriva¬ 
tive  of  density  and  by  proceeding  to  the  logarithmic  variables 
and  o 

The  first  two  equations  of  system  (5)  are  the  conditions  of 
compatability  on  wave  surfaces  and  the  derivatives  in  these  equations 
in  the  meridianal  cross  sections  are  written  as  in  a  work  [9],  accord¬ 
ing  to  a  four  point  scheme.  The  last  three  equations  are  equations 
on  the  flow  surfaces  and  their  corresponding  derivatives  are  written 
according  to  a  six  point  scheme  [10]. 

The  obtained  difference  equations  are  solved  in  each  meridional 
plane  by  the  trial  and  error  method.  The  change  in  derivatives 
according  to  the  meridional  angle  ip  is  taken  into  account  in  the 
interative  fashion  from  layer  to  layer.  The  given  scheme  was  used 
in  a  work  [11]  for  calculating  two  dimensional  flow  around  bodies 
by  a  perfect  gas  using  the  establishment  method.  Geometrical  inter¬ 
pretation  of  the  condition  of  the  scheme  stability  is  a  requirement 
that  in  the  meridional .  planes  the  grid  line  5  *  const  proceeding 
from  each  point  of  the  calculated  region  be  located  in  a 
characteristic  triangle,  which  is  usually  satisfied  throughout  the 
entire  flow  field. 

We  shall  now  assume  that  all  flow  parameters  are  known  in  layer 
nodes  .  Calculation  of  the  gas  dynamic  functions  in 

the  nodes,  of  a  layer  n  +  1  is  made  according  to  the  method  described 
above  using  the  known  values  of  */  and  $  from  a  layer  n. 

We  move  elements  of  the  stream  line  before  the  intersection 
with  the  surface  *  *,  backward  according  to  the  found  gas  dynamic 
functions  and  the  nodes  of  a  layer  n  +  1  in  the  plane  of  symmetry. 
Values  of  the  gas  dynamic  functions  and  chemical  components  in  the 
points  of  interesection  are  found  with  the  aid  of  quadratic  interpo¬ 
lation  according  to  known  values  in  nodes  of  these  n  layer  in  the 
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flow  symmetry  plane. 

The  equations  of  the  physico-chemical  processes  (the  last  rela¬ 
tionships  of  system  (1))  shall  be  solved  according  to  the  non-expli¬ 
cit  Newton  method  (  similarly  to  work  [12])  using  the  eauation  of 
state  (2). 

Thus,  in  the  plane  of  symmetry  of  flow,  concentrations  y.^  and 
oscillation  energy  .'r;  and  then  *,  and  <t>  can  be  found  in  all 

nodes  of  the  layer  n  +  1.  The  values  *•  and  $  in  other  meridional 
planes  are  determined  according  to  the  trigonometric  approximation 

«.l"’  ■  m  M  ,0}  t.  «  '•  • 

;  C«,s  (6) 
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The  obtained  values  of  and  4  in  all  nodes  of  the  n  +  1 

layer  are  the  original  ones  for  calculating  the  gas  dynamic 
functions  in  nodes  of  the  n  +  2  layer,  etc. 

As  was  noted  above,  the  effect  of  physico-chemical  processes 
on  gas  dynamic  parameters  in  the  supersonic  flow  region  can  be  quite 
well  estimated  according  to  the  freezing  principle.  In  the  method 
of  this  work,  only  1  assumption  of  the  smoothness  of  change 

in  *.  and  $  according  to  the  meridional  angle,  and  the  kinetic  equa¬ 
tions  are  completely  calculated  in  the  planes  of  symmetry.  For  the 
case  of  flow  around  segmental  bodies,  function  $  in  the  supersonic 
region  decreases  very  rapidly  and  becomes  practically  equal  to  0  and 
naturally  has  a  very  weak  effect  on  the  gas  dynamic  parameters.  The 
relationship  of  the  frozen  heat  capacities  «,  on  the  surface  of  the 
shock  wave  is  equal  to  ,  and  on  the  surface  of  the  body  remains 
approximately  the  same  in  all  points.  The  change  in  •,  according  to 
the  meridional  angle  ip  with  a  fixed  5  is  smooth,  and,  therefore,  can 
be  approximated  according  to  formula  (6>  with  sufficient  accuracy. 
Because  of  the  point  weak  effect  of  the  change  **  on  the  gas  dynamic 
parameters  of  flow,  one  can  hope  that  for  quite  short  bodies  (of  the 


104 


segmental  type)  and  with  comparatively  small  angles  of  attack  the 
developed  method  of  calculation  will  give  good  results. 

We  shall  examine  the  results  of  calculations  of  the  supersonic 
region  of  nonequilibrium  flow  around  a  segmentally  shaped  body,  made 
according  to  the  given  method,  with  angles  of  attack  a  »  15°  (solid 
lines  in  figures  6-25),  and  25°  (dash  lines).  The  calculations  were 
made  with  the  following  parameters  of  the  body  and  the  free-stream 
flow: 

9**  3D*  .  »*■  lO*.  *.  -  l  CM.  ao,  -  250 ’k , ?-  to'*  atm 

The  initial  data  for  calculating  the  supersonic  region  are  assigned 
on  a  surface  «>«.  (figure  6,  the  apex  of  the  conical  surface 
is  situated  in  a  point  with  coordinates  x  *  0.3547  and  y  *  -0.022 
for  a  *  15°  and  x  *  0.3573,  y  *  -0.04165  for  a  *  25°)  from  calcula¬ 
tions  of  nonequilibrium  flow  in  the  shock  layer  near  the  frontal 
surface  of  the  body  (1).  Calculation  of  rates  of  the  physico-chemical 
processes  and  the  thermodynamic  functions  in  the  supersonic  regions 
was  made  according  to  precisely  the  same  program  as  was  used  in  the 
subsonic  region  of  the  shock  layer. 

Figure  6  shows  the  picture  of  flow  in  the  symmetry  plane  cal¬ 
culated  for  an  angle  of  attack  a  *  15°.  Here,  as  in  the  other  figures 
as  well,  the  dash-dot  lines  give  the  results  of  calculating  flow 
around  the  given  body  by  a  perfect  gas  with  an  adiabatic  exponent 
*  *1.4  when  a  *  15°. 

From  the  cited  flow  picture  and  distribution  of  stream  lines 
between  the  surface  of  the  body  and  the  shock  wave  in  the  supersonic 
region,  a  different  character  of  flow  with  the  leeward  (^  »  0)  and 
windward  (^  ■  II  )  sides  of  flow  is  obvious. 

« 

The  shock  wave  is  situated  very  close  to  the  surface  of  the  body 
on  the  windward  side  of  flow  and  the  entire  stream  of  gas  in  the 
shock  layer  near  the  point  of  rounding  passes  the  region  of  the  fan 
of  the  rarefaction  wave ,  running  from  the  contour  rounding  in  the 
region  of  the  middle.  In  this  case,  flow  temperature  sharply 
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decreases  (Figure  7),  which  leads  to  freezing 
of  the  physico-chemical  processes  and  a  decrease 
in  their  energy  contribution  to  practically  0 
(Figure  8).  And  although  the  values  of  T  and 
on  the  initial  surface  are  significantly 

higher  on  the  winward  side  (Figures  7,  8,  and  9, 
10)  after  a  turn  to  a- ?  (calculation  is  made  in 
this  region  in  the  spherical  coordinate  system), 
they  differ  little  from  the  corresponding  values 
on  the  windward  side  of  flow. 


I 


A  significant  drop  in  temperature  along  the  shock  wave  leads 
to  a  situation  such  that  the  physico-chemical  processes  in  the  gas 
which  has  entered  the  shock  layer  in  the  supersonic  region  of  flow 
are  practically  not  excited.  It  is  evident  from  Figure  6  that  the 
stream  line  running  from  the  subsonic  flow  region  on  the  windward 
side  completely  transits  the  fan  of  rarefaction  waves  and  proceeds 
in  parallel  to  the  surface  of  the  body,  proportional  to  the  increase 
in  the  x  coordinate  (calculation  after  turn  to  * •  i  is  made  in  the 
cylindrical  coordinate  system).  The  relative  region  occupied  by 
these  stream  lines  in  the  shock  layer  rapidly  decreases  and  comprises 
approximately  402  for  a  *  15° ,  when  x  ■  1.5#  The  entire  remaining 
shock  region  on  the  windward  side  of  flow  is  occupied  by  a  gas  in 
which  the  physico-chemical  transformations  do  not  practically  arise. 

On  the  leeward  side  of  flow,  the  stream  lines  along  which  chemical 
reactions  arise  in  the  flow  led  to  a  noticeable  change  in  composition 
on  occupy  a  significantly  greater  part  of  the  shock  layer.  Although 
the  gas  composition  remains  constant  along  the  stream  lines,  it  proves 
to  be  variable  with  respect  to  £  in  the  greater  part  of  the  shock 
layer  even  at  the  end  of  the  calculated  region. 
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Figure  11 


figure  1 2 


Figure  1 


figure  14 


Figures  11-16  give  the  corresponding  distribution  of  f».,  r— r> 
as  an  illustration.  Relationships  of  the  frozen  heat  capacities 
«,  transverse  to  the  shock  layer  (Figures  17  and  18)  also  prac¬ 
tically  remain  constant  along  the  stream  lines  everywhere  in  the 
region  x  >  0.3  and  the  character  of  Its  distribution  according  to 
5  is  actually  determined  by  the  flow  picture. 

The  effect  of  increasing  the  angle  of  attack  from  a  value 
a  ■  15°  to  a  *  25°  on  the  distribution  of  flow  parameters  is  evi¬ 
dent  from  the  provided  figures.  Specifically,  the  increase  in  the 
angle  of  attack  leads  to  a  relative  decrease  in  the  region  occupied 
by  the  reacted  gas  on  the  windward  side,  which  is  due  to  a  large 
influx  in  the  given  cross  section  *  *  const  of  the  gas  which 
passed  through  the  wave  region  with  lower  intensity. 

It  was  shown  in  work  with  the  aid  of  an  approximate  method 
that  the  effect  of  physico-chemical  processes  on  the  distribution 
of  pressure  along  the  lateral  surface  is  significant  in  the  case  of 


near-equilibrium  flow  around  the  frontal  surface  of  segmented  bodies. 
In  the  examples  of  calculations  cited  in  this  article,  flow  around 
the  frontal  surface  occurs  under  conditions  far  from  equilibrium 
conditions.  In  this  case,  the  effect  of  the  physico-chemical  pro¬ 
cesses  on  the  distribution  of  gas  dynamic  parameters  in  the  super¬ 
sonic  flow  region  should  be  significantly  smaller. 

Figure  19  shows  the  pressure  distribution  along  the  surface  of 
the  body  in  the  symmetry  plane  of  flow.  The  sharp  increase  in 
pressure  on  the  leeward  side  of  flow  is  due  to  overflow  [1],  The 
distribution  of  pressure  across  the  shock  layer  (Figures  20  and  22) 
also  differs  little  from  the  case  of  a  perfect  gas.  The  distribution 
of  density  (Figures  21  and  23)  on  the  initial  surface  (0  *  1.12 
for  oi  *  15°)  before  cases  of  flow  by  a  nonequilibrium  and  perfect 
gas  is  significantly  different.  In  the  case  of  nonequilibrium  flow 
near  the  surface  of  the  body,  there  is  a  region  of  increased  density 
due  to  the  presense  of  physico-chemical  processes.  However,  after 
the  turn  (0  ■  1.16),  the  character  of  density  distribution  already 
corresponds  to  density  distribution  in  the  perfect  gas,  and  when 
*•  t  ,  the  corresponding  density  values  are  practically  identical 
in  those  cases. 

A  sharp  drop  in  density  in  the  vicinity  of  the  middle  of  the 
body  leads  to  a  characteristic  change  in  the  number  of  electrons 
per  unit  of  volume  ", •  */<*.  (Figures  24  and  25),  where  n^  -  total 
number  of  particles  per  unit  of  volume  of  the  free-stream  flow. 

The  increase  in  n1  near  the  surface  of  the  body  on  the  leeward  side 
at  the  end  of  the  calculated  region  is  due  to  flow  stagnation. 

When  *»«4  for  a  ■  15°  and  x  >  1.3,  for  a  *  25°  the  decrease 
in  flow  velocity  on  the  leeward  side  of  flow  leads  to  a  situation 
such  that 
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CALCULATING  VISCOUS  FLOW  AROUND  A  SPHERE  BY  A  HYPERSONIC 
FLOW  OF  A  GAS  MIXTURE  CONTAINING  CARBON  DIOXIDE  GAS 

by 

V.  G.  Gromov 

The  work  gives  the  method  and  certain  results  of  calculating 
the  viscous  flow  of  a  gas  mixture  of  C02  -  N2  near  the  leading 
critical  line  of  flow  around  a  sphere.  Calculation  was  made  on 
the  basis  of  a  system  of  equations  which  is  the  first  approximation 
of  a  complete  system  Navler  -Stokes  equations  obtained  by  the 
method  "truncated"  series.  A  9-component  model  of  the  gas  mixture 

*.o;  was  used,  taking  into  account  13  gas 
phase  chemical  reactions  and  with  the  assumption  that  the  internal 
degrees  of  freedom  of  the  mixture’s  components  are  equally  excited. 


By  applying  the  method  of  "truncated"  series  [2]  to  the  complete 
system  of  Navier  -Stokes  equations  [1],  and  by  limiting  oneself 
to  the  first  approximation,  one  obtains  the  following  equation  system 
which  describes  flow  of  v-component  of  the  reacting  gas  near  the 
critical  stream  line  of  the  sphere:  The  equation  of  continuity  is 
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the  equations  of  continuity  for  mixture  components 
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the  equations  of  continuity  for  chemical  elements 
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the  energy  equation 


(5). 


the  equation  of  state 
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Equations  ( 2 ) — ( 6 )  can  be  rewritten  in  the  vectorial  form 
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In  these  equations  and  below,  s  and  n  are  coordinates  which  run 
along  the  surface  of  the  body  and  along  the  normal  to  it,  u  and  v 
are  the  corresponding  velocity  vector  components,  p  is  density,  p 
is  pressure,  T  Is  temperature,  /. is  total  enthalpy, 
is  enthalpy  of  i-th  component,  Is  the  gravimetric  concentration 
of  the  i-th  component,  lt  Is  the  gravimetric  concentration  of  the 
1-th  chemical  element,  ^  is  the  rate  of  formation  of  the  i-th 
component  in  the  chemical  reactions,  *«  is  the  number  of  chemical 
elements  that  enter  into  the  composition  of  the  mixture  components. 
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e-norming  factor; 
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*t.  -  gravimetric  proportion  of  the  1-th  element  in  the  i-th  com¬ 
ponent  ; 

;  *  i-i  ..  r  :  < 

Here  ^ -  coefficients  of  viscosity,  X  -  coefficient  of  heat 
conductivity,  Dlk  -generalized  diffusion  coefficients. 

In  the  cited  expressions  and  below,  dimensionless  variables 
are  used  and  all  lengths  are  related  to  the  radius  of  the  sphere  L, 
velocity  components  to  V^,  density  to  p^,  pressure  to 
enthalpy  to  v-  ,  temperature  to  vj <'•/<-..> *;  /  -  absolute  gas 

constant,  the  transfer  coefficients  -  respectively  to 
The  index  "infinity  symbol"  pertains  to  conditions  in  the  unperturbed 
flow.  The  equations  include  dimensionless  complexes 


A. 


characteristic  reaction  time. 


It  Is  assumed  during  statement  of  the  problem  that  the  surface 
of  the  sphere  is  impermeable  ,  has  an  assigned  temperature  /r-“r- / 
and  has  ideal  catalytic  properties.  Moreover,  the  absence  of  slipping 
Is  postulated  which  is  a  good  approximation  in  the  case  of  a  cold  wall. 
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Boundary  conditions  on  the  surface  of  the  body  have  the  follow¬ 
ing  form  with  such  assumptions: 

U  •  J  #  *  •  c 


-  V  .  •  •  CO  '  HO  ,CS  'Ct  H  mO 


i 


(11) 


The  flow  parameters  tend  toward  their  values  in  the  unperturbed 
flow  proportional  to  distance  from  the  body  surface. 
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(12) 


2.  The  system  of  equations  ( i)-(9) ,  (11),  (12)  is  solved  by 
the  method  of  finite  differences.  Equations  (l)-(2),  (10),  (7)-(9) 
are  approximated  with  the  aid  of  a  3-point  difference  scneme  of  the 
following  system  of  difference  equations: 
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(17) 

(18) 

The  system  of  equations  (13)-(l8)  is  closed  by  difference  equations 
which  approximate  the  boundary  conditions 
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Here  the  following  grid  symbols  are  used: 
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s ■*’  -  gravimetric  functions  whose  specific  selection  will  be 

stated  below. 


3.  The  system  of  difference  equations  shall  soived  by  the 
method  of  successsive  approximations  according  to  the  following 
scheme.  Initially,  using  the  (n-1)  approximation  for  the  vector 
of  the  sought  functions  *' from  equations  ( 1 8 )  and  (21) 
we  determine  the  n-th  approximation  for  the  function  Pg.  We  then  find 
the  value  from  equations  (13)-(17),  (19),  (20).  In  order  to 
calculate  *f"’  ,  we  expand  the  left  hand  side  of  equations  (13)-(17), 

(19)  and  (20)  in  a  series  in  powers  of  ax. 71*1- t<m"' 

Limiting  ourselves  to  the  linear  terms,  we  obtain  a  system  of  linear 
equations  relative  to  tT  ,  u 
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We  shall  separately  write  the  equation  which  corresponds  to  the 
difference  approximation  of  the  continuity  equation  in  the  point 
m  +  0.5 

*  v.,f  ■ 


(24) 


(25) 


In  order  to  reduce  the  volume  of  calculations  during  linearization 
of  the  difference  equations,  the  transport  coefficients  can  be  viewed 
as  values  independent  of  ,  whose  values  are  calculated  according 
to  the  (n-l)-th  approximation  . 


We  shall  employ  the  vector,  trial  and  error  method  for  solving 
the  system  of  equations  (22)-(24).  With  the  direct  trial  and  error 
method,  we  exclude  */•«  from  equation  (23)  with  the  aid  of  (25)  and 
calculate  the  trial  run  coefficients 
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Here  .  *».«  are  the  corresponding  matrix  coefficients 

of  equation  (23)  after  excluding  S'/m„, 

Using  boundary  conditions  (24)  and  relationships  (25)  and 
(26)  when  m**-f  we  express  s?+  through  *4 

(27) 

and  carry  out  the  first  inverse  trial  run 
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*  m 't"  C  m.m.i*  m,wkt  ‘  C*M  (28) 

*r- 

f  •  0  f  •  4  A 

From  the  supplementary  condition 

tm»'T'  (29) 

we  determine  the  value  of  J*  ,  and  during  the  second  inverse  trial 
run  we  find  the  values  *•  .**  and 

x"’-  t'Z" 

m •  a ,f , .... «*/  ( 30 ) 

e“> .  t**’ .  it 


4.  The  method  examined  above  was  employed  for  calculating  flow 
of  the  C02-N2  gas  mixture.  The  model  of  a  rigid  rotor  -  a  harmonic 
oscillator  with  the  following  values  of  the  characteristic  vibratory 
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temperature  was  used  during  the  description  of  the  thermodynamic 
mixture  properties: 


CO 

V. 

■U 

*2 

1  CH  ] 

r  •< 

MO,  19X) 

iUX* 

J.TJ 

J2JO 

-J>7« 

t* 

:JMO 

1 

It  was  assumed  that  the  atomic  components  only  have  sequential 
degrees  of  freedom.  The  heats  of  formation  of  components  under  steady- 
state  conditions  were  taken  from  a  work  [3]. 


During  calculation  of  the  transport  characteristics  of  the 
mixture,  the  following  simplified  assumptions  were  used: 


1/  the  viscosity  coefficient  does  not  depend  on  the  composition 
of  the  mixture  and  is  determined  by  the  formula 

2/ 

j* ]IT/Tmfn 

3/  the  binary  diffusion  model  is  used 


c/C- 


$  » 4.  *  0, 7 

<  all 


4/  the  heat  conductivity  coefficient  is  determined  from  the 
condition 


1 0.7 


13  gas  phase  reactions  were  taken  into  account  in  the  calculations: 


l. 

C  (?t  .  M  s*  CO  *0  *  M  , 

7. 

*rcss 

CH*H 

3. 

8 

♦ 

:* 

» 

c  *0  *  n  ; 

5. 

St*0  — 

HO*  H 

3. 

/V  *  /I  ; 

e. 

CO*  H  ss 

CH  *0  , 

4. 

HO  * 

M*0*M  ; 

10. 

CH  *0  9 

C  *  NO 

5. 

0,  ♦ 

0*0  *  M  ; 

11. 

0  *  HO* 

0t  *  H  , 

9. 

C*  *  M  * 

13. 

CO  *  It  * 

HO  *  C 

13. 

C  0f  *  0  sc  CO  *  0 1  . 
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The  rate  constants  of  these  reactions,  given  in  Table  1,  were 
calculated  on  the  basis  of  recommendations  made  in  the  works  of  a 
number  of  authors  who  studied  chemical  processes  in  CC^-^  mixtures. 

Calculations  were  made  on  a  uniform  grid  with  M  *  80.  Values 
of  the  gravimetric  functions  were  determined  as  follows: 

*“’•  o  when  f  » f  * 
a'*’. f  when 

6"1’  i  ,  K  •  I  i  2,  t  i  3,  t  ,4 

Such  a  choice  of  insures  stability  of  the  trial  and  run  method, 
good  stability  properties  of  the  scheme,  and  an  error  of  approximation 
on  the  order  of  h  throughout  the  entire  examined  range  of  values 
of  the  R  number.  Calculations  were  made  whenj»’**. 

oo 

Convergence  of  the  iteration  process  employed  for  solving  the 
difference  equations  depends  significantly  on  the  choice  of  the 
initial  approximation.  In  the  cited  calculations,  the  initial  dis¬ 
tribution  of  the  sought  functions  was  determined  in  the  process 
of  motion  according  to  the  value  of  the  number.  The  calculation 
started  with  finding  a  solution  when  %  -  io*'.  Quite  arbitrary 
profiles  of  these  sought  functions  were  used  in  this  case  as  the 
zero  approximation,  but  when  making  the  first  two  iterations,  for  *r 
damper  <5  *  0.5  was  used. 

X  •  X  •  0,5  •* 

The  motion  step  according  to  was  chosen  experimentally  and 

changed  within  limits  of  from  10^*“^.  to  10. 

Figures  1-13  give  certain  calculated  characteristics  of  flow 
obtained  with  the  following  values  of  exterior  parameters: 

r_  -  iso* k,  v.  -  aoom/sec (  f/taole*.  -  33oo*k) 

(t.+ )  m  -  oa  (  c„f  )_  •  ojx,  Tw  -  iooo’k, 
p  (7^)  -  7to  •  io’’  9/cm  •  sec 

*.  -  JO*1*  -  IO**  .  *  10'*  mm 
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Figures  5-11  show  the  distribution  of  concentrations  of  com¬ 
ponents  of  the  cot  ,co  ,nt,  «o ,  o,  ,  m,o.  mixture.  The  concentration 

-5 

of  atomic  carbon  and  cyanogen  does  not  exceed  0.5*10  throughout 
the  entire  examined  range  of  parameters. 


It  is  evident  from  the  cited  graphs,  under  the  given  conditions 
intensive  processes  began  in  the  gas  already  following  the  establish¬ 
ment  of  the  dynamic  flow  structure  characteristic  for  high  Reynolds 
numbers — the  "shock  wave-nonviscous  gas  layer — boundary  layer”.  In 
this  case,  one  should  note  that  the  equilibrium  gas  composition  does 
not  established  before  *«•»**.  on  the  "outer  boundary"  of  the  boundary 

layer. 


Figure  12  illustrates  the  change  in  "regression 

of  the  shock  wave  e  with  the  change  of  R^.  The 

value  of  e  decreases  monotonlcally  with  the  increase 

4 

in  R^  and  up  to  R^  ^  10  ,  the  decrease  in  e  is  bas¬ 
ically  associated  with  the  change  in  the  dynamic 

structure  of  flow,  but  with  large  R^  values,  the  decrease  e  is  due 
to  an  Increase  in  density  in  the  shock  layer  because  of  gas  dissocia¬ 
tion. 


Figure  13  shows  the  calculation  relationship  between 
the  heat  exchange  parameter  stiff  •»•/*- and  the 
value  of  R^.  The  cited  data  on  heat  exchange  which 
correspond  to  high  ('vlO^)  evidently  need  refining 
inasmuch  as  the  values  were  obtained  with  comparatively 
small  number  of  points  in  the  region  of  the  boundary  layer 
(*5-10). 
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III.  LOSS  OF  STABILITY  AND  SELF  EXCITED 
VIBRATIONS  IN  THE  FLOW 

The  Interaction  of  Plane  and  Cylindrical  Supersonic 
Streams  with  Acoustical  Fields 

by 

M.  G,  Lebedev  and  G.  F.  Telenin 

The  flow  of  gas  streams  into  the  external  medium  is  accompanied 
by  the  emission  of  noise  from  the  stream  surface.  When  supersonic 
stream  flow,  the  spectrum  of  their  noise  contains  discrete  compo¬ 
nents  of  significant  Intensity  under  certain  conditions.  In  this 
case,  the  stream  pulsates  at  a  certain  frequency  which  is  equal  to 
the  frequency  of  the  discrete  component,  and  emits  acoustical  waves 
into  the  ambient  medium  which  are  similar  in  shape  to  spherical  waves 
and  are  in  counter  phase  on  the  opposite  sides  of  the  stream.  The 
emission  of  these  waves  chiefly  occurs  from  the  region  of  stream 
disintegration. 

Analysis  of  the  experimental  data  makes  it  possible  to  conclude 
that  the  stream  emitting  sound  at  a  discrete  frequency  is  a  self- 
vibratory  system  in  which  the  inverse  relationship  is  accomplished 
by  acoustical  waves  propagating  through  the  ambient  medium,  and, 
in  their  turn,  acting  upon  the  stream.  A.  Powell  [10],  who  first 
observed  the  examined  phenomenon  [19531,  supports  this  idea. 
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In  works  [1,2]  the  authors  suggested  a  possible  physical  scheme 
which  explains  the  appearance  and  maintenance  of  pulsations  of  super¬ 
sonic  jets  at  a  discrete  frequency,  based  upon  the  idea  mentioned 
above  and  analysis  of  the  available  experimental  data.  According  to 
this  scheme,  the  supersonic  stream  can  play  the  role  of  a  resonator 
because  of  its  periodic  (cellular)  structure.  Actually,  if  a  periodic 
wave  runs  along  the  surface  of  such  a  stream,  then  the  perturbations 
in  the  stream  caused  by  this  wave. will  also  be  quasiperiodic  in  the 
direction  along  the  stream.  Then,  having  chosen  frequency  of  the 
exterior  acoustical  wave  in  a  certain  fashion,  one  can  so  synchronize 
this  wave  with  the  propagation  of  perturbations  in  ther stream  that 
perturbations  increase,  with  respect  to  absolute  value,  along  the  stream 
to  the  region  0f  ±^s  disintegration,  exceeding  surface  perturba¬ 

tions.  Thus,  the  stream  will  intensify  the  perturbations  coming  in 
from  without  at  a  certain  frequency  (or  frequencies)  w*;  this  critical 
frequency  w*  is  comparable  with  the  natural  resonator  frequency. 

It  is  natural  to  assume  the  presence  of  a  relationship  between 
the  intensity  of  pulsations  in  the  stream  and  the  intensity  of 
acoustical  waves  emitted  from  the  region  of  its  disintegration.  Then 
the  system  which  Includes  this  region,  which  is  the  source  of  the 
acoustical  omission,  and  the  acoustical  waves  themselves  which  propagate 
through  the  ambient  medium  and  the  perturbed  stream  will  be  a  free 
vibratory  one. 

[1,2]  examine  the  problem  of  interaction  of  a  plane  supersonic 
stream  flowing  into  a  flooded  space  with  plane  acoustical  waves  of 
assigned  frequency  propagating  through  a  circular  stream  in  space 
in  the  direction  of  the  nozzle,  corresponding  to  the  suggested  diagram 
of  the  phenomenon.  In  solving  this  problem  (for  the  calculated  stream 
in  the  linear  approximation),  critical  values  of  frequencies  of  the 
exterior  acoustical  waves  were  detected.  At  these  frequencies,  the 
amplitude  of  perturbations  within  the  stream  caused  by  the  periodic 
wave  on  its  surface  increases  infinitely  with  distance  from  the  nozzle, 
while  at . the  same  time  the  relationship  between  perturbations  and 
the  coordinates  x  at  all  other  frequencies  along  the  stream  will  have 


a  periodic  character.  The  intensification  of  perturbations  within 
stream  in  comparison  with  the  assigned  perturbations  on  its  boundar¬ 
ies  at  the  critical  frequency  is  maximum ;  in  this  case,  a  sharp 
increase  in  perturbations  in  the  stream  in  the  vicinity  of  this 
frequency  takes  place  in  a  narrow  range  of  frequencies,  outside 
which  the  perturbations  have  precisely  the  same  order  as  in  the 
exterior  acoustical  wave.  The  obtained  results  thus  confirm  the 
suggested  physical  scheme  of  the  phenomenon  of  pulsations  of  super¬ 
sonic  streams  with  a  discrete  frequency. 

We  note  that  the  problem  of  interaction  of  the  stream  with  the 
external  acoustical  waves  was  solved  in  [1,2]  in  a  simplified  pre¬ 
sentation,  when  the  removal  of  acoustical  energy  from  the  surface 
of  the  pulsating  stream  into  the  surrounding  space  was  ignored. 

This  made  it  possible  to  obtain  the  solution  to  the  problem  in  the 
analytical  form,  and  specifically,  to  derive  formulas  for  the  critical 
frequencies  of  external  acoustical  waves.  Precisely  the  same  problem 
is  solved  in  this  work  without  the  simplifying  approximations.  In 
this  case,  we  shall  examine  cases  'for  both  the  plane  and  the  cylindri¬ 
cal  supersonic  jet  with  the  symmetrical  and  antisymmetrical  character 
of  perturbations  in  the  exterior  acoustical  field. 

1.  Statement  of  the  problem  and  method  of  solution.  Let  a 
supersonic  cylindrical  or  plane  Jet  with  a  Mach  number  M  flow  from 
a  rigid  cylinder  2  ft  in  diameter  or  a  flat  slot  with  •  the  width  2R 
into  a  flooded  space C£j/)»  pressure  in  the  Jet  and  the  ambient  medium 
is  p.  Flow  in  the  underperturbed  jet  is  uniform;  the  problem  of 
interaction  of  such  a  jet  with  the  exterior  acoustical  field  will  be 
solved  in  the  linear  approximation,  whose  bases  are  given  in  [1,2]. 


Figure  1 
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In  accordance  with  the  physical  phenomenon  of  emission  of  sound 
at  a  discrete  frequency  by  a  jet  simulated  by  us,  we  shall  consider 
that  the  jet  disintegrates  at  a  distance  xq  from  the  nozzle's  cross- 
section  and  that  the  region  of  disintegration  of  the  jet  is  a  source 
of  emission  of  acoustical  waves  with  a  cyclic  frequency  u.  We  shall 
term  this  emission  basic  emission;  let  it  be  described  by  the  poten¬ 
tial  : 

<■***)] .  (1) 

where 


In  (1),  x,  y  are  cylindrical  coordinates  (Figure  1),  0  -  meri¬ 
dional  angle,  j  *  0  or  1  in  the  case  of  symmetrical  or  antisymmetrical 
perturbations.  We  shall  give  the  entire  presentation  for  the  cylindri¬ 
cal  case  (m  *  1),  and  slight  changes  which  should  be  inserted  into  the 
initial  equations  and  calculation  formulas  for  the  plane  case  (m  *  0) 
will  be  indicated  at  the  end  of  the  section. 


Propagating  through  the  space  which  surrounds  the  jet,  the 
acoustical  waves  (1)  act  on  the  Jet  and  bring  it  into  the  perturbed 
state.  Let  perturbed  flow  in  the  jet  be  described  by  a  potential  * 

( 2 ) 

+  (t.y.9.  t)-  r.yWitM-  9 

Oscillations  of  the  perturbed  jet  lead  to  the  emission  of  acous 
tical  waves  from  the  surface  of  the  jet  into  the  exterior  medium. 

We  shall  term  this  emission  secondary  emission;  let  it  be  described 
by  the  potential  ip: 

(■3) 

•.£)•*»  [Y(  t,  If)  -  »tp(-  icjt  )]  cot  f9 


Potential  +  of  the  perturbed  flow  in  the  jet  and  potential 
tp  of  secondary  emission  satisfy  the  following  wave  equations: 


i»'+  )‘+ 

7P  'iy* 


i  a*  t  a’*  i‘t  a 
1  iy  jr‘  it1  it  it  IP 


i’f  .  ,  <  i*  J_  i'f 

TP  *Tp  7  iy  '  7  ~7ei  77*  *  0 


(4) 

(5) 


In  these  equations  and  everywhere  subsequently,  coordinates  x,  y 
relate  to  the  radius  of  the  Jet  R;  moreover,  dimensionless  time  values 
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are  inserted  r,aif,t  . h'at*"t.  where  a  and  a  are  the  velocities  of  sound 

e 

in  the  jet  and  in  the  exterior  medium.  We  shall  also  use  dimension¬ 
less  values  of  the  cyclic  frequency  and  *, such  that 

and  «<■«,»,  .  Finally,  in  equation  (4)  a  **'*■  > ' 

Velocity  perturbations  v‘  and  7;  and  pressure  perturbations  p* 
and  p^  in  the  Jet  and  in  the  exterior  medium  (from  the  secondary 
emission)  are  expressed  through  the  potentials  *  and  tp  in  the  follow- 
way  : 


v 


i  ft  ad 


vr°' 


’  _  K  *+  ,-a-f  a*  i 

P'-TPTTt'-TPfii  •*rJ 


jtfd  f 

Jt 


l  p'*'t'p!e7i  "h?Ht 


(6) 

(7) 


where  p  is  unperturbed  pressure,  identical  in  the  jet  and  in  the 
exterior  medium,  rr>  -  adiabatic  exponents. 

The  potential  of  basic  emission  ipQ  and  the  potential  of  total 
emission  *x  *  *.  •  *  together  with  satisfy  the  wave  equation  (5). 

The  corresponding  velocity  and  pressure  perturbations  can  be  deter¬ 
mined  according  to  formulas  similar  to  (6)  and  (7). 

The  basic  purpose  of  this  work  is  to  study  the  amplitude-frequency 
characteristic  of  the  supersonic  jet  interacting  with  the  exterior 
acoustical  field.  During  the  presentation  of  this  problem  we  shall 
pose,  as  it  were,  the  inverse  relationship  in  the  "jet-exterior  acous¬ 
tical  waves"  self-oscillatory  system  and  shall  hypothesize  the  pre¬ 
sence  of  a  generator  of  basic  acoustic  emission  ^  independent  of_ 
of  the  Jet.  Thus,  we  shall  investigate  a  perturbed  flow  in  the  jet 
and  secondary  emission  from  the  Jet  under  the  effect  of  assigned 
acoustical  waves  which  have  a  potential  ij;  . 

The  basic  acoustical  field  ipo,  whose  source  lies  in  the  region 
of  Jet  disintegration,  has  not  been  sufficiently  studied  experimentally 
at  present.  Therefore,  we  shall  make  assumptions  relative  to  its 
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character  which  simplify  the  calculations  and  at  the  same  time  do 
not  contradict  known  experimental  facts.  We  shall  consider  that 
derivative  **./**«».  on  the  surface  of  the  jet  (y  *  1),  i.e., 
the  effect  of  the  basic  emission  on  the  jet  occurs  only  via  pressure, 
while  the  pressure  perturbation  itself  is  a  pressure  wave  running  in 
the  direction  of  the  nozzle  when  y  =  1; 


p‘(*.  9,t)  «  Pem  (it*t 


(8) 


One  can  assume  that  a  number  of  solutions  to  the  wave  equation 
exists  which  precisely  or  approximately  satisfy  the  indicated  condi¬ 
tion  and  the  condition  of  limitation  at  infinity.  In  the  symmetrical 
case,  for  example,  a  plane  wave  propagating  from  x  8  “  to  a  nozzle 
is  such  a  solution: 

In  the  antisyrnmetrical  case,  the  following  can  serve  as  an  example 
of  such  a  solution: 

•HT*7t)c-s .  ( 10  ) 

Here,  k  >  0  -  a  small  constant  value,  and  -  Bessel  functions 
and  Neumann  functions  of  the  first  order.  We  note  that  the  distri¬ 
bution  *.(*)  qualitatively  corresponds  to  distribution  from  a  dipcle 
whose  axis  is  normal  to  the  axis  of  the  Jet;  when 

When  k  -*■  0,  the  wave  equation  corresponding  to  the  potential  (10)  can 
be  made  as  close  to  (8)  as  one  wants.  (However,  when  ‘-°  — 

when  *-•  ).  Subsequently,  we  shall  need  an  expression  for  pressure 
corresponding  to  the  basic  emission  only  on  the  surface  of  the  jet, 
and  we  shall  use  formula  (8)  in  both  the  symmetrical  and  the  anti- 
symmetrical  cases;  moreover,  it  will  everywhere  be  assumed  that 

The  boundary  conditions  which  uniquely  determine  the  perturbed 
flow  in  the  Jet  and  the  external  acoustical  field  should  be  specified 
applicable  to  the  physical  character  of  the  examined  problem  [1,2], 

When  the  jet  is  producing  sound  at  a  discrete  frequency,  effective 
interaction  between  the  jet  and  the  exterior  medium  only  exists  in 
the  region  from  the  cross-section  of  the  nozzle  to  the  region  of  jet 
disintegration  to*  i<  tj  ,  Therefore,  the  conditions  of  equality  of 
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pressures  and  displacements  on  the  boundary  between  the  jeo  and  the 
exterior  acoustical  field  should  only  be  imposed  when  »****,,  and 
*  **.  •  0,16  can  select  the  condition  of  rapid  damping  of  vibrations 
of  the  boundary  of  the  jet  and  the  exterior  acoustical  field  with  the 
growth  of  X  ( c[l,2 ]).  When  **o  on  the  surface  of  a  rigid  cylinder, 
the  ordinary  condition  of  non  -flow  is  imposed. 


Let 


(ID 


be  perturbation  of  the  boundary  surface,  calculated  according  to  the 
parameters  of  perturbed  flow  in  the  jet,  and 

y#  {*.».*)•  (12) 

be  precisely  the  same  value,  but  one  calculated  according  to  para¬ 
meters  of  the  perturbed  state  at  the  exterior  medium.  The  boundary 
conditions  when  y  *  1,  which  impose  limitations  on  mixing  of  the  jet 
and  the  exteror  medium,  can  be  written  in  the  following  form: 


j  yfx.a.t 1 
i  sfx.e,*' 


when 

when 

when  *.«<«• 


where  e  is  a  function  which  rapidly  vanishes 
in  x. 


(13) 


with  the  increase 


Values  Yg  and  Y  are  determined  according  to  the  potentials  •* 
and  with  the  aid  of  differential  equations: 


(li4) 

ST?  *  <9y  /  y.  r  » 

*  * '*?'»•<  (15) 


The  condition  of  equality  of  pressure  perturbations  on  the 
boundary  surface  when  0  *  *  ‘  *.  can  be  written  in  the  form: 

P  '(*.  8,  t)  •  Pf  ('» ,  9 ,  t .)  ♦  Pen  (if  t )  co*  j  B 


(16) 
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In  order  that  the  solution  be  completely  determined,  it  is 
necessary  to  add  the  condition  of  emission  for  the  exterior  acous¬ 
tical  field  at  infinity  to  the  described  condition,  and  for  the 
perturbations  in  the  jet,  a  condition  on  the  axis,  *  ■  ^  when  J  =  1 
or  when  j  =  0 ,  and  a  condition  on  the  nozzle  cross-section 

'  .  jia'd  ^  =  0  when  «*0,0iy«'l,0i»(2J 


Thus,  the  problem  consists  in  solving  the  wave  equations  (4) 
and  (5)  with  boundary  conditions  (13)-(l6),  as  well  as  on  the  nozzle's 
cross-section,  on  the  axis,  and  at  infinity.  We  shall  make  one  remark 
here.  We  shall  note  equation  (5)  which  describes  tne  perturbed  state 
of  the  exterior  medium,  with  equation 

3  V  3* <4/  1  ^  V  ^  „ 

’  ^y1  *  r  lx)  '  *  7te  "  (17) 

in  which  the  term  <t>o),  }  which  describes  the  dispersion  of 

energy  is  present.  Such  a  substitution  provides'  certain  advantages 

during  the  solution  of  the  wave  equation  by  means  of  the  Fourier 

transformation.  It  is  ordinary,  for  example,  when  using  the  Wiener- 

Hopf  method  [4],  where  we  pass  to  the  limit  e  =  0  in  the 

final  solution.  In  this  work,  during  the  numerical  solution  of  our 

problem,  we  shall  assume  that  e  is  small,  but  differs 

from  0, considering  that  the  solution  to  equation  (17)  with  e  is 

similar  to  the  solution  to  equation  (5).  The  investigation  of  the 

relationship  between  the  solution  and  the  parameter  e  confirms  this 

conclusion. 


Having  substituted  expressions  (2),  (3),  (11)  and  (12)  in  equa 
tion  (4),  (17),  (14),  and  (15)  and  having  used  the  relation  (7) 
between  perturbations  of  pressure  and  the  potentials,  we  obtain  the 
equations  and  boundary  conditions  for  the  complex  potentials  0  and  ¥ : 
when 


.»**♦  t  i<p  »  t. 

*  *'0 

0  y\  y  *  '  i  *  «*  *  x  *  m 


(IS) 


(19) 
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wher* 


when  y » - , o *  x *  «* 


.i  -  *  *.  .f  i 
i  *  f 


>:•;*(-** -{f.t ’*')*)]  *  -</r 


when 


when 


where 


*  -  o  when  /• 1  , 

a +/»i-o  when  jm°  , 

«  -  0  ,  ff  «  y*  / 

♦  •  ytaJ  +  •  O 

If  ml  %  O  %  t  l  m 

j  (it*  -  *  a+/n).  n(t)  •  r  ttp  (-  n,  i)  , 


when  y  = 


when 


_ l  *ro-«  .  *. *< 

/•‘Ty7  •«  . 

/*/«•  - 


We  shall  determine  the  solution  to  equation  ( 1 8 )  with  boundary 
conditon  (20),  (22)  and  equation  (19)  with  boundary  condition  (23), 

(24),  considering  that  functions  "(•*  ,9(*) .  that  enter  into  (22),  (23), 
are  known.  We  seek  the  solution  to'  equation  (18)  in  the  form 


j'a  If  tj. i*)f] i(* •  t, )*]' *(•)•* (■'*•) 


We  apply  the  Laplace  transformation  to  equation  (26)  and  (27) 

G(?,  t)  ■  /  f  (x,t)e  dt  (2  8) 

where  p  is  a  complex  value,  p>c  .  Taking  (21)  into  account 
for  determining  G,  we  have  the  equation 

IP  *7  S -(S'-tp)6-0  >*)  (29) 

with  a  boundary  condition  with  y  *  1 


where 


fi  (p)  ■  /  n(t)*  '*‘  e" 


The  solution  to  equation  (29)  which  satisfies  the  condition  F  *  0 
on  the  axis  for  j  =  1  or  ar/a^o  for  j  =*  0, 

G(p.i)'  (m)  .  (32) 

where  Ij  is  the  Bessel  function  of  an  imaginary  argument  of  the  order 
j.  The  constant  A(p)  is  determined  from  condition  (30).  Having 
completely  determined  ,  we  find  F(x,y)  with  the  aid  of  the 

inverse  Laplace  transformation 


f(  1, s') •  Hi  f  6 (p. *)*"  dP  (t’°)  (33) 

»-*• 

In  order  to  find  the  solution  to  equation  (19),  we  apply  a 
Fourier  transformation  to  it 

Y(‘.r)'  **  fff. /««.*,).  (34) 

•m 

Function  satisfies  the  equation 


where 


iy*  y  rfy  y*  ' 

f*  9  m*  -  ir#  •  *•  ^ 


The  solution  to  equation  (35),  which  satisfies  the  condition  of 
limitation  at  infinity  is 

H(~.  i)*  *(+)  (1  v)  .  (37) 

where  Kj  -  a  MacDonald  function  on  the  order  of  j.  In  this  case, 
the  following  branch  of  the  function  was  chosen: 


where 


i,  •[</*("  o' 

t,  •  *[ '/*(■  s  >(i', 

s,»,'  *  »,* 

r  c  i  /  t 

*  **,»,  •(*•«)■(!-«)-  tiji  .  (38) 

Having  applied  the  Fourier  transformation  to  (23),  we  obtained  the 
boundary  condition  for  equation  (35)  when  y  *  1 

*  6  •  (39) 


where 


W  i.  -  x  •  „  .  .««  7  v 

5  ~  • .  -»<*'•  rf. •  f 3 («v tf.  . 

.ft'',  >« 

•|  -*«/rw  .  (4o) 

Function  I(<C)  and  (40)  is  associated  with  the  specific  damping  law 
of  vibrations  of  the  boundary  surface  when  *"• 

Using  (40),  we  find  the  final  expression  for  the  solution  *(».■/> : 

■  (41) 

With  the  aid  of  the  inverse  Fourier  transformation,  we  obtain 


J  "0*.  y>  ‘“‘  <«  f*  *  *  ‘  (  42 ) 

Using  the  results  obtained  above,  one  can  construct  an  iteration 
process  of  solving  equations  (18),  (19)  with  boundary  conditions  (20) 
-(24).  It  consists  in  the  following: 

1.  When  « *  o  the  zero  approximation  n*(t)  0f  function  -ait)  is 

assigned. 

2.  One  determines  the  solution  which  describes  the  perturbed 

flow  in  the  jet  according  to  (25)— ( 33) • 

3.  Integration  of  the  equation 

.  (43) 

which,  associates  perturbation  of  the  boundary  of  a  jet  Y  with  the 
potential  0  is  determined  by  the  right  hand  part  in  the  boundary  con¬ 
dition  (23). 
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4.  The  solution  Yi«.w  which  describes  the  perturbed  state  of 
the  exterior  medium  is  determined  according  to  (34),  (42).  The 
first  approximation  for  n(x)  is  determined  according  to  the  obtained 
solution : 

5.  The  described  process  is  multiply  repeated  so  long  as  the 

sequential  approximations  />*(•)  and  n“' (•)  are  not  sufficiently  close 
to  in  the  entire  range.  We  note  that  it  is  desirable  to  refine 

the  successive  approximations  of  nn  according  to  forumulas  of  the 
following  type  in  order  to  insure  and  speed  up  convergence  of  the 
iteration  cycle: 

In  this  work,  the  method  of  the  false  position  [5]  was  used  for  such 
refinement.  This  method  is  a  generalization  of  the  method  of  secants 
to  the  case  of  the  operator  equations. 

The  described  algorithm  was  performed  on  the  computer.  In  this 
case,  the  Laplace  and  Fourier  transformations  were  made  numerically 
according  to  the  Simpson  formula  on  segments  of  a  finite  line  divided 
into  95  intervals.  With  suffificiently  long  lengths  L,  of  these  seg¬ 
ments,  the  solutions  practically  did  not  depend  upon  L. 


In  order  to  test  the  accuracy  of  the  obtained  solutions,  the 
latter  were  placed  in  the  equations  and  boundary  conditions  (18)- 
(23).  Calculations  showed  that  these  relationships  are  satisfied 
with  good  accuracy. 


In  the  plane  case,  the  original  equations  (4)  and  equation  (17) 
should  be  replaced  by  equations: 


*  IP  '  iy‘ 


in- 


it  it 


*  it1' 

»*  n 
—  «  0 


(45) 

(46) 


Equations  (18),  (19),  (26),  (29)  and  (35)  also  correspondingly  change; 
the  relationship  with  angle  0  in  all  formulas  is  naturally  absent. 
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The  solutions  to  the  differential  equations  for  the  Images 
G  and  H  of  potentials  F  and  ¥  in  the  plane  case  are  written  via 
the  Poiner  and  not  the  Bessel  functions.  Namely,  instead  of 

(32)  we  have 

*(*><* when  J  -  0  (47) 

a (p)tijty  v*ien  ;  ■  1 


and  instead  of  (37)  we  have 

•  a(«)#  *r  viien  >  *  0  «  1 


(48) 


We  now  note  that,  if  one  sets  in  the  iteration  cycle 

described  above,  and  one  limits  oneself  to  the  zero  approximation, 
then  we  obtain  a  solution  to  the  problem  in  the  simplified  statement. 
In  such  a  solution,  one  ignores  the  secondary  acoustical  emission 
from  the  surface  of  the  jet. 


The  inverse  Laplace  transformation  (33)  when  n(*)%o 

can  be  made  analytically,  using  the  Jordan  lemma  and  Cauchy  theorem 
of  residues  ,  similar  to  the  way  this  was  done  for  the  plane  case 

in  Oj  2-J*  We  here  write  the  single  solution  for  the  plane  (m  *  C)  and 
cylindrical  (m  ■  1)  solution  for  the  antisymmetrical  (J  *  1)  perturba¬ 
tions  in  the  exterior  acoustical  waves.*/ 


Ff*.y)  •  *  if,  *  Q(t, ttp (..*,) 

F'  ’  A  /rV/4 >(  v)  u”  7  -  7/  - 

.f  _ ±1  A’"' _  *m,  *  1 

-■  >-  } 

c  .  fl  J  1  f /')  .  d,  /•)  ct  ] 

1  ( ihn 7  -/„  - 

~  L  — 1 


(49) 


*/  -  The  solution  to  the  problem  in  the  simplified  presentation 
for  the  cylindrical  case  was  obtained  by  a  student  at  Moscow 
State  University,  V.  V.  Labutina. 
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where 


• 

*  i  -  i  «  *  *  *  ■ ,  • 

-  - 

*  —  * 1  m  * »  • 

j  liC*.’ 

-•  '  -.*6 

--  .  '  t*' 

*  *  •  *11  c  * 

.  ■-»  Cv ! 

y,(c)  ’ 

y  - 1 

•  * 

un  -  the  n-th  root  of  the  Bessel  function  T^. 

We  note  that  the  solutions  obtained  from  (49)  in  the  cylindrical 
case  for  pressure  perturbation  p'  and  the  velocity  components  $ 

X 

have  infinite  discontinuities  on  characteristics  BC,  CD,  FG,  GK 
(Figure  1).  A  similar  result  was  earlier  obtained  in  a  work  [6] 
during  investigations  of  perturbed,  supersonic  steady-state  flow  in 
cylindrical  tubes  and  Jets  during  an  investigation.  One  can  assume 
that  this  effect,  associated  with  the  description  of  the  reflection 
of  discontinuities  from  the  axis  by  linearized  equations  of  gas  dynamics, 
has  an  effect  on  the  solution  only  in  the  vicinity  of  the  discontinuity 
lines.  This  is  also  confirmed  by  calculations  that  were  made  espec¬ 
ially,  in  which  the  condition  t  *  0  or  >+a¥.a  was  set  not  on  the  axis, 
but  on  a  certain  cylindrical  surface  of  small  radius  6  surrounding 
the  axis. 

It  is  also  important  to  note  that  the  solution  for  the  boundary 
of  the  perturbed  Jet  and  the  distribution  of  pressure  on  the  Jet 
boundary  do  not  have  discontinuities. 


2.  Investigation  of  obtained  solutions.  We  shall  initially 
refer  to  the  solution  to  the  problem  in  the  simple  statement  (49), 

By  investigating  the  singular  points  of  the  solution  (49),  one  can 
as  was  done  in  [1,2],  derive  a  formula  for  critical  frequencies  of 
the  exterior  acoustical  waves,  with  which,  as  will  be  shown  below, 
one  deals  with  maximum  intensification  of  perturbations  in  the  Jet. 

We  shall  write  the  general  formula  for  the  plane  (m  *  0)  and  cylindrical 
(m  *  1)  cases,  both  for  the  symmetrical  (J  ■  0),  and  for  the 
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antisymmetrical  (J  =  1)  perturbations  in  the  exterior  wave: 


Here,  Strouhal  number,  calculated  according  to  the  linear  fre¬ 

quency  ,  velocities  of  sound  in  the  exterior  medium  a^,  and  diameter 
width  of  jet  2R;  n  *  1,  2,  3,...,  and  functions  h(n)  is 


defined  as: 


* 4J.  /  *  • Ji 


/  *  ■;) 


By  ,  we  have  in  mind,  respectively,  n-th  roots  of  the 

Bessel  functions  T  and  T. .  We  shall  give  the  first  three  values 

o  1 

of  the  functions  h(n)'  for  the  cylindrical  case: 

when  j  “  °  *-7o,  j.-r> 

when  j  *  1  h(‘-  ”  ‘-a.  3.J3 

We  note  that  the  relationship  of  critical  frequencies  for  the  plane 
and  cylindrical  cases  which  correspond  to  precisely  the  same  value 
of  n  does  not  depend  upon  the  parameters  which  determine  flow  in  the 
jet  and  in  the  exterior  medium,  and  is  always  greater  than  unity. 


In  Figure  2,  according  to  the  formula  (50),  the  relationship 
is  plotted  when  n  *  1  for  a  cylindrical  cold  jet  (stagnation 
temperature  TQ  which  is  equal  to  temperature  of  the  exterior  medium 
T  )in  the  case  of  an  antisymmetrical  acoustical  field.  Here,  too, 
a  semi-empirical  relationship  between  the  frequency  of  the  discrete 
component  (the  first  tone)  and  parameters  of  the  circular  jet  and 
the  exterior  medium 


1Ml-  1  *  Og/o 


are  obtained  in  [7]  when  processing  a  large  number  of  experimental 


data  obtained  from  different  authors  is  given.  The  dash  line  in 
Figure  2  plots  the  scatter  boundaries  of  the  experimental  data.  The 
coincidence  of  our  results  with  the  experimental  data  can  be  considered 


satisfactory,  particularly  at  low  M  nubmers  of  the  jet.  For  comparison, 
Figure  2  also  plots  the  relationship  s\ r for  the  plane,  anti¬ 
symmetrical  case. 


1 — critical  frequency  of  the  antisymmetrical , 
exterior  acoustical  wave  for  the  cylindrical 
case;  2 — the  same,  for  the  plane  case;  3 — 
empirical  relationship  of  the  frequency  of  the 
discrete  component  in  the  noise  spectrum  of  the 
circular  jet  [7];  —  -  scatter  boundaries  of 
experimental  data 

We  shall  now  proceed  to  examine  solutions  to  the  non-simplified 
problem  */.  We  shall  chiefly  investigate  the  behavior  of  the  value 
K  ,  which  is  the  ratio  of  the  amplitude  of  pressure  perturbations 
in  any  particular  point  of  the  jet  or  exterior  medium  to  the  ampli¬ 
tude  P  of  the  exterior  acoustical  wave  (8).  For  points  within  the 
jet,  Kp  is  the  coefficient  of  intensification  of  pressure  perturba¬ 
tions  in  the  jet  in  comparison  with  the  assigned  perturbations  (8) 

which  characterize  basic  emission.  For  the  exterior  medium,  K  is 

P 

the  amplitude  of  pressure  pulsations  caused  both  by  the  basic  emission 
(8)  and  by  secondary  emissions  from  the  edge  of  the  oscillating  jet. 

We  shall  also  examine  perturbations  of  pressure  in  the  exterior 
medium  due  solely  to  secondary  emissions;  their  amplitude,  related 
to  P,  shall  be  designated  as  K* .  Subsequently,  inasmuch  as  the  char- 
acteristic  dimension  of  our  problem  is  the  length  of  a  unit  of  the  jet 
L.ji*  ,  rr  (Figure  1),  we  shall  use  the  coordinate  »**/•*  •» 
during  the  analysis. 

Figure  3 


*  "  0,315 

“  8 
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- -  1 2 

•  *  I 


*  / 

Everywhere,  unless  otherwise  stated,  the  case  of  the  anti¬ 
symmetrical  acoustical  field  (j  *  1)  is  investigated. 
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It  is  primarily  necessary  to  explain  what  effect  the  boundary 

condition  ( 1 8 )  of  rapid  damping  of  boundary  oscillations  of  the  jet 

behind  the  cross  section  z  =  z^  set  by  the  authors  has  on  the  solu- 

o 

tion  in  the  region  z  <  zQ,  inasmuch  as  this  condition,  of  course, 
quite  crudely  simulates  the  termination  of  effective  interaction  of 
the  jet  and  the  exterior  medium  existing  in  reality.  Figure  3  plots 
relationships  with  the  coordinate  z  of  the  value  Kp  on  the  boundary 
of  the  cylindrical  jet  for  calculation  with  m-3,  s* .  s*. .  ojis. 

Comparing  the  solutions  obtained  when  zq  *  8,  10  and  12,  one  can  see 
that  a  common  curve  exists  with  which  the  curve  *,(•.*,)•  merge  which 
corresponds  to  various  z  .  The  divergence  of  curves  *>/*,*.;  from  the 
common  curve  begins  approximately  when  z  *  zq-2  ^or  each  zQ.  The 
curves  **(•.•.)  behave  similarly  when  y  <  1,  i.e.,  within  the  jet. 

In  y  *  1,  the  relationship  K  (z)  in  the  vicinity  of  the  point  z 
*  zq  has  the  character  of  a  rapidly  dampning  flare£^  The  value 
of  this  flare  depends  upon  the  selected  law  of  dampning  of  vibrations 
of  the  boundary  surface  when  z  >  zQ.  Cases  were  examined  in  the 
work  when  the  displacement  value  of  Y  when  z  >  zQ  either  identically 
equals  0,  or  decreases  to  0  with  the  Increase  in  z  exponentially  or 
according  to  a  parabolic  law.  However,  the  solution  when  z  <  zo 
outside  the  small  vicinity  of  the  point  zQ  is  practically  not  affected 
by  the  selected  damping  law. 

Thus,  one  can  conclude  that  if  at  least  for  z„  $  12,  the  effect 

o 

of  the  boundary  condition  when  z  ^  zq  Is  only  felt  on  the  solution 
to  our  problem  when  z  <  zQ  in  the  vicinity  of  a  point  z  =  z  ,  with 
a  length  on  the  order  of  1  jet  unit  (for  pressure).  In  this  case, 
the  case  examined  in  Figure  3,  when  amplitude  of  pressure  perturba¬ 
tions  K  on  the  boundary  of  the  jet  rapidly  increases  with  distance 
from  the  nozzle  is  the  most  unfavorable  one.  In  those  cases 

when  a  noticeable  Increase  in  K  on  the  boundary  of  the  jet  with 

P 

the  increase  z  does  not  occur,  the  effect  of  the  boundary  condition 
when  z  >  z  on  the  solution  is  still  weaker. 
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We  finally  note  that  the  solution  which  describes  vibrations 

of  the  jet  boundaries  when  z  <  zQ  does  not  practically  depend  on 

the  choice  of  z  ,  and  the  condition  set  when  z  >  z  . 

o  o 


The  figure  4  compares  the  simplified  solution  (49)  and  the 
complete  solution  obtained  numerically  for  the  plane  case.  The 
distribution  of  the  amplitude  of  pressure  perturbations  Kp  within 
the  jet  (y  *  0.5)  along  the  z  coordinate  for  a  version  "  *  3,  «/«,  ■  i, 
sh  $h.  -  0.26.  is  plotted. 


Figure  4 
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-  complete  solution 

- -  simplified  solution 


Comparison  of  the  data  in  Figure  4  shows  that  even  in  a  case  when 
one  takes  the  secondary  acoustical  emission  from  the  boundary  of  the 
jet  into  account,  the  value  of  Kp  at  the  critical  frequency  increases, 
on  the  average,  with  distance  from  the  nozzle  cross-section;  in  this 
case,  during  the  transition  through  each  discontinuity  surface,  value 
Kp  increases.  Doubtlessly,  pressure  perturbations  are  smaller  here 
(1.5  -  2  times)  than  in  the  simplified  solution  which  does  not  take 

into  account  the  removal  of  acoustical  energy  from  the  jet’s  surface. 

/ 

We  note  that  the  segments  of  the  jet  on  which  Kp  diminishes  pertain 
to  characteristic  triangles  of  the  type  C  D  E  in  Figure  1,  one  of 
whose  sides  is  the  boundary  of  the  jet  from  which  the  acoustical  energy 
is  removed.  On  the  other  hand,  the  segments  on  which  Kp  conserves 
a  constant  value  or  increases  pertain  to  triangles  of  the  D  E  F  type, 
not  associated  with  the  jet  boundary. 


A  comparison  of  the  complete  and  simplified  solutions  is  also 
made  in  Figure  4  for  the  frequency  &h.o.ss\  .  Both  solutions  are 
similar  and  are  periodic  solutions. 
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Figure  5 
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simplified  solution 


A  similar  comparison  was  made  in  Figure  5  For  the  cylindrical 

case.  Here  the  value  K  once  again  increases  with  distance  from 

P 

the  nozzle  at  the  critical  frequency  ,  on  the  average, 

and  at  the  frequency  relationship  is  periodic. 

One  can  note  that  the  effect  of  secondary  emission  on  the  perturbed 
flow  in  the  jet  is  weaker  in  the  cylindrical  case  than  in  the  plane 
case. 


m  -  a,  •/•(  -  i,  *  -  e , 

i  •  ofi 
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**  •  I,  -  m  •  1, 

«-0.  - m  -  O. 

Figures  6  and  7  plot  the  relationships  of  the  Kq  value  within 
the  jet  (y  *  0.5)  for  cross  section  z  -  6  and  7  with  the  ratio  a 
of  the  frequency  of  exterior  acoustical  wave  Sh  to  the  critical 
frequency  Sh#  calculated  according  to  (50),  when  n  »  1.  Taking 
into  account  the  removal  of  acoustical  energy  from  the  boundary  of 
the  jet  does  not  lead  to  a  noticeable  change  in  the  value  of  the 
critical  frequency,  whose  formula  was  derived  when  solving  the  problem 
in  the  simplified  statement,  for  the  cylindrical  jet.  For  cold  jets, 
the  maximum  Kp  is  reached  as  it  were  when  a  *  1,  l.e., 

When  the  jet  is  heated,  one  observes  a  slight  decrease  in  the  fre¬ 
quency  which  corresponds  to  the  peak  of  the  value  K  ,  as  in  Figures 

P 

6  and  7;  however,  even  for  an  infinitely  heated  jet  (■•/«,—>  ,  dis¬ 

placement  according  to  the  peak  frequency  of  the  value  K  relative  to 
a  *  1  does  not  exceed  5-7$. 

Curves  *’(m)  have  an  absolute  maximum  when  *  •  for  the 

plane  Jets  when  z  -  6  and  other  whole  even  values  of  z.  Additionally, 
one  has  a  local  maximum  when  «•«■>.«.  whose  value  is  comparable  to 
the  value  of  the  absolute  maximum.  It  is  interesting  that  when  -.0.3 
the  curve  for  the  whole,  odd  z  (for  z  «  7  in  Figure  3)  has  an 


V*;- 


143 


absolute  maximum,  but  when  a  ■  1  the  maximum  does  not  exist  at  all, 
and  instead  of  it,  one  has  a  point  of  inflection.  Thus,  in  the  plane 
case  secondary  acoustical  emission  from  the  jet  surface  leads  to  the 
appearance  of  still  another  frequency  at  which  perturbations  develop 
in  the  jet  which  are  comparable  to  perturbations  at  the  critical  fre¬ 
quency. 

We  note  that  significantly  greater  perturbations  develop  in  the 
cylindrical  jet  at  the  critical  frequency  than  in  the  plane  one; 
relationships  k,i»)  are  characterized  by  more  clearly  pronounced 
peaks. 

W.e  return  to  the  behavior  of  the  relationship  on  che  boundary 

of  the  jet  (Figures  8,9).  In  the  simplified  solution,  I<o  is  the 
amplitude  of  the  exterior  wave  (8)  when  y  *  1,  such  that 
The  relationships  in  the  complete  solution  are  depicted  by  the 

curves  with  a  large  number  of  maxima  and  minima.  In  this  case,  the 
values  of  the  extrema  increase  with  the  increase  in  the  frequency 
a,  and  a  particularly  sharp  increase  in  the  value  X„  is  observed  when 
a  >  1;  this  increase  is  the  more  intensive,  the  further  we  are  from 
the  nozzle. 

Figure  10  investigates  the  effect  of  the  Mach  number  of  the  jet 
M  and  the  ratio  -  of  velocities  of  sound  in  the  jet  and  in  the 
exterior  medium  on  the  distribution  on  the  boundary  of  the  jet 

for  critical  frequencies  in  the  cylindrical  case.  Value  Kp  does  net 
strongly  differ  from  the  amplitude  of  exterior  acoustical  waves  (i.e., 
from  1)  .  When  z  >  6,  the  contribution  of  secondary 

emission  from  the  surface  of  the  jet  to  total  emission  is  the  basic 
one,  and  the  value  of  K  increases  very  rapidly  with  the  increase  in  z, 
exceeding  the  amplitude  of  the  base  wave  5-10  times  when  z  *  8, 
and  10-20  times  when  z  *  9  (the  effect  of  the  boundary  condition  when 

*'».  *io  j_s  already  felt  in  the  data  when  z  *  9,  to  the  point). 

For  cold  jets,  the  intensive  increase  of  X  begins  earlier  relative 
to  z  (i.e.,  closer  to  the  nozzle)  with  a  decrease  in  M  number  of  the 
jet;  we  note  that  when  M  *  4,  such  an  increase  is  generally  not 
observed.  In  the  case  of  the  jet  with  different  temperatures, 

Ihk 


a  mere  raoid  increase  in  K  (z)  is  observed  for  the  hotter  jets, 

P 

although  the  differences  here  are  not  very  significant. 


The  value  K 
plane  jet. 


P 


behaves  similarly  as  well  on  the  surface  of  the 


The  behavior  of  K  within  the  jet  at  the  critical  frequency 

- 

is  characterized  by  the  fact  that  prior  to  *  \ 3  - * ,  the  relationship 
between  K  and  the  coordinate  z  is  close  to  monotonic  increase 
(Figure  5).  In  this  case,  the  complete  and  simplified  solutions  are 
similar;  the  relationships  nearly  do  not  differ  until  z  =  5, 

corresponding  to  various  .  When  z  >  6  in  regions  of  conti¬ 

nuity,  K  ,  one  has  sharp  drop  and  increases  in  K  .  Such  a  character 
of  the  behavior  of  Kp  is  clearly  defined  for  a  jet  with  small  M. 

Figure  11  plots  the  relationship  between  coordinate  z  and  the 
value  jt/*;  for  the  case  of  the  critical  frequency: 


where  -  amplitude  of  vibrations  of  the  jet  boundary.  A  comparison 
was  made  for  the  relationships  for  the  simplified  and  complete 

solutions.  These  two  solutions,  as  in  the  case  with  the  pressure 
distributions,  are  similar  when  z  <  5,  but  a  sharp  increase  in  the 
scale  of  vibrations  of  the  jet  occurs  for* the  complete  solution  at 
higher  values  of  z.  In  both  cases,  the  jet  vibrates,  significantly 
deforming  (bending).  There  are  points  on  the  boundary  of  the  jet 
which  remain  stationary  for  the  simplified  solution.  In  the  complete 
solution  these  points  correspond  to  points  of  the  minimum  range  of 
vibrations  of  the  stream  boundary.  We  note  that  the  value 
for  the  simplified  solution  with  the  given  «•**/»*•  and  *•*/.* 
does  not  practically  depend  upon  the  parameters  .  The  most 

rapid  increase  ’«,*  with  the  increase  in  z  is  observed  for  cold  Jets 
with  small  Mach  numbers  for  the  complete  solution. 

It  is  interesting  to  examine  a  perturbed  flow  in  the  exterior 
medium  due  to  secondary  emission  from  the  jet  surface.  Figure  12 


146 


plots  lines  of  equal  phase  (through  90°)  of  the  secondary  acoustical 
emission  (for  pressure)  in  the  plane  case  when  i,  s*  .  &*,  .  oje. 

It  is  evident  from  the  graph  that  pressure  perturbations  in  the 
exterior  medium  caused  by  pulsations  of  the  jet  in  it,  propagate 
approximately  like  cylindrical  waves  with  a  center  in  the  mid  point 
of  the  third  unit  of  the  jet. 


*  •  1,  M  *  3,  -  0.8,  ill  *  sfc.  -  0.28  . 

«.  - a*-  O  - »<f  -  180* 

•.  1.  Af  -  22*  ,5.  2.  *♦-«*.  3.  -  87*.  5. 

4.  »*  -  90* .  3.  A*  •  210* . 

The  picture  of  secondary  emission  in  the  case  of  the  cylindrical 
Jet  is  somewhat  more  complicated  (in  Figure  13,  for  M  *  3,  */«r  -  o.a 

(cold  jet),  »*•**.- ojb).  .  in  this  case,  we  once  again  observe 
the  presence  of  an  effective  source  of  emission  situated  on  the 
surface  of  the  Jet  within  the  confines  of  its  thi^d  unit  (when  z 
*  4.5).  In  the  vicinity  of  the  point  z  =  4.5,  the  equiphase  lines 
are  closed  lines  with  a  pole  at  this  point  (Figure  11B).  However, 
far  from  the  source  one  has  points  in  which  the  amplitude  of  the 
secondary  emission  vanishes.  Not  one,  but  several  sets  of 

phase  lines  run  through  such  points,  and  with  the  transition  through 
such  a  point  the  phase  changes  in  Jump  fashion  by  180°.  In  this  case, 
within  the  confines  of  the  calculated  region  (up  to  *«•  ),  one  does 

not  detect  a  tendency  toward  formation  of  closed  lines  with  a  common 
pole.  Thus,  unlike  the  plane  case,  perturbations  in  the  exterior 
medium  propagate  as  waves  from  a  certain  effective  source  only  in  a 
comparatively  small  vicinity  of  the  source,  and  subsequently  the  pic¬ 
ture  is  more  complex  in  the  far  acoustical  field. 


Qualitatively,  the  picture  of  phase  lines  is  precisely  the  same 
with  other''  .  «/«,  as  well.  The  position  of  the  effective  source 
depends  weakly  on  M;  for  heated  jets,  the  source  is  displaced  further 
from  the  nozzle  cross  section,  but  remains  within  the  confines  of  the 
third  cell  even  for  an  infinitely  heated  jet. 

The  Figure  14  plots  the  relationship  between  the  amplitude  of 
pressure  perturbations  in  the  exterior  medium  K^,  due  to  secondary 
emission  with  the  polar  angle  (the  pole  in  the  effective  source  on 

the  surface  of  the  jet,  i.e.,  when  * «  i,  ,0°  corresponds 

to  the  direction  from  the  source  to  the  nozzle)  at  distances  from 

the  source  j  *  0,5 ^  and  for  a  case  m  -  3,  «/«,  -  0.0, »*•»*.-  om. 


Figure  14 
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Key: 

(1)  pole  and  the  point 

(2)  direction  from  pole  to  nozzle 


As  is  evident  from  the  constructed  diagram,  secondary  emission  from 
the  surface  of  a  pulsating  jet  at  the  critical  frequency  has  the 
form  of  emission  from  a  dipole  with  an  axis  parallel  to  the  jet  axis, 


We  shall  interpret  the  obtained  results  according  to  secondary 
emission  and  its  effect  on  flow  in  the  Jet: 


At  the  critical  frequency,  secondary  emission  from  the  surface 
of  the  jet  or  at  least  near  this  surface  has  the  form  of  emission 
from  an  effective  source  situated  within  the  confines  of  the  third 
Jet  cell; 

at  basic  emission  frequencies  comparable  with  the  critical  fre¬ 
quency  and  exceeding  it,  the  amplitude  of  secondary  emission  can 
noticeably  exceed  the  amplitude  of  the  base  emission; 

a  sharp  Increase  in  the  amplitude  of  the  secondary  emission 
«;  on  the  surface  of  the  jet  .and  in  the  scale  of  vibrations  of  the 
boundary  of  the  jet  begins  beyond  the  third-fourth  cells  of  the  jet; 


it  also  leads  to  a  significant  change  in  the  character  of  per¬ 
turbed  flow  in  the  Jet; 

Tile  increase  in  the  value  «,  observed  at  jet  Mach  numbers  M 
<  3  is  particularly  sharp  for  cold  jets;  on  the  other  hand,  when 
secondary  emission  does  not  make  an  important  contribution  to  total 
emission;  for  Jets  with  identical  M  and  different  temperatures, 
differences  are  not  very  great; 

Ihe  conclusions  that  were  drawn  are  identically  valid  for  plane 
and  cylindrical  jets  and  in  the  cylindrical  case  the  noted  characteris¬ 
tics  of  secondary  emission  appear  more  sharply. 

It  is  interesting  to  compare  the  obtained  results  with  the  data 
of  an  experiment,  photography  of  plane  and  circular  jets,  which  study 
the  sound  of  the  discrete  frequency  cited  in  works  [8-14],  which 
make  it  possible  to  draw  .a  schematic  picture  of  jet  pulsations  (Figure 
15). 


Fiqure  15 


In  all  the  photographs  one  can  clearly  see  the  first  two  cells 
of  the  jet,  whose  displacement  relative  to  the  steady-state  position 
is  practically  absent.  This  is  in  agreement  with  our  data,  according 
to  which  the  amplitude  of  secondary  emission  in  small  in  this  region, 
as  are  displacements  of  the  jet  boundary.  The  two  succeeding  cells 
are  also  quite  clearly  evident  in  the  photographs;  however,  their 
vibrations  are  already  extremely  noticeable;  vortices  -  break  away 
from  the  surface  of  the  jet  into  the  surrounding  medium.  According 
to  our  own  results,  a  noticeable  Increase  in  the  amplitude  of  second¬ 
ary  emission  and  the  scale  of  jet  vibration  begins  within  limits  of 
the  third-fourth  cells.  This  increase  in  the  linear  solution  becomes 
exceedingly  rapid  behind  the  fourth  cell  of  the  jet.  In  the  full 
scale  model,  the  presence  of  non-linear  and  dissipated  mechanisms 
limits  the  increase  in  pulsations.  Actually,  it  is  evident,  according 
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to  the  photographs  that  a  powerful  formation  of  vortices  begins 
behind  the  fourth  cell  and  the  jet  practically  disintegrates. 

The  jet  cells  behind  the  fourth  cell  in  certain  photographs  are  in¬ 
distinguishable,  while  on  others  they  are  visible  with  great  diffi¬ 
culty — this  depends  upon  the  conditions  of  the  experiments,  specifically, 
on  the  presence  of  reflecting  surfaces  in  the  acoustical  field  of 
the  jet,  for  example,  the  end  of  the  nozzle  which  is  not  covered  with 
sound  absorbing  material. 

The  Figure  16  gives  a  qualitative  comaprison  of  our  results 
with  experimental  data  [14]  for  a  sonic,  super- expanded  jet.  In 
the  experiment,  the  pressure  .drop  in  the  nozzle  receiver  and  the 
surrounding  medium  was  o./p,  -3.07,  which  corresponds  to  an  average 
M  number  of  the  jet  of  ^  1.5.  The  points  in  Figure  16  plot  the  values 
of  the  amplitude  of  pressure  fluctuations  near  the  surface  of 

the  jet,  due  both  to  the  emission  of  noise  by  the  jet  and  to  the 
presence  of  pulsations  with  a  discrete  frequency,  taken  from  [14]. 

The  maxima  and  minima  in  the  distribution  of  are  in  fact  due  to 

the  presence  of  discrete  pulsations.  One  can  note  that  the  positions 
of  these  extrema  coincide  well  with  the  positions  of  the  extrema 
and  the  distribution  of  K  when  y  3  1,  plotted  according  to  the 

s' 

result  of  our  calculations,  for  a  cold  jet  with  M  *  1.5,  affected  by 
an  exterior  acoustical  wave  with  a  frequency  *'1  .  It  is 

evident  according  to  the  photograph 
[14],  that  when  z  >  8,  i.e.,  where  a 
rapid  Increase  in  Kp  takes  place  in 
the  linear  solution,  the  jet  has  totally 
disintegrated;  the  noise  spectrum  of 
the  jet  Is  even  in  this  region,  without 
discrete  components,  and  Ap  monoton- 
ically  diminishes. 


Thus,  there  is  a  qualitative 
coincidence  between  our  results  and  the 
experimental  data,  and  the  deterioration 
of  the  Jet  corresponds  to  a  rapid  growth  in  perturbations  on  the  surface 
of  the  Jet  in  the  linear  solution. 
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We  note  that  the  picture  of  acoustical  waves  emitted  by  jets 
can  be  reduced  to  a  picture  of  spherical  waves  emitted  by  a  certain 
effective  source t  with  good  accuracy.  The  position  of  the  source  can 
be  determined  according  to  photographic  data.  The  measurement  of 
photographs  made  in  a  number  of  works  shows  that  that  the  effective 
source  lies  within  the  region  of  the  sixth-eighth  cells  of  the  jet. 
It  is  interesting  that  one  other  source  of  emission  was  detected  in 
an  experiment  [10]  with  plane  jets  besides  the  basic  one.  The  addi¬ 
tional  source  is  located  at  the  beginning  of  the  third  jet  cell. 

This  source  can  evidently  be  compared  with  the  effective  source  of 
secondary  emission,  which  we  discussed  above  on  the  basis  of  the 
results  of  our  calculations. 


We  shall  now  turn  to  an  investigation  of  the  problem  of  inter¬ 
action  of  the  cylindrical  jet  with  a  symmetrical  acoustical  field 
(j  »  0).  Figure  17  plots  the  relationship  between  frequencies  Sh  of 
the  exterior  waves  of  the  amplitude  of  pressure  perturbations  K  for  * 
a  cross  section  z  -  6  when  y  *  0  (axis  of  the  jet)  and  y  *  1  (bound¬ 
ary  closed)  for  a  case  !  i  .  There  Is  a  maximum  within  the 

jet  Kq  at  a  frequency  equal  to  the  critical  frequency  [50],  When 
Y  3  x,  secondary  emission  does  not  make  a  noticeable  contribution  to 
total  emission. 


Figure  18  compares  the  amplitudes  of  vibrations  of  the  boundary 
of  the  jet  *r  in  the  symmetrical  and  antisymmetrical  cases  with 
critical  frequencies  of  the  corresponding  perturbations.  In  the 
symmetrical  case,  the  value  of  "J  increases  on  the  average  linearly 
with  the  increase  in  z  unlike  the  antisymmetrical  case,  where  the 
scale  of  vibration  of  the  Jet  sharply  builds  up  with  large  z. 
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The  cited  data  indicate  that  the  jet  in  the  symmetrical, 
exterior  acoustical  field  is  more  stable  than  in  the  antisymmetrical 
one.  This  result  is  in  accord  with  the  fact  that  no  jet  emission  of 
sound  at  a  discrete  frequency  with  a  symmetrical  character  of  the 
acoustical  field  is  observed  in  experiments. 

Figure  19  plots  phase  lines  of  the  secondary  emission  of  the 
jet,  confirmed  by  the  effect  of  the  symmetrical  pressure  wave  at  a 
frequency  equivalent  to  the  critical  frequency.  It  is  curious  that 
these  lines  are  parallel  straight  lines  throughout  nearly  the  entire 
calculated  region,  such  that  in  this  case  secondary  emission  is  a 
plane  wave  inclined  toward  the  jet  axis  at  an  angle  of  about^O0. 


3.  The  effect  of  a  subsonic  satellite  flow.  One  can  solve  the 
problem  of  the  interaction  of  the  supersonic  jet  with  the  exterior 
acoustical  field  with  the  presence  of  a  subsonic  satellite  flow  in 
the  exterior  medium  with  the  Mach  number  .  In  this  case,  the 

basic  equations  and  calculation  formulas  only  change  slightly.  Cer¬ 
tain  results  of  solving  this  problem  for  the  antisymmetrical  case  will 
be  given  below. 


By  using  the  simplified  statement  of  the  problem,  one  can  derive 
the  formula  for  critical  frequencies  of  exterior  acoustical  waves  in 
the  examined  case: 


K  (*) 


(53) 


where  the  function  h(n)  is  defined  in  (51). 


Figure  20  plots  the  relationship  between  the  critical  frequency 
s»,  and  the  number  Mg  of  the  satellite  flow  for  a  cold  cylindrical 
jet  with  M  »  2.5.  The  same  figure  plots  A.  G,  Yepifanskiy’ s  experi¬ 
mental  data  on  the  frequency  of  the  discrete  component  of  the  cold 
jet  in  the  satellite  flow.  In  the  experiment,  the  discrete  component 
was  registered  with  numbers  m,«  c.0-0.7  at  2  close  frequencies  (1  and 
11  in  Figure  20).  The  coincidence  of  our  results  with  the  experimental 
data  when  m^o.8-o.7  is  satisfactory  with  respect  to  the  frequency  value. 
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Figure  2C 
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-  critical  frequency  of  anti- 
symmetrical  ,  exterior  acoustical 
wave  for  the  cylindrical  case; 

- experimental  relationship 

of  frequencies  of  the  discrete 
components  in  the  noise  spectrum 
of  the  circular  jet. 
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Relationships  between  the  value  Kp  within  the  Jet  and  the 
ratio  -  of  the  frequency  of  exterior  waves  to  the  critical  fre¬ 
quency  (53)  for  n  ■  1  appear  in  Figure  21  for  several  values  for  M  . 

In  the  cylindrical  case,  formula  (53)  derived  when  solving  the  problem 
in  the  simplified  presentation,  determines  the  frequency  value  at  which 
the  Jet  maximally  intensifies  the  perturbations  transmitted  to  it,  and 
then,  when  one  takes  secondary  emission  into  account.  The  presence  of 
satellite  flow  slightly  alters  the  value  of  Kp  within  the  Jet.  In  the 
plane  case,  one  observes  displacement  of  the  basic  maximum  of  the  curve 
with  the  critical  frequency  to  a  frequency  -  •  cj?  .  with 


the  increase  in  the  Mg  number;  in  the  vicinity  of  the  latter  frequency  , 
the  value  of  Kp  weakly  depends  on 

Figure  22  gives  the  distribution  of  K  along  the  boundary  of  the 
jet  for  -a  -  3,  *;•,  -  i. m.-o ■  oj.  and  0.5.  All  data  pertain  to  cases 

of  critical  frequencies  of  exterior  waves.  With  the  presence  of  a 
satellite  flow  on  the  surface  of  the  jet  in  the  entire  calculated 
segment  (up  to  z  *  10),  secondary  emission  proves  to  be  slight  and  the 
amplitude  Kp  of  total  emission  is  nearly  1.  Thus,  the  presence  of  a 
satellite  flow  has  a  stabilizing  effect  on  the  jet.  These  data  are  in 
accord  with  the  results  of  an  experiment,  according  to  which  jets  in 
a  satellite,  subsonic  flow  do  not  emit  sound  of  a  discrete  frequency 
when  m,  »'  o,o  -  o.r 

Figure  23  plots  the  phase  lines  of  secondary  emission  from  the 
surface  of  a  cylindrical  jet  situated  in  a  satellite  flow  (m  -’3. m. •  oi. «/«,  - 1 
with  a  critical  frequency  of  the  exterior  wave  s*  •  *".  -o.ie.  .  it  is 
interesting  that  the  picture  of  the  phase  lines  in  this  case  does  not 
have  anything  in  common  with  the  corresponding  picture  for  the  case  of 
a  jet  flowing  into  a  flooded  space  With  the  presence  of  the 

satellite  flow,  perturbations  caused  in  the  exterior  medium  by  jet 
pulsations  propogate  like  plane  waves,  inclined  toward  the  jet  axis. 

This  picture  of  phase  lines  is  very  similar  to  the  picture  for  the 
case  of  symmetrical  Jet  pusations£^<  /£)  These  data  and  those  given  above 
indicate  that  both  the  presence  of  a  satellite  flow  during  anti  symmetrical 
perturbations  in  the  exterior  acoustical  field  and  "symmetrization"  of 
perturbations  have  an  identical,  stabilizing  effect  on  the  jet. 

4.  Basic  conclusions.  On  the  basis  of  the  results  of  the 
investigation  of  solutions  to  the  problem  of  the  interaction  of  a 
supersonic  Jet  with  the  exterior  acoustical  field  (in  the  linear  approx¬ 
imation  and  the  simplified  presentation  for  the  plane  case  in  [1,2]  >  in 
the  linear  approximation  and  the  complete  presentation  for  the  plane 
and  the  cylindrical  cases  in  this  work,  and  the  nonlinear  solution  in 
the  simplified  presentation  for  the  plane  case  in  [1,2])  one  can  draw 
the  following  conclusions: 


/.The  authors*  hypothesis  that  the  supersonic  jet,  because  of  ics 
periodic  (cellular)  structure  has  properties  of  a  unique  resonator  is 
confirmed:  During  interaction  with  the  exterior  acoustical  field,  rhe 
jet  intensifies  the  perturbations  transmitted  to  it  from  the  outside 
at  certain,  entirely  determined  frequencies.  In  accordance  with  the 
model  of  the  phenomenon,  accepted  in  the  work,  the  emission  of  sound 
at  a  discrete  frequency  by  supersonic  jets  should,  in  fact,  occur  at 
these  critical  frequencies.  This  is  confirmed  by  comparing  values  of 
critical  frequencies  calculated  according  to  the  formulas  cited  in  the 
work  with  the  experimental  values  of  frequencies  of  the  discrete  com¬ 
ponents  of  the  noise  spectrum  of  the  supersonic  jets. 

2.  During  the  interaction  of  a  supersonic  flowing  into  a  flooded 
space  with  the  antisymmetrical  acoustical  field,  perturbations  on  the 
jet's  surface  rapidly  increase  with  distance  from  the  nozzle  cross- 
section;  the  development  of  perturbations  is  particularly  rapid  behind 
the  fourth  jet  cell.  This  effect  does  not  exist  if  the  jet  flows  in  a 
subsonic  satellite  flow  with  a  sufficiently  high  number,  as  well  as 
during  the  interaction  of  the  jet  with  the  symmetrical  acoustical  field. 
In  all  cases  when  the  results  of  the  calculations  could  be  compared  with 
the  experimental  results,  it  proved  that  the  presence  of  emission  at  a 
discrete  frequency  in  the  experiment  corresponds  to  the  rapid  development 
of  perturbations  on  the  jet  surface  in  the  calculation. 

3.  The  investigation  of  the  nonlinear  solutions  shows  that  the 
jet  ceases  to  intensify  the  perturbations  arriving  from  the  outside  with 
sufficiently  high  amplitudes  of  the  exterior  acoustical  waves.  This 
Indicates  the  possible  mechanism  of  limiting  the  growth  of  pulsations 

in  the  Jet,  which  leads  to  the  establishment  of  self-excited  vibrations. 
It  is  in  fact  the  jet  which  intensifies  the  perturbations  transmitted 
to  it  by  the  exterior  wave.  In  this  case,  perturbations  of  already  high 
intensity  are  emitted  in  the  region  of  its  disintegration,  which  will 
be  intensified  in  their  turn  by  the  jet,  etc.  However,  when  sufficiently 
high  perturbations  develop  in  the  jet,  this  leads  to  the  distortion  of 
its  periodic  structure  and  deterioration  of  its  resonator  properties, 
as  a  result  of  which  the  increase  in  perturbations  ceases. 
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J'  In  orde1,  t0  oul">lnate  the  study  of  the  examined  problem 
la  the  Unear  approximation.  It  Is  necessary  to  Investigate  the 
effect  of  the  fora  of  basic  emission  *  as  well  as  the  waves' ~fieo-ed 
by  the  various  screens. 


The  linear  solutions  make  It  possible  to  obtain  the  statistical 
characteristics  of  the  acoustical  field,  which  take  into  account 
the  random  character  of  jet  pulsations  [1,2].  Great  Interest  also 
a  .acnes  to  the  further  study  of  this  problem  in  the  linear  oresent- 
atlon,  specifically  the  mechanism  of  acoustical  wave  emission  from 
the  region  of  jet  disintegration. 
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STABILITY  OF  AN  AXISYMMETEICAL 
COMPRESSED,  N0N-VI3C0US  WAKE 


by 

S.  Ya.  Gertsenshteyn  and  A.  V.  Kashko 

The  article  gives  certain  results  of  calculating  the  stability 
of  flow  in  the  wake  behind  a  body  moving  (flying)  at  hypersonic 
velocity.  Wave  numbers,  rates  of  propagation  and  amplification 
coefficients  were  obtained  which  are  characteristic  for  the  vibra¬ 
tions  which  arise  in  the  wake.  The  relationships  between  these 
values  and  the  Mach  number  and  the  temperature  drop  on  the  axis  of 
the  wake  and  the  periphery  are  examined. 

The  article  also  gives  a  brief  survey  of  theoretical  investiga¬ 
tions  on  stability  of  flow  in  the  far  wake. 

1.  It  is  well  known  from  many  experimental  works  [1,2]  that 
sinusoidal  oscillations  appear  in  the  supersonic  wake  behind  the 
"throat"  (at  a  distance  of  10-15  diameters  of  the  body).  These  oscill 
ations  develop  quite  rapidly  and  at  a  distance  of  approximately  25 
diameters  downstream  from  the  "throat"  of  the  wake,  flow  becomes  tur¬ 
bulent.  A  study  of  the  given  oscillations  is  of  certain  interest 
and  a  number  of  theoretical  works  have  been  devoted  to  .this  question. 

Basically,  results  were  obtained  for  plane  wakes.  The  purpose 
of  this  work  is  to  study  axisymr.etrical  compressible  wakes. 

We  note  that  the  analogy  between  plane  and  axisymmetrical  wakes 
is  far  from  complete.  For  example,  it  has  been  shown  in  the  axisymmet 
cal  case  for  quite  "smooth"  profiles  [11],  that  2- dimensional, 
increasing  oscillations  do  not  exist  and  only  a  3- dimensional ,  sinu¬ 
soidal  wave  can  develop  in  the  flow.  In  the  plane  case,  however, 
according  to  the  Squire  theorem,  plane  perturbations  are  the  most 
unstable  ones.  As  tne  conducted  calculations  show,  the  corresponding 
numerical  data  for  axisymmetrical  and  plane  wakes  also  significantly 
differ. 
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Evidently,  Helmholtz  first  explained  the  appearance  of 
oscillations  in  the  wakes  with  instability  at  the  vortex  layers 
which  form  in  the  wakes.  Raley  [3]  examined  and  solved  the  prcb= 
lem  of  the  stability  of  the  vortex  sheet.  He  also  solved  the  prob¬ 
lem  of  stability  of  the  vortex  layer  of  finite  thickness.  The 
problem  of  flow  stability  in  a  wake  with  vortex  distribution  of 
finite  intensity  was  first  examined  by  G.  I.  Petrov  [4],  In  this 
work  the  velocity  profile  has  been  replaced  by  a  broken  line  con¬ 
sisting  of  5  components  (J+i-  i)\  the  solution  was  obtained  by  analogy 
with  the  above-mentioned  Rayleigh  solution.  With  the  aid  of  com¬ 
paratively  simple  calculations,  it  was  specifically  demonstrated 
in  a  work  [4]  that  the  antisymmetrical  perturbations  developed  more 
intensively  than  the  symmetrical  ones.  In  this  work  [4]  with  the 
aid  of  the  generalized  Rayleigh  method,  an  examination  was  made  of 
stability  of  an  axisymmetrical  wake  relative  to  zonal  oscillations. 
Certain  results  of  a  work  [4]  were  duplicated  and  refined  In  the  work 
of  Michalke  and  Schade  [12].  The  work  [12]  also  employed  the  Rayleign 
method,  but  the  velocity  profile  was  replaced  by  a  broken  line  with 
a  significantly  larger  number  of  components  than  in  the  work  [4]. 

We  note  that  the  overwhelming  majority  of  works  on  wake  stability 
relate  to  the  plane  wake.  Basically,  three  works  [11,  4,  13]  contain 
results  on  the  stability  of  axisymmetrical  wakes.  The  stability  of 
flow  relative  to  the  zonal  oscillations  was  examined  in  a  work  [4], 
as  has  already  been  mentioned  above.  A  work  [11]  derives  certain 
estimates  and  criteria  which  are  necessary  and  sufficient  for  the 
existence  of  neutral  and' increasing  oscillations,  by  analogy  with 
the  known  criteria  for  the  plane  case.  Greatest  interest  attaches 
to  the  necessary  condition  for  the  existence  of  increasing  perturba¬ 
tions,  obtained  in  a  work  [11].  It  will  be  briefly  explained  below 
how  this  condition  can  be  obtained  and  what  conclusions  can  easily 
be  drawn  from  it. 

The  problem  of  the  stability  of  an  axisymmetrical  wake  in  an 
ideal,  incompressible  fluid  can  pertain  to  a  problem  for  eigenvalues  ■ 
for  a  second  order  equation 
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a  l  T  ■ 

dz 


j-il’-c'G-rS—  l- 


(1) 


Here  the  boundary  conditions  -  conditions  of  solution  limitation 
when  x  *  0  and  t  -*•  ®.  &(*>  -  basic  flow  velocity  profile,  ct  and  n  - 

assigned  parameters  (perturbation  was  presented  in  the  following 
form) : 


$'■  5 (t)  .  l*ct 

By  multiplying  equation  (1)  by  ^rc  (here  e  is  a  function  con¬ 
jugate  of  G)  and  integrating  the  obtained  expression  from  t  =  0  to 
t  -*•  00 ,  having  used  the  known  boundary  conditions,  integrating 

by  terms,  and  by  separating  the  real  and  imaginary  part, 
one  easily  can  obtained: 

c.  //*/  tf'3'*  x  0  ( 2 ) 


here 


If  c. •  £?,  ,  then  it  follows  from  (2)  that  Q’  can  change  sign  in  the 

interval  i ,  i.e.,  when  c„*a  should  vanish  at  a  certain 

point  of  the  interval  (0,°°).  For  example,  it  is  easy  to  show  from 
this  necessary  condition  that  9'*o  for  “>.)*  when  n  *  C, 


It  was  shown  in  work  [11]  that  increasing  perturbations  can  only 
exist  for  this  profile  when  n  3  1,  i.e.,  only  the  "sinusoidal" 
oscillations  can  be  unstable. 


It  was  noted  in  a  work  [11]  that  this  conclusion  can  be  related 
solely  to  the  "smoothed"  velocity  profile. 

The  given  conclusion  is  not  true  for  profiles  with  a 

sharp  velocity  change.  No  specific  calculations  relative  to  increas¬ 
ing  perturbations  are  given  in  work  [11]. 

Work  [13]  is  chiefly  devoted  to  plane  wakes  in  a  compressible, 
non-viscous  fluid.  Only  individual  calculations  are  given  for  the 
axlsymmetrical  case  in  the  work  for  a  Gaussian  distribution  of  the 
velocity  profile  when  n  ■  1,  i.e.,  for  the  "sinusoidal" 

waves.  A  work  [13]  not  only  gives  specific  calculations,  but  also 
derives  certain  necessary  criteria  for  the  existence  of  neutral 
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perturbations  and  a  certain  sufficient  condition  for  the  existence 
of  increasing  perturbations.  The  conducted  investigation  is 
similar  to  investigations  made  in  the  case  of  the  incompressible 
fluid  [5,6],  The  authors  also  used  the  investigations  of  Lynn  and 
L.  Liz,  which  are  unpublished  in  the  Soviet  Union,  and  which  deal 
with  the  stability  of  a  compressible  boundary  layer.  In  addition  to 
the  work  referenced  above  [13],  stability  of  the  compressible  plane 
wake  is  investigated  in  works  [7, 14,15],  Works  [7,14]  investigate 
the  plane  case.  A  work  [15]  examines  the  plane  and  axisymmetrical 
cases.  The  basic  flow  in  this  case  is  assumed  to  be  constant  in  the 
vicinity  of  the  axis,  i.e.,  stability  of  an  axisymmetrical  vortex 
sheet  in  a  compressible  gas  was  examined. 

2.  In  order  to  investigate  oscillations  in  an  axisymmetrical 
wake,  one  can  take  the  original  equation  system  in  the  cylindrical 
coordinate  system  [8]: 
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4T  p  dp 
P  ■ pKT . 


(3) 


We  shall  present  the  sought  solution  to  the  system  in  the  form  of 
the  sum  of  the  basic  flow 
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and  a  certain,  sufficiently  small  perturbation: 
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T  .  T  •  T' 


Having  substituted  (4)  in  (3),  and  ignoring  the  terms  which 
are  quadratic  relative  to  the  perturbation  components,  we  obtain 
a  system  consisting  of  6  linear  equations  relative  to  6  unknown 
functions : 
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(5) 


If  all  perturbation  components  are  sought  in  the  form, 

then  it  is  easy  to  obtain  a  system  of  ordinary  differential  equations 
from  the  system  (5)  for  determining  the  function  {<*’  : 


v»*#  ♦  v’«  **"**  *  rs 

1  '■  ft)  ^ 

*  ...  .  . 
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We  shall  Insert  the  dimensionless  variables:  •'  . 

„  p  }  l  » 

<  •  r  ?  o.  i  i  -  ?  - 

.»  1  •  ■  p  '  •  j,  •  *  ^ 

System  (6),  in  the  dimensionless  variables  has  the  form: 
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Here  M  is  the  Mach  number:  «*•-«  ■ 

T.K  '• 

It  is  subsequently  convenient  to  reduce  the  system  (8)  tc  1  equa¬ 
tion.  We  shall  write  the  equation  for  pressure  4 .  It  is  prelim¬ 
inarily  expedient  to  obtain  the  expression  for  through  y  and 
.  For  this  purpose,  it  is  sufficient  to  substitute  the  expression 
for  $,  from  equation  (8)  in  equation  (10).  Having  substituted  the 
expression  for  t  from  equation  (13)  in  equation  (12),  we  then 
obtain  the  relationship  between  p ,  r, ,  +,  ,  obtained  from 

equation  (11),  having  substituted  the  expression 

equation  (13)  in  equation  (12),  instead  of  ?  its  expression  through 
V#,  and  .  If  one  now  uses  the  expressions  already  found  for 

,^,and  via  ?  and  *'*  ,  then  we  obtain  the  following  equation 

for  f  : 


r 


(14) 
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Here 


3.  Boundary  conditions  for  Equation  (14)  may  be  respresented 
simply:  at  z  =  0  p(0)  =  0  and  at  «-•  the  solution  *C*) 

must  be  bounded  in  absolute  value.  However,  close  to  the  point 

*  -o  it  is  advantageous  to  integrate  Equation  (14)  numerically, 
since  expressions  of  the  form  ■  ~  and  d*  are 

included  in  the  coefficients  of  the  equation.  Therefore  it  makes 
sense  to  use  the  asymptotic  representation  of  the  solution  in 
the  vicinity  of  zero. 

To  perform  a  numerical  calculation,  it  is  necessary  also  to 
replace  the  boundary  condition  at  infinity  by  the  boundary 
condition  at  a  finite  distance.  This  may  be  done  by  using  the 
asymptotic  expansions.  Below  we  give  the  corresponding 
calculations . 


If  at 


» » 1  vt «  1 


,  the  Equation 


(14)  assumes  the  form 
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The  asymptotic  behaviour  at  infinity  of  the  bounded  solution  of 
Equation  (16)  is  well  known  [9,10]: 


where 


and 


,  *  .it 

x  •  /-  ♦/)  '  “ 

<#  x  >  0 


(13) 


We  shall  not  consider  the  case  ?#i  -o  »  since  at 

*  0  the  boundary  condition  at  infinity  is 

automatically  satisfied  for  any  8  (case  of  a  continuous  spectrum) 
Using  (17),  we  may  readily  see  that  for  very  large  values  of  x 
the  following  relations  must  be  satisifed  with  great  accuracy: 


/>.’  1  p^P-.Q 


(19) 


This  relationship  (19)  will  replace  the  boundary  condition  at 
infinity  below. 


Let  us  write  the  asymptotic  solution  of  Equation  (14)  in 
the  vicinity  of  zero.  Equation  (14)  may  be  rewritten  in  the 
following  form: 


^1'  ■ 


(20) 


Let  us  give  the  solution  of  Equation  (20)  in  the  form  of  the  series 


p  .  r  c 


Collecting  terms  of  the  form  * *"*.  ,  we  have 


\(\- f)c,-  j'c. 


from  which  we  have 


The  remaining  coefficients  of  the  series  may  be  readily 


determined  by  means  of  C  .  For  example,  collecting  terms  of  the 
form  ,J'\  ,  we  obtain  : 


./So,)’  a",  7 

C,  X  •)-€,*  T'C,*iC'[  IjtI.J 


Thus,  in  the  vicinity  of  zero  the  unknown  solution  behaves 
like  *r  .  It  is  assumed  that  t*o.  .  For  Tm° 
we  may  obtain  the  corresponding  asymptotic  representation  by 
analogy  with  the  previous  one  [27]. 
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4.  These  calculations  make  it  possible  to  obtain  the 
representation  of  the  nature  of  vibrations  in  a  supersonic 
medium.  As  an  example,  we  examined  flow  in  the  medium  behind 
a  body  in  the  range  of  Mach  numbers  from  10  to  30.  The  flow 
instability  is  considered  in  the  wake  behind  the  "throat", 
at  a  distance  10  -  15  times  greater  than  the  body  dimensions. 

Data  on  the  velocity  and  temperature  distribution  were  taken 
from  [1,2,13]: 

f  4 

u  •  t  •  *  «4  e 

f  •  1  ♦  A  f  *  * 

; 

t  •  !  fry. 

i  'n, 

* 

Here  is  the  dimensionless  velocity;  *  ,>  is  the  tem¬ 
perature  distribution;  -  enthalpy  in  the  cross  section 

we  selected  on  the  wake  axis;  *«  -  enthalpy  in  the  advancing 

stream. 
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Figure  1 
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Figure  2 
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Figure  3 
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Figure  4 


The  results  of  the  calculations  are  partially  given  on 
graphs  1-6.  As  shown  in  Figures  1,2,3,^,  with  an  increase 
in  the  Mach  number  the  amplification  coefficient  decreases, 
but  very  little.  There  the  well  known  sharp  increase  in 
stability  of  flow  with  an  increase  in  the  Mach  number  can  be 
better  explained  by  the  large  phase  propagation  rate  of 
perturbations  in  the  compressed  wake  (Figure  5)  -  (the  per¬ 
turbations  very  rapidly  drift  along  the  stream,  without 
becoming  developed).  It  is  of  interest  to  compare  the  results 
obtained  with  the  well  known  results  [13]  on  the  stability  of 
plane  streams.  The  comparison  shows  that  the  amplif ication 
coefficients  for  plane  streams  are  approximately  four  tines 
greater  than  for  axisymmetric  streams. 


These  calculations  also  show  that  with  an  increase  in  the 
Mach  number,  the  extremal  wave  number  a*,  which  corresponds 
to  the  maximum  amplification  coefficient,  decreases,  and 
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Figure  6 


the  phase  velocity  « *  remains  almost  unchanged.  We  should 
also  note  that  the  phase  velocity  changes  little  with  a  change  in 
the  wave  number  a  (for  fixed  Mach  number  )  -  Figure  5. 

Separate  calculations  were  made  to  clarify  the  effect  of 
a  hot  "core"  for  an  axisymmetric  wake. 

With  an  increase  in  AT  (from  0.1  to  2.0),  the  amplification 
coefficient  increases  by  approximately  a  factor  of 

1.5  (Figure  6).  The  extremal  wave  number  «'  ,  with  an 

increase  in  AT  from  0.1  to  2.0,  decreases  by  approximately  1/10 
of  the  initial  value.  The  perturbation  propagation  rate  in  the 
range  examined  changes  very  little.  We  should  note  that  the 
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for  the  observer,  who  moves  at  the  perturbation  propagation 

rate  c,  ,  ,  the  local  flow  velocity,  whose  stability 

is  being  studied,  is  less  than  the  speed  of  sound.  This  occurs 
in  the  entire  range  in  which  AT  and  M  change. 
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